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ABSTRACT 

This  thesis  evaluates  several  bias  minimizing  and 
variance  minimizing  experimental  designs  in  terms  of 
their  effectiveness  and  efficiency  in  constructing 
response  equations  for  a  deterministic,  nuclear  exchnage 
problem.  The  criteria  which  is  used  to  evaluate  these 
designs  includes:  1)  number  of  required  design  points;  2) 
number  of  terms  in  the  response  equation;  3)  accuracy  of 
fit  of  the  response  equation;  4)  orthogonality  of  the 
design;  and  5)  rotatability  of  the  design. 

In  addition,  the  response  surface  equations  are 
evaluated  in  terms  of  their  predictive  power  and  their 
explanatory  power.  The  predictive  power  addresses  the 
equation's  ability  to  adequately  estimate  the  true 
surface  and  to  accurately  predict  a  future  response  for  a 
given  set  of  inputs.  The  explanatory  power  addresses  the 
equation's  ability  to  present  a  response  equation  which 
is  simple  to  interpret  so  that  the  true  surface  can  be 
easily  evaluated  and  the  results  can  be  easily  explained. 


I. OVERVIEW 


BACKGROUND 

Response  Surface  Methodology  (RSM)  is  a  method  for 
-fitting  mathematical  models  to  surfaces  generated  by 
experiment.  RSM  is  concerned  with  the  functional 
relationship 

V  =  h(U,,  Uj,.  .  .  ,  U,^)  (1) 

between  the  response  variable,  V,  and  the  n  independent 
variables  Ui,  U2  The  response  relationship  is 

usually  expressed  as  a  low  order  polynomial  or  some  other 
non-linear  function  which  approximates  the  true  surface. 
Graphically,  the  response  relationship  is  represented  as  a 
surface  in  a  space  whose  coordinates  are  the  n+l  variables, 
U|,  Uj  1---1  Un  and  V.  For  example,  the  response  surface 
can  be  expressed  as  a  series  of  contours,  similar  to  the 
rise  and  fall  of  land  masses  on  a  topographic  map. 

In  applying  RSM,  the  analyst  selects  an  experimental 
design  <i.e.,  a  theory  which  requires  the  least  number  of 
points  to  accurately  fit  the  model)  such  that  the  functional 
rel ati onship 

y  =  f(x^,  X2,>->i  x^)  (2) 

fitted  by  the  method  of  least  squares  closely  approximates 
the  true  function 

y  —  Xj,...,  Xp)  (3) 

over  a  specified  region.  When  f  inadequately  represents 


g,  there  are  two  errors  associated  with  this  problem 


1)  variance  error  -  which  is  error  due  to  sampling; 
and  2)  bias  error  -  which  is  error  due  to  the  lack  o-f 

fit  of  f  to  g 

Initially,  criteria  for  evaluating  the  adequacy  of 
designs  have  focused  primarily  on  minimizing  the  variance 
error  while  the  question  of  bias  error  has  been  given 
somewhat  secondary  considerations.  However,  a  bias 
minimizing  design  takes  on  added  significance  when  RSli  is 
applied  to  a  deterministic  model  because  there  is  no 
sampling  error  associated  with  this  model.  In  this  context, 
a  deterministic  model  is  a  model  in  which  the  response 
variables  are  free  from  stochastic  variation. 

For  example,  non-determini Stic  models,  such  as 
computer  simulations,  yield  different  results  (due  to  random 
sampling  of  distributions)  each  time  the  test  is  conducted. 
However,  deterministic  models,  such  as  linear  programming 
problems,  yield  the  same  results  each  time  they  are  applied 
because  sampling  is  not  involved  for  the  fixed  output  level. 
Therefore,  any  error  associated  with  the  deterministic  model 
is  solely  attributable  to  bias  error  (lack  of  fit)  and  none 
of  it  to  variance  error  (sampling). 

In  addressing  the  concept  of  experimental  designs  for 
response  surfaces,  most  of  the  RSH  literature  has  focused  on 
those  designs  which  minimize  variance.  These  designs  have 
received  a  large  share  of  attention  because  RSH  has  been 
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almost  exclusively  applied  to  non-determini stic  models. 
Examples  of  variance  minimizing  designs  are  provided  in  Box 
and  Behnken's  work  <1:455)  on  three  level  factorial  designs. 

Thus  far,  only  a  limited  amount  of  literature  has  been 
published  on  bias  minimizing  designs.  Only  the  basic  theory 
of  bias  minimizing  designs  exists  in  the  literature  with  the 
applications  of  the  theory  still  to  be  fully  documented.  In 
general,  the  applications  of  RSli  have  focused  on  generating 
response  surface  equations  for  two  purposes:  to  predict 
future  responses  and  to  determine  the  input  values  of  the 
optimum  response.  In  determining  the  optimum  operating 
conditions,  the  experimenter  employs  various  search 
techniques,  such  as  gradient  search  methods,  to  solve  the 
problem.  However,  in  examining  the  response  surface  of  a 
deterministic  model,  the  experimenter  may  not  be  interested 
in  determining  the  optimum  operating  conditions.  Rather,  as 
in  this  thesis,  the  experimenter  may  be  interested  in 
adequately  representing  the  true  surface  with  his  predicted 
surface  so  that  he  can  explore  the  relationships  between  the 
input  variables  and  accurately  predict  the  response  values. 
As  a  result,  the  experimenter  is  interested  in  selecting  a 
design  which  allows  him  to  employ  simple  analytical 
techniques  to  the  response  equation  to  examine  these 
relationships. 

The  basic  theory  of  bias  minimizing  designs  proposes 
that  a  design  can  be  constructed  which  requires  fewer  number 
of  design  points  and  which  achieves  a  better  fit  of  the 


function  to  the  response  surface  than  a  variance  minimizing 
design  because  a  bias  minimizing  design  does  not  consider 
any  variance  in  its  formulation.  Manacapilli  (17:6.1)  has 
stated  that  for  experiments  with  four  or  more  variables,  a 
bias  minimizing  design  provides  a  function  with  a  better  fit 
to  a  surface  than  a  variance  minimizing  design.  However, 
the  bias  minimizing  design  requires  more  design  points  to 
achieve  the  fit.  Thus,  these  results  indicate  that  not  all 
bias  minimizing  designs  achieve  a  better  fit  with  fewer 
design  points. 

problem  ilAIEMENI 

In  developing  experimental  designs  for  response 
surfaces,  researchers  have  focused  on  variance  minimizing 
designs  as  the  basic  criterion  in  order  to  capture  the 
response  surface.  However,  deterministic  models  possess 
characteristics  whereby  a  bias  minimizing  design  may  be  a 
more  appropriate  choice.  Thus,  the  problem  is  to  find  a 
bias  minimizing  design  which  best  fits  the  response  surface 
of  a  deterministic  model  and  to  compare  the  efficiency  and 
effectiveness  of  this  design  to  a  variance  minimizing 
design. 

BiSiARCH  QUESIION 

What  experimental  designs  produce  effective  and 
efficient  response  surfaces  of  deterministic  models? 
SUBSIDIARY  QUESIIONS 

RSM  when  applied  to  a  deterministic  model  is 


investigated  to  find  bias  minimizing  experimental  designs. 
The  results  of  this  design  are  evaluated  over  various  factor 
levels  and  its  effectiveness  is  compared  to  the  results  of 
variance  minimizing  designs.  The  primary  measures  of 
effectiveness  (MOE)  are:  orthogonality  of  the  design, 
accuracy  of  fit  to  the  true  response  surface,  and  number  of 
required  design  points. 

Some  of  the  questions  to  be  answered  from  these  MOEs 

are: 

1)  Mhat  are  the  tradeoffs  of  using  a  bias  minimizing  design 
instead  of  a  variance  minimizing  design  to  fit  the  response 
surface  of  a  deterministic  model? 

2)  When  is  this  design  more  efficient  or  more  effective 
than  a  variance  minimizing  design? 

3)  Is  the  analysis  of  response  surfaces  made  easier  with  a 
bias  minimizing  design? 

4)  Is  the  bias  minimizing  design  a  practical  method  (in 
terms  of  the  required  number  of  design  points)  to  fit  a 
response  surface  of  a  deterministic  model? 

5)  Do  these  designs  produce  equivalent  results? 

6)  Can  an  orthogonal  bias  minimizing  design  be  constructed 


to  facilitate  the  analysis  of  the  surface? 


LiliBAiyOi  Biyiiw 

The  concepts  of  Response  Surface  Methodology  were  first 
developed  by  Box  and  Wilson  in  the  early  195E)'s  (5:1).  This 
methodology  is  designed  to  be  a  practical  and  efficient 
manner  in  dealing  with  the  problem  of  determining  the 
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optimum  operating  conditions  -for  a  specific  process.  For 
example,  RSM  is  ideally  suited  for  agricultural  and  chemical 
experiments  in  which  the  experimenter  seeks  to  maximize  the 
yield  or  purity  of  the  product.  An  extensive  list  of  the 
applications  of  RSM  in  chemical,  agricultural,  and  other 
related  fields  can  be  found  in  Hill  and  Hunter's  article 
(11:591) . 

More  recent  applications  of  RSM  have  been  in  the  area 
of  animal  husbandry.  Articles  by  Toyomizu,  Akiba, 

Horiguchi,  and  Matsumoto,  (24:886)  and  by  Roush,  Petersen, 
and  Arscott,  (21:1504)  have  identified  the  uses  of  RSM  in 
determining  the  optimum  operating  conditions  for  raising 
chickens. 

With  regards  to  military  applications.  Smith  and 
Mellichamp  (22)  used  RSM  to  perform  multidimensional 
parametric  analysis  to  study  mathematical  programming 
models.  In  addition,  several  recent  Air  Force  Institute  of 
Technology  (AFIT)  thesis  efforts  have  applied  Smith  and 
Mellichamp 's  approach  and  techniques  to  various  military 
related  problems.  For  example,  Manacapilli  (17)  applied  RSM 
along  with  economic  production  functions  to  a  nuclear 
exchange  linear  programming  model.  Manacapilli  showed  that 
economic  production  functions  could  be  fitted  to  the 
response  surface  of  a  deterministic  model  so  that  basic 
economic  theory  could  be  used  to  analyze  the  surface. 

Another  application  of  RSM  to  a  nuclear  exchange  model 
is  the  work  by  Graney  (19).  In  his  thesis,  Graney  combined 


two  different  response  surfaces  (representing  two  different 
objectives)  as  a  means  of  evaluating  a  multiple  objective 
problem. 

Further  use  of  Smith  and  lielli champ's  paper  in  this 
research  effort  is  defined  in  Chapter  II  under  the  proposed 
methodology  section. 

iiNiRBL  meihodology 

In  general,  the  basic  approach  in  answering  the 
research  question  is  to  apply  various  experimental  designs 
to  a  response  surface  of  a  deterministic  model  and  to 
evaluate  the  differences  in  the  results. 

Specifically,  the  first  step  is  to  identify  bias 
minimizing  designs,  such  as  the  Box-Draper  designs  and  the 
Koshal  designs. 

The  second  step  is  to  identify  variance  minimizing 
designs.  Several  variance  minimizing  designs  are  considered 
in  this  study.  These  designs  include,  but  are  not  limited 
to,  second  order  rotatable  designs  and  central  composite 
rotatable  designs.  These  types  of  designs  are  selected 
because  of  their  frequent  use  in  RSM. 

Next,  both  the  bias  minimizing  and  variance  minimizing 
designs  are  applied  to  a  deterministic  model.  The 
deterministic  model  to  be  used  in  this  research  effort  is  a 
nuclear  exchange  model.  The  nuclear  exchange  model  is  a 
simple  linear  programming  model  designed  to  optimally 
allocate  the  number  of  strategic  weapons  per  weapon  type 


(input  variables)  to  targets  so  that  maximum  damage  occurs 
(response  variable).  This  model  was  selected  because  o-f  its 
basic  simplicity  and  because  past  AFIT  thesis  e-fforts,  by 
Graney  (19:51)  and  by  Manacapilli  (17:6.1)  have  used  a 
similar  model  with  -Favorable  results. 

Finally,  the  Functions  obtained  from  the  two  types  of 
experimental  designs  are  evaluated  in  terms  of  their 
effectiveness  and  efficiency. 

In  terms  of  effectiveness,  the  fit  of  the  postulated 
surfaces  to  the  actual  surface  is  measured.  To  measure 
the  fit  of  a  postulated  surface,  a  random  set  of  data  points 
from  the  actual  surface  is  collected  and  is  applied  to  the 
function  which  is  used  to  generate  the  postulated  surface. 
The  differences  between  the  expected  response  and  the  actual 
response  are  then  computed.  The  sum  of  the  deviations 
squared  (SDS)  is  then  obtained  and  a  comparison  is  made 
between  each  of  the  functions  estimated  by  different  designs 
based  on  their  SDS  value.  The  designs  which  generate  models 
which  most  closely  represent  the  true  surface  are  the 
model  which  minimizes  SDS. 

In  terms  of  efficiency,  the  minimum  number  of  design 
points  required  to  generate  the  surface  are  evaluated.  As 
stated  earlier,  Manacapilli  showed  that  a  bias  minimizing 
design  required  more  design  points  than  a  variance 
minimizing  design  while  gaining  only  a  small  improvement  in 
the  fit  of  the  postulated  surface.  If  this  result  holds  for 
the  generalized  bias  minimizing  design,  then  the 
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ef -Feet  1  ven ess  and  ef-Ficiency  trade-o-Ffs  between  bias  and 
variance  minimizing  designs  are  evaluated. 

One  o-F  the  measures  of  this  trade-off  is  an 
effectiveness-efficiency  factor.  This  factor  is  given  by 

E=<error  in  fit  of  design) * (number  of  design  points)  (4) 

where  the  error  in  the  fit  of  the  design  is  the  mean  of  the 
differences  between  the  predicted  and  the  actual  responses 
divided  by  the  actual  responses. 

This  factor  simultaneously  takes  into  consideration 
the  precision  with  which  the  design  estimates  the  true 
surface  and  the  number  of  required  experimental  design 
points  (18:135). 

The  next  chapter  expands  the  methodology  and  discusses 
the  specific  concepts  involved  in  applying  RSM  to  a 
deterministic  model. 


II.  METHODOLOGY 


INIRODUQIION 

This  chapter  details  the  methodology  used  in  applying 
Response  Surface  Methodology  (RSM)  to  a  deterministic  model. 
It  builds  on  the  general  methodology  described  in  Chapter  I 
by  explaining  the  various  experimental  designs  which  are 
used  in  this  study,  by  detailing  the  nuclear  exchange 
linear  programming  (LP)  models  and  by  outlining  the 
application  of  experimental  designs  in  generating  the 
response  surfaces. 

SECOND  ORDER  EXPERIMENIAL  DESIGNS  FOR  RESPONSE  SURFACES 

A  major  goal  of  Response  Surface  Methodology  is  to 
answer  the  question  of  what  happens  to  the  response 
variable  when  input  variables  are  allowed  to  vary  over  a 
specified  range.  Instead  of  utilizing  a  one  variable  at  a 
time  procedure  (which  changes  the  level  of  one  factor  while 
fixing  all  other  input  variables  at  a  specified  level),  the 
theory  of  experimental  design  is  utilized.  By  applying  the 
appropriate  experimental  designs,  a  response  function  is 
generated  to  estimate  the  response  value  for  any  given 
combination  of  input  variables.  Two  books  which  are 
excellent  references  on  the  theory  of  experimental  design 
have  been  written  by  Davies  (8)  and  by  Hicks  (113).  An 
extensive  list  of  applications  of  experimental  design  in  RSM 
can  be  found  in  Steinberg  and  Hunter's  article  (23:71). 

Furthermore,  a  significant  portion  of  RSM  literature 


has  -focused  on  second  order  experimental  designs.  That  is, 
designs  specifically  constructed  for  fitting  a  second  order 
polynomial  equation  to  the  data.  This  type  of  design 
produces  a  second  order  polynomial  equation  of  the  form  : 

y  =  bo+EbiXi  +  Ebii  xf  +EbijX,Xj  (5) 

This  equation  includes  the  first  few  terms  of  a  Taylor 
series  expansion.  Therefore,  a  second  order  design  usually 
produces  an  adequate  response  equation  in  estimating 
the  actual  surface  because  the  response  equation  depicting 
the  actual  surface  can  usually  be  represented  by  a  Taylor 
series  expansion.  The  difference  between  the  actual  and 
predicted  surfaces  is  due  to  the  higher  order  terms  of  the 
Taylor  series  expansion  which  are  not  included  in  the 
equation  for  the  predicted  surface,  but  are  included  in  the 
equation  for  the  actual  surface. 

QRIbQGQNALITY.  One  of  the  characteristics  of  the 
experimental  design  which  is  sought  by  the  experimenter  is 
the  mutual  orthogonality  of  coefficients.  Orthogonality  is 
that  property  of  an  experiment  which  ensures  that  the 
different  classes  of  effects  shall  be  capable  of  direct  and 
seperate  estimation  without  any  entanglement.  The  sum  of 
squares  of  all  the  effects  are  then  independent  and  additive 
(8:587) . 

Orthogonality  permits  surer  assessment  of  those  areas  on 
which  process  design  and  control  efforts  should  be 
concentrated.  It  also  provides  means  for  unambiguous 
simplification  and  improvement  of  response  models  along  with 


the  potential  for  uncovering  basic  mechanisms  (20:419). 

That  is,  those  terms  which  are  most  important  in  predicting 
the  surface  can  be  easily  identified. 

RQTATABILITY.  A  second  desirable  characteristic  of  an 
experimental  design  is  that  the  design  be  rotatable.  An 
experimental  design  is  said  to  be  rotatable  if  the  variance 
of  the  estimated  response  some  point  (x^, 

depends  on  the  distance  from  the  point  to  the  design  center 
and  not  on  direction.  In  other  words,  as  far  as  the  design 
is  concerned,  points  in  the  factor  space  which  are  the  same 
distance  from  the  origin  are  treated  as  being  equally 
important  (18: 165). 

Box  and  Draper  (3:339)  indicate  that  a  rotatable  design 
will  minimize  bias  provided  that  certain  conditions  are  met. 
These  conditions  are:  the  region  of  interest  is  spherical 
and  the  fitted  model  is  an  mth  order  polynomial  and  the  true 
model  is  an  nth  order  polynomial,  such  that  if  m+n=2r ,  then 
the  appropriate  design  is  an  rth  order  rotatable  design;  if 
m+n=2r+l,  then  the  appropriate  design  is  a  rth  order 
rotatable  design  with  moments  of  order  2r+l  all  zero.  Thus, 
in  fitting  a  second  order  polynomial  to  a  true  surface  of 
order  three,  the  fifth  moments  of  the  design  must  all  be 
zero  for  the  design  to  be  rotatable. 

NUMIiB  QE  design  PQINIS.  a  third  desirable 
characteristic  of  an  experimental  design  is  that  the  fewest 
possible  number  of  design  points  be  used  to  estimate  the 
response  surface.  The  fewest  number  of  design  points  is  a 


function  of  the  number  of  factors  and  the  type  of  design 
being  used. 

The  minimum  possible  number  of  design  points  that  can 
be  used  to  estimate  any  surface  is  defined  by  the  expression 
(k  +  1 )  ■•Kk+2) /2,  where  k  defines  the  number  of  factors  in  the 
experiment.  This  quantity  is  in  stark  contrast  to  the 
number  of  data  runs  which  are  required  if  the  experimenter 
chose  to  utilize  a  one  variable  at  a  time  procedure  to 
estimate  the  surface.  For  example,  in  a  simple  two  variable 
problem  in  which  the  first  variable  is  allowed  to  vary 
between  0  and  450  and  the  second  variable  is  allowed  to 
vary  between  0  and  750,  Manacapilli  calculated  that  338,701 
(451  *  751)  runs  are  required  to  produce  the  exact  surface 
of  a  two  factor  experiment  if  the  variables  are  varied  one 
at  a  time  (17:3.8).  On  the  otherhand,  certain  experimental 
designs  are  able  to  estimate  this  surface  with  only  six  data 
runs.  The  savings  in  time  and  cost  are  readily  apparent. 

ilAS  MINIMIZING  DESIGNS 

Since  the  only  error  associated  with  deterministic 
models  is  bias  error,  bias  minimizing  designs  are  very 
effective  in  fitting  response  surfaces  to  these  models.  The 
bias  minimizing  designs  which  are  utilized  in  this  project 
are:  Box -Draper  designs  and  Koshal  designs. 

iQX-DRAPER  BIAS  f11NI?1iZING  DESIGNS.  The  Box -Draper 
bias  minimizing  design  is  constructed  in  a  similar  manner  as 
a  central  composite  design.  That  is,  the  design  is  composed 
of  a  fractional  factorial  portion,  an  axial  point  portion 


and  center  points.  Like  the  central  composite  design,  the 
bias  minimizing  design  is  nearly  orthogonal  and  nearly 
rotatable. 

For  the  case  of  three  independent  variables  the  design 
matrix  is  given  by 
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The  specific  values  of  p  and  q  are  dependent  on  the 
second  and  fourth  moments  of  the  design  matrix.  The 
specific  values  which  minimize  bias  and  a  detailed 
derivation  of  the  design  matrix  can  be  found  in  Box  and 
Draper's  article  <2:622)  and  Myers'  book  (18:196). 

!£Q§b6L  designs.  The  Koshal  design  is  a  simple  design  used 
by  Koshal  in  fitting  frequency  distributions  by  the  method 
of  maximum  likelihood  (12:577).  The  Koshal  design  requires 
the  minimum  number  of  design  points  to  fit  a  quadratic 
response  surface,  and  the  major  advantage  of  using  this 
design  is  its  simplicity  and  economy. 

For  the  case  of  three  independent  variables  the  design 


matrix  is  given  by 


designs  which  are  utilized  in  this  project  are:  Central 
Composite  Rotatable  Designs,  Box-Behnken  Designs,  Hybrid 
Designs,  and  Minimum  Point  Designs. 


CENTRAL  COMPOSITE  ROTATABLE  DESIGNS.  The  central 


composite  rotatable  design  was  first  developed  by  Box  and 
Wilson  as  a  practical  alternative  to  the  3 factorial 
design  for  estimating  a  second  order  response  equation. 
This  design  is  probably  the  most  widely  known  and 
recommended  design  for  estimating  quadratic  response 
surfaces  (16:412). 

The  central  composite  design  is  the  2  factorial  or 
fractional  factorial,  augmented  by  axial  points  and  center 
points.  For  the  case  of  three  independent  variables  the 
design  matrix  is  given  by 
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(8) 


With  the  proper  selection  of  q,  the  central  composite  design 
is  orthogonal  and  rotatable.  A  detailed  derivation  of  the 
q-values  can  be  found  in  Box  and  Hunter’s  article  (4:195). 

IQX-BEHNKEN  designs.  In  1960,  Box  and  Behnken 
developed  a  new  class  of  three  level  factorial  designs  which 
were  useful  for  estimating  the  coefficients  in  a  second 
degree  graduating  polynomial  (1:455).  These  designs  are 


a  small 


nearly  rotatable  and  nearly  orthogonal  (i.e., 
correlation  exists  between  the  constant  term  and  the  squared 
terms  and  between  the  squared  terms  themselves) .  A 
practical  advantage  of  the  Box-Behnken  design  over  the  full 
central  composite  design  is  that  when  evaluating  problems 
with  large  number  of  factors  (i.e.,  greater  than  five),  the 
Box-Behnken  design  requires  fewer  design  points  to  estimate 
the  surface. 

The  designs  are  formed  by  combining  factorial  designs 
with  incomplete  block  designs  (1:457).  For  the  case  of 
three  independent  variables,  the  design  matrix  is  given  by 
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A  detailed  derivation  of  the  design  matrix  can  be  found  in 
Box  and  Behnken’s  article  (1:455). 

Hybrid  designs  were  created  to  achieve 
the  same  degree  of  orthogonality  as  central  composite 
designs  or  regular  polyhedral  designs,  to  be  nearly  minimum 
point,  to  be  nearly  rotatable,  and  to  possess  some  ease  in 
coding  (20:419). 

The  hybrid  designs  are  constructed  in  a  similar  manner 
as  central  composite  designs.  That  is,  the  design  is 
composed  of  a  fractional  factorial  portion,  an  axial  point 
portion,  and  center  points.  The  major  difference  between 
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the  hybrid  and  central  composite  designs  is  that  hybrid 
designs  are  augmented  with  an  extra  variable  column  which 
resembles  a  cross  polytope  design.  For  the  case  of  three 
independent  variables  the  design  matrix  is  given  by 
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A  detailed  derivation  of  the  design  matrix  can  be  found  in 
Roquemore's  article  (20:419). 

mNIMUri  PgiNI  second  order  designs.  The  purpose  of 

second  order  designs  with  minimum  number  of  points  is  to 
provide  a  low  cost,  practical  design  to  estimate  a  quadratic 
response  surface.  For  a  k— factor  design,  the  minimum  number 
of  points  is  <k+l ) * (k+2) /2.  The  design  should  give  rise  to 
least  squares  estimates  with  minimum  generalized  variance 
(6:613) . 

For  the  case  of  three  independent  variables,  the  design 


matrix  is  given  by 
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(11) 


The  specific  values  of  p  and  q,  along  with  a  detailed 
derivation  of  the  design  matrix  can  be  found  in  Box  and 
Draper's  article  (6:613). 


MISCELLANEOUS 

Several  articles  have  been  published  which  have  offered 
an  alternate  approach  to  minimizing  bias  error  in  response 
surfaces.  The  articles  by:  1)  Karson,  Manson ,  and 
Hader ( 14: 461 ) ,  2)  Karson  (13:1565),  3)  Karson  and  Spruill 
(15:329),  and  4)  Cote,  Manson  and  Hader  (7:633),  have 
applied  similar  techniques  to  address  the  problem  of 
minimizing  bias.  That  is,  their  methods  consist  of 
developing  an  estimator  for  a  given  design  which  minimizes 
bias. 

Although  these  works  will  not  be  used  in  this  project, 
their  availability  is  noted  here  for  any  future  work  on  this 
topic. 
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NUCLEAR  EXCHANGE  MODEL 


The  deterministic  model  used  to  generate  the  response 
values  for  this  study  is  a  nuclear  exchange  linear 
programming  model.  In  this  rather  simplistic  model  of  a 
nuclear  force  structure,  the  single  objective  of  the  model 
is  to  maximize  the  damage  expectancy  of  the  force  structure 
The  damage  expectancy  is  maximized  given  the  effectiveness 
of  each  weapon  type  against  each  target  type. 

For  example,  consider  the  hypothetical  strategic  force 
of  three  weapon  systems  (ICBM,  Bomber,  SLBM)  and  two  target 
types.  Table  1  shows  the  amount  of  damage  expected  (Dl,  D2 
D3,...,  D6)  given  a  particular  weapon  type  and  target  type. 

TABLE  1 

SINGLE  SHOT  DAMAGE  EXPECTANCY  TABLE 

WiBEQN  lYPE  IFACIQRSl 

ICBM  (Ml)  BOMBER  (Bl)  SLBM  (SI) 


TGT  TYPE 

TARGET  1 

Dl 

D2 

D3 

TARGET  2 

D4 

D5 

D6 

The  damage  expectancy  values  in  Table  1  represent  single 
shot  damage  expectancv,  that  is,  the  damage  expectancy  from 
a  single  warhead  of  a  particular  weapon  against  a  single 
target.  In  this  model,  an  additional  warhead  could  be 
launched  at  the  same  target.  Therefore,  it  is  possible  to 
increase  the  damage  expectancy  against  a  single  target  by 


launching  two  warheads  <of  the  same  or  di-F-ferent  weapon 
type)  at  the  target.  In  those  cases  where  two  warheads  are 
launched  at  the  same  target,  the  new  damage  expectancy  value 
■for  that  target  can  be  computed  in  the  following  manner 

°^NEW=  1  “  C<1-DE,  )*(1-DE2  )]  <12) 

Thus,  there  are  six  additional  combinations  of  damage 
expectancy  which  must  be  considered  for  each  target  type 
(e.g.,  M1*S1,  Bl-M-Sl,  S1*S1).  After 

considering  all  possible  combinations,  the  LP  model  now 
maximizes  the  total  damage  expectancy  subject  to  the 
limitations  on  the  number  of  weapons,  number  of  targets,  and 
other  targeting  constraints.  Further  details  of  the 
specific  model  used  to  generate  the  data  along  with  the  LP 
formulation  are  available  in  Appendix  A. 

APPLICAIIQN  of  IHE  EXPERIMENIAL  design 

As  mentioned  earlier,  there  is  a  practical  advantage  in 
employing  the  theory  of  experimental  design  in  estimating  a 
response  surface.  For  example,  in  the  simple  three  factor 
problem  used  in  this  study,  if  each  factor  were  allowed  to 
vary  one  at  a  time,  the  total  number  of  runs  required  to 
generate  the  surface  is  14,601,339,651  (  1351'N-18O1*6001  )  .  By 
constructing  the  appropriate  experimental  design,  the 
surface  could  be  estimated  with  as  few  as  ten  runs. 

The  experimental  designs  which  are  used  in  this  study 
are  second  order  designs.  The  purpose  of  a  second  order 
design  is  to  fit  a  second  order  response  equation  to  the 


data.  The  second  order  response  equation  -for  the  three 
■factor  case  is  of  the  form 
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(13) 

where,  Ml  =  number  of  ICBM  warheads,  B1  =  number  of  bomber 
warheads,  and  SI  =  number  of  SLBM  warheads. 

One  of  the  requirements  to  produce  a  second  order  design 
is  that  at  least  three  levels  of  each  factor  must  be 
selected  so  that  the  coefficients  of  the  response  equation 
can  be  estimated.  The  three  levels  which  are  usually 
selected  are  the  values  at  the  upper  end  (having  a  coded 
value  of  •♦•1)  and  lower  end  (having  a  coded  value  of  -1)  of 
the  range  for  each  input  parameter  plus  the  mid-point 
(having  a  coded  value  of  0)  between  the  extreme  points. 
These  points  are  selected  and  coded  in  this  manner  so  that 
the  coefficients  of  the  response  equation  are  not  highly 
correlated. 

Table  2  indicates  the  uncoded  and  coded  values  of  the 
three  factor  problem  used  in  this  project. 


TABLE  2 


FACTOR  LEVELS 


WEAPON  TYPES  (FACTORS) 

LEVELS 

(CODED  VALUES) 

Ml 

1650(-1> 

2325(0) 

3000  (-Hi) 

Bi 

1200(-1> 

2100(0) 

3000(^1) 

SI 

6000 (-1) 

9000(0) 

12000  (-Hi) 

In  many  o-f  the  central  composite  designs  which  are  used 
in  this  study,  the  coded  values  are  greater  than  +1  and  less 
than  -1.  The  appearance  of  these  values  can  best  be 
explained  by  examining  the  central  composite  design  for  the 
three  factor  case.  In  this  situation,  the  central  composite 
design  selects  its  initial  design  points  (factorial  portion) 
at  the  corners  of  a  cube  centered  at  the  origin.  Because 
the  region  of  interest  is  taken  to  be  spherical,  the  design 
positions  additional  points  along  the  axes  so  that 
representative  points  can  be  sampled  over  the  entire  region. 
These  axial  design  points  lie  near  the  boundary  of  the 
sphere  and  along  one  of  the  axes.  Therefore,  the  axial 
points  have  values  greater  than  +1  and  less  than  -1  due  to 
the  radius  of  the  sphere. 

There  is  a  possibility  of  significant  error  in  the 
estimated  response  values  because  the  axial  points  are 
positioned  outside  the  original  feasible  region.  That  is, 
if  the  true  surface  is  radically  different  outside  the 
feasible  region  than  it  is  within  the  feasible  region,  a 
significant  error  occurs  due  to  the  positioning  of  the 
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axial  points.  In  that  case,  a  di-f-Ferent  means  of  coding  or 
a  different  design  <which  does  not  utilize  points  outside 
this  region)  should  be  employed.  However,  if  the  true 
surface  is  not  radically  different  outside  the  feasible 
region  than  it  is  within  the  feasible  region,  then  the  axial 
design  points  do  not  contribute  a  significant  error  in 
estimating  the  predicted  surface. 

Furthermore,  many  of  the  designs  which  have  been 
selected  in  this  study  have  been  originally  constructed  with 
multiple  center  points.  Often  times,  multiple  center  points 
are  used  to  minimize  variance,  to  achieve  rotatability  or  to 
achieve  orthogonality.  However,  the  replication  of  the 
center  points  is  not  required  while  analyzing  deterministic 
models  because  there  is  no  sampling  error  associated  with 
this  type  of  model.  Although  pure  rotatability  and 
orthogonality  may  be  sacrificed  in  these  cases,  the  designs 
still  maintain  a  high  degree  of  rotatability  and 
orthogonality.  Therefore,  only  one  center  point  is  included 
in  any  design  which  has  been  originally  constructed  with 
multiple  center  points. 

The  designs  which  are  used  in  this  study  for  the  three 
factor  case  have  already  been  outlined.  The  designs  which 
are  used  in  the  three,  four,  five,  and  six  factor  problems 
are  provided  in  Appendix  B. 

After  selecting  the  appropriate  experimental  design  and 
after  defining  the  LP,  the  next  step  in  the  methodology  is 
to  execute  the  appropriate  number  of  force  allocation  runs 


using  the  LP  model.  One  run  is  accomplished  for  each  of  the 
identified  combinations  in  order  to  obtain  the  optimum 
damage  expectancy  for  each  set  of  input  factor  values. 

The  final  steps  in  obtaining  a  response  equation  is  to 
use  the  results  of  a  particular  design  as  input  to  a 
multiple  linear  regression  program.  The  coded  levels  of 
each  of  the  factors  are  input  to  the  program  as  the 
independent  variables,  and  the  damage  expectancy  is  input  as 
the  dependent  variable.  Using  multiple  stepwise  linear 
regression,  the  coefficients  of  the  response  equation  are 
determined.  The  coded  values  are  then  uncoded  in  order  to 
produce  a  relationship  in  which  the  actual  number  of 
warheads  of  the  weapon  systems  can  be  used.  The  response 
equation  for  uncoded  values  portrays  the  relationship  of 
each  factor  to  the  response  variable.  A  sample  of  the  input 
and  the  output  for  the  stepwise  linear  regression  is  found 
in  Appendix  C. 

The  next  chapter  portrays  the  results  and  details  the 
analysis  of  the  effectiveness  and  efficiency  of  the  various 
designs  used  in  this  study. 
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Ill.  RESULTS  AND  ANALYSIS 


INIRODUCTION 

This  chapter  presents  the  analysis  of  the  response 
equations  for  the  various  designs  which  are  used  in  this 
study.  The  purpose  of  this  analysis  is  not  to  evaluate  the 
LP  but  rather,  to  evaluate  the  designs  and  the  response 
equations. 

The  analysis  is  divided  into  two  parts:  the  predictive 
power  of  the  response  equations  and  the  explanatory  power  of 
the  response  equations.  The  predictive  power  of  the 
response  equation  addresses  the  equation's  ability  to 
adequately  estimate  the  true  surface  and  to  accurately 
predict  a  response  for  a  given  set  of  inputs.  The 
explanatory  power  addresses  the  equation's  ability  to 
present  a  response  equation  which  is  simple  to  interpret  so 
that  the  true  surface  is  easily  evaluated  and  the  results 
are  easily  explained. 

PRiDICIlVE  POWER  OF  RESPONSE  SURFACE  EQUATIONS 

As  stated  in  Chapter  I,  the  basis  for  evaluating  the 
predictive  capabilites  of  an  equation  of  a  response  surface 
design  is  Bq(4),  the  effectiveness  -  efficiency  factor. 

E=(error  in  fit  of  design) * (number  of  design  points)  (4) 
This  factor  generates  a  value  so  that  various  designs  are 


compared  in  terms  of  their  precision  in  estimating  the 


true  surface  and  the  number  of  required  experimental  design 
poi nts. 

In  order  to  obtain  the  error  in  fit  of  the  design,  a 
random  set  of  points  is  used  as  inputs  into  the  response 
surface  equation  and  the  differences  between  the  actual  and 
predicted  responses  are  recorded.  In  comparing  these 
response  equations,  the  random  data  set  is  not  used  to 
estimate  any  response  equation,  but  rather  to  compare  the 
predicted  values  of  the  estimated  equations  against  each 
other.  The  benefit  of  the  random  data  set  is  that  all 
equations  are  compared  against  a  common  data  set  so  that 
prediction  bias  is  avoided.  Prediction  bias  arises  when  the 
choice  of  a  regression  equation  is  uniquely  related  to  the 
observations  from  which  it  was  created. 

The  response  equations  which  are  used  to  evaluate  the 
error  in  fit  of  the  design  are  modified  versions  of  the  full 
regression  equations.  That  is,  selected  terms  from  the  full 
regression  equations  have  been  omitted  from  the  modified 
regression  equaitons  because  these  terms  account  for  less 
than  two  percent  of  the  total  sum  of  squares. 

The  primary  reason  for  excluding  these  terms  is  to 
obtain  a  relatively  simple  equation  (i.e. ,  an  equation  with 
as  few  terms  as  possible)  which  is  easy  to  interpret,  but 
still  produces  a  good  fit.  If  the  full  regression  equation 
is  used  instead  of  the  modified  equation  to  estimate  the 
true  surface,  the  ability  of  the  response  equation  to 
predict  future  response  values  is  improved;  however,  the 
ability  to  interpret  the  equation  is  significantly 


degraded.  That  is,  the  interpretation  of  the  equation  is 
very  difficult  due  to  the  complexity  of  the  equation.  Thus, 
by  selecting  two  percent  of  the  total  sum  of  squares  as  the 
criteria  for  determining  what  terms  to  include  in  the 
response  equation,  a  reasonably  good  fit  is  obtained  while 
still  preserving  the  interpretative  value  of  the  response 
equation. 

The  response  surface  equations  along  with  the  random 
data  sets  which  are  used  in  the  analysis  of  the  predictive 
power  can  be  found  in  Appendix  D. 

Tables  3  through  6  depict  the  results  of  the  predictive 
power  of  the  various  designs  which  are  used  in  this  study. 

TABLE  3. 

PREDICTIVE  POWER  OF  THREE  FACTOR  DESIGNS 

PERCENT  ERROR  IN  FIT  NUMBER  OF 
design  ig  RANDOM  DATA  SET  DESIGN  POINTS  E  r  yBLUE 

VARIANCE  DESIGNS 


BOX-BEHNKEN 

* 

1.08 

13 

14.04 

HYBRID  310+CP 

* 

1.28 

11 

14.08 

HYBRID  311 A 

* 

1.63 

1 1 

17.93 

CENTRAL  COMP 

* 

1.56 

15 

23.40 

MIN  POINT 

6.  12 

10 

61.20 

BIAS  DESIGNS 

BOX-DRAPER 

* 

1.38 

15 

20.70 

KOSHAL 

2.74 

10 

27.40 

*  indicates  nearly  rotatable  design 


TABLE  4 


PREDICTIVE  POWER  OF  FOUR  FACTOR  DESIGNS 

PERCENT  ERROR  IN  FIT  NUMBER  OF 

design  IQ  RANDOM  DAIA  SEX  DESIGN  PQINIS  E  z 

VARIANCE  DESIGNS 


HYBRID  416A+CP 

* 

1.06 

17 

HYBRID  416C 

* 

1.80 

16 

BOX-BEHNKEN 

* 

1.18 

25 

CENTRAL  COMP 

* 

1.31 

25 

MIN  POINT 

4.37 

15 

BIAS  DESIGNS 

BOX-DRAPER 

* 

0.96 

25 

KOSHAL 

1.80 

15 

*  indicates  nearly 

rotatable  design 

VALUE 

18.02 
28.80 
29.50 
32.75 
65.  55 

24.00 

27.00 


TABLE  5 


PREDICTIVE  POWER  OF  FIVE  FACTOR  DESIGNS 


PERCENT  ERROR  IN  FIT 

NUMBER  OF 

design  IQ  RANDOM 

daia  sex 

DESIGN  POINTS 

i  =  VALUE 

VARIANCE  DESIGNS 

CENTRAL  COMP  * 

1.73 

27 

46.71 

(HALF  REPLICATE) 
BOX-BEHNKEN  * 

1.41 

41 

57.81 

MIN  POINT 

5.88 

21 

123.48 

BIAS  DESIGNS 

KOSHAL 

1.15 

21 

24.  15 

BOX -DRAPER  * 

1.37 

27 

36.99 

(HALF-REPLICATE) 
BOX -DRAPER  * 

1.31 

43 

56.33 

(FULL) 

*  indicates  nearly  rotatable  design 


TABLE  6. 


PREDICTIVE  POWER  OF  SIX  FACTOR  DESIGNS 

PERCENT  ERROR  IN  FIT  NUMBER  OF 
design  IQ  RANeOtJ  DAIA  SET  DESIGN  PQINIS  E  z  VALUi 

VARIANCE  DESIGNS 


HYBRID  628 A 

* 

1.33 

28 

37.24 

HYBRID  628B 

1.35 

28 

37.80 

MIN  POINT 

1.76 

28 

49.28 

BOX-BEHNKEN 

* 

1.15 

49 

56.35 

CENTRAL  COMP  * 
(HALF-REPLICATE) 

1.33 

45 

59.85 

BIAS  DESIGNS 

KOSHAL 

1.72 

28 

48.  16 

BOX-DRAPER 

* 

1.39 

45 

62.55 

(HALF-REPLICATE) 

*  indicates  nearly  rotatable  design 

ACCURACY.  Tables  3  through  7  indicate  that  the 
accuracy  o-f  a  response  sur-face  equation  is  dependent  on  the 
number  o-f  design  points  and  on  the  degree  of  rotatability  of 
the  design.  That  is,  the  accuracy  of  a  response  surface 
equation  is  highly  correlated  to  the  number  of  design 
points,  but  this  relationship  is  not  a  perfect  correlation 
and  their  differences  in  fit  may  be  explained  by  the  degree 
of  the  rotatability  of  the  designs. 

For  example,  it  may  be  hypothesized  that  the  design 


with  the  most  points  produces  the  most  accurate  fit.  That 
is,  each  design  selects  its  design  points  under  the 
assumption  that  the  design  points  are  representative  of 


their  surrounding  area.  In  theory,  by  selecting  more  design 
points,  the  area  for  which  each  point  is  an  estimate  is 
decreased..  Thus,  the  total  error  in  the  estimate  of  the 
true  surface  is  reduced.  But  as  Tables  3  through  7 
indicate,  more  design  points  does  not  guarantee  a  more 
accurate  fit.  This  result  is  rather  surprising,  but  it  does 
indicate  that  other  contributing  factors  influence  the 
degree  of  accuracy  of  the  response  equation  other  than  the 
number  of  design  points. 

A  major  contributing  factor  which  influences  the 
accuracy  of  an  equation  of  a  response  surface  design  may  be 
the  degree  of  rotatability  of  the  design.  In  general,  the 
variance  and  bias  minimizing  designs  are  able  to  generate 
similar  degrees  of  fit  in  estimating  a  response  surface 


equation  to  the  true  surface,  provided  that  the  variance 
minimizing  designs  are  nearly  rotatable.  Rotatability  is  an 
important  characteristic  for  variance  minimizing  designs 
when  estimating  deterministic  models  because  the  selection 
of  the  design  points  is  not  strictly  based  on  minimizing 
bias.  But,  as  stated  in  Chapter  I,  if  the  region  of 
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interest  is  spherical,  rotatable  designs  minimize  bias.  The 
importance  of  rotatability  in  minimizing  bias  is  highlighted 
by  the  fact  that  the  only  variance  minimizing  design  which 
fails  to  consistently  produce  a  relatively  accurate  fit  is 
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the  minimum  point  design,  which  is  a  non-rotatable  design. 

However,  none  of  the  designs  mentioned  thus  far  is 
perfectly  rotatable;  thus,  bias  is  not  totally  minimised  for 
these  designs.  In  fact,  all  of  the  designs,  except  the 
minimum  point  and  the  Koshal  designs,  are  nearly  rotatable. 
That  is,  these  nearly  rotatable  designs  have  varying  degrees 
of  rotatability  and  thus,  have  varying  degrees  of 
effectiveness  in  minimzing  bias.  Therefore,  it  might  be 
hypothesized  that  the  differences  in  accuracy  between  the 
designs  which  are  not  accounted  for  by  the  number  of  design 
points,  are  accounted  for  by  the  degrees  of  rotatability  of 
the  designs. 

For  example,  the  response  surface  equation  of  the 
hybrid  design  <416A+CP)  achieves  a  more  accurate  fit  to  the 
true  surface  than  the  response  equation  of  the  Box-Behnken 
design  for  the  four  factor  problem.  In  this  example,  the 
hybrid  design  requires  eight  fewer  points  than  the  Box- 
Behnken  design,  but  is  able  to  achieve  a  higher  degree  of 
accuracy.  By  the  hypothesis,  this  indicates  that  the  hybrid 
design  (416A+CP)  has  a  higher  degree  of  rotatability  than 
the  Box-Behnken  design  for  four  factors.  This  hypothesis  is 
confirmed  by  the  fact  that  the  Box-Behnken  design  is  not 
highly  rotatable  for  four  factors,  while  the  hybrid  design 
maintains  a  high  degree  of  rotatability. 

A  similar  argument  can  be  made  to  explain  why  no 
single  design  is  able  to  achieve  the  most  accurate  fit  for 
all  factors.  That  is,  for  the  different  number  of  factors. 


these  designs  have  di-F-ferent  number  of  design  points  and 
different  degrees  of  rotatabi 1 i ty.  Thus,  the  degrees  of 
accuracy  for  the  response  equations  are  affected. 

iEEieiiVENESS-EFFICIENCY  FACIQR.  Since  the  accuracy 
of  fit  of  the  design  is  highly  correlated  to  the  number  of 
design  points,  the  results  which  are  obtained  concerning  the 
ef f ecti venss-ef f iciency  factor  are  not  surprising.  That  is, 
a  design  which  provides  the  most  accurate  fit  with  the 
fewest  number  of  design  points  does  not  exist.  Even  a 
perfectly  rotatable  design  does  not  provide  the  most 
accurate  fit,  if  a  sufficient  number  of  design  points  are 
not  specified. 

Thus,  in  attempting  to  minimize  the  effectiveness- 
efficiency  factor,  it  remains  for  the  experimenter  to 
evaluate  the  problem  to  determine  what  degree  of  accuracy  is 
required  and  how  many  design  points  are  available.  In 
evaluating  the  trade-off  between  accuracy  and  the  number  of 
design  points,  one  major  consideration  is  the  cost  of 
performing  additional  data  runs.  In  those  problems  in  which 
cost  of  executing  additional  data  runs  is  minimal,  the 
experimenter  may  concentrate  on  achieving  a  higher  degree  of 
accuracy.  On  the  other  hand,  in  those  cases  in  which  set-up 
costs  and  run  time  are  prohibitive,  the  experimenter  may 
select  the  most  economical  design  (in  terms  of  the  fewest 
design  points)  which  achieves  an  adequate  level  of  accuracy. 

The  trade-offs  between  accuracy  (level  of  fit)  and  cost 
(number  of  design  points)  are  further  complicated  when  the 


experimenter  desires  a  response  equation  which  has  superior 
explanatory  powers. 

iXPLANATORY  POWER  OF  RESPONSE  SURFACE  EQUAJIDN 

The  basis  for  evaluating  the  explanatory  capabilities 
of  an  equation  of  a  response  surface  design  is  the 
interpretation  of  the  coefficients  of  the  response  equation. 
The  coefficients  provide  answers  to  such  questions  as: 

“What  is  the  overall  effect  to  the  dependent  variable, 
if  one  independent  variable  is  increased  by  an  amount  X 
and  all  other  independent  variables  are  held  constant?" 
By  applying  the  analysis  techniques  of  multiple  regression, 
similar  “what. . . if . . . “  questions  can  be  answered  quickly  and 
accurately.  The  answers  to  these  questions  are  quick  and 
easy  because  the  basic  interpretation  of  the  coefficients  in 
multiple  regression  analysis  is  that  the  coefficients  of  the 
main  effect  terms  represent  the  marginal  effects  of  change 
of  the  dependent  variable  when  the  values  of  the  independent 
variables  are  altered. 

iEEiCI  OF  NULIICOLLINEARIIY  ON  REGRESSION 
COEFFICIENTS.  As  stated  earlier,  Box  and  Hunter  have 
identified  a  small  correlation  that  exists  between  the 
constant  term  and  the  squared  term  and  between  the  squared 
terms  themselves  (4:201).  This  correlation  is  commonly 
referred  to  as  multicol 1 inearity  in  statistical  literature. 
Since  Box  and  Hunter  concentrated  on  the  predictive  power 
of  the  response  equation,  multicol 1 inearity  did  not  affect 


their  results  because  of  the  fact  that  when  some  or  all 
independent  variables  are  correlated  among  themselves,  the 
ability  to  obtain  a  good  fit  is  not  generally  inhibited  nor 
does  it  tend  to  affect  prediction  of  new  observations 
<19:384).  That  is,  the  predictive  power  of  the  response 
equation  is  not  generally  affected  when  multicoll inearity 
exists. 

However,  the  common  interpretation  of  regression 
coefficients  as  measuring  the  change  in  the  expected  value 
of  the  dependent  variable  when  the  corresponding  independent 
variable  is  increased  by  one  unit  while  all  independent 
variables  are  held  constant  is  not  fully  applicable  when  a 
high  degree  of  multi col 1 inearity  exists.  While  it  may  be 
conceptually  possible  to  vary  one  independent  variable  and 
hold  the  others  constant,  it  may  not  be  possible  in  practice 
to  do  so  for  independent  variables  that  are  highly 
correlated  (19s385).  Because  a  high  degree  of  correlation 
among  the  independent  variables  exists  in  non-orthogonal 
designs,  these  designs  do  not  possess  a  high  degree  of 
explanatory  power  by  the  fact  that  the  regression 
coefficients  of  the  response  equations  do  not  accurately 
measure  the  marginal  effects  of  change  of  the  dependent 
variable  when  the  values  of  the  independent  variables  are 
altered.  Thus,  only  those  designs  which  have  a  high  degree 
of  orthogonality  will  be  considered  in  the  explanatory  power 
section  of  this  study.  The  designs  which  are  to  be 
evaluated  are:  the  Central  Composite  Rotatable  designs,  the 


Box-Behnken  designs,  the  Hybrid  designs,  and  the  Box-Draper 
Bias  Minimizing  designs. 


IHREE  factor  RESULIS.  Equation  (14)  is  the  basic 
response  equation  -for  the  three  factor  problem. 

DE  =  bg  +  b,(Ml)  +  bjCBl)  +  bjCSl) 

-  b33  (S1*S1)  -  bjj  (B1*S1>  (14) 

The  terms  which  have  been  omitted  from  the  equation 
account  for  less  than  two  percent  of  the  total  sum  of 
squares,  and  therefore  are  dropped  from  the  final  equation. 

Table  7  lists  the  coefficients  of  the  response  equation 
for  the  various  designs  for  the  three  factor  problem. 

TABLE  7. 

COEFFICIENTS  OF  THE  RESPONSE  EQUATION 
FOR  THREE  FACTOR  DESIGNS 

DESIGN  COEFFICIENTS 


bo 

bi 

bz 

b3 

b33 

b  23 

(1)  Central 
Comp 

124.9285 

.4848 

1.  153 

1.004 

. 0000276 

. 000062 

(2)  Box- 

Behnken 

-648.3319 

.4965 

1.225 

1.  180 

. 0000375 

. 0000665 

(3)  Hybrid 
(311A) 

68.4118 

.5243 

1.  195 

.961 

.0000251 

. 0000602 

(4)  Hybrid 
(310+CP) 

-626. 1627 

.5002 

1.327 

1.  163 

. 0000352 

. 0000774 

(5)  Box- 
Draper 

-285.6179 

.4798 

1.223 

1.085 

.0000318 

. 0000667 

EQUB  EeQIQq  BESULTS.  Equation  (15)  is  the  basic 
response  equation  for  the  four  factor  problem. 

DE  =  b„  +  b,(Ml)  +  b  2(M2)  +  b3(Bl)  +  b4<Sl) 

-  b4  4  (S1*S1)  -  b24  (M2*S1)  -  b34  (B1*S1)  (15) 

The  terms  which  have  been  omitted  from  the  equation 
account  for  less  than  two  percent  of  the  total  sum  of 
squares,  and  therefore  are  dropped  from  the  final  equation. 

Table  8  lists  the  coefficients  of  the  response  equation 
for  the  various  designs  for  the  four  factor  problem. 


DESIGN 


TABLE  8. 

COEFFICIENTS  OF  THE  RESPONSE  EQUATION 
FOR  FOUR  FACTOR  DESIGNS 

COEFFICIENTS 


(1) 

Central 

Comp 

939.3 

(2) 

Box- 

Behnken 

1320.0 

(3) 

Hybrid 

(4160 

1129. 1 

(4) 

Hybrid  >1060.1 
(416A+CP) 

(5) 

Box- 

Draper 

1332.3 

“l 

*53 

C4 

O44 

>^24 

“34 

.416 

.961 

1.02 

.907 

. 000023 

. 000054 

.000051  ■ 

.383 

.946 

.987 

.877 

. 000023 

. 000054 

. 000051 

.434 

1.12 

1.08 

.793 

.000015 

. 000068 

. 0000645 

.439 

1.09 

1.11 

.823 

. 000017 

. 000064 

.0000601 

.395 

.944 

.988 

.887 

. 000024 

. 000053 

. 0000503 

FWE  FAQXQR  results.  Equation  (16)  is  the  basic 


response  equation  for  the  five  factor  problem. 

DE  =  bo  +  b,(Ml)  +  b2<M2)  +  boCBl) 

+  b4(B2)  +  bg(Sl)  -  bggCSl^Sl)  (16) 

The  terms  which  have  been  omitted  from  the  equation 
account  for  less  than  two  percent  of  the  total  sum  of 
squares,  and  therefore  are  dropped  from  the  final  equation. 
Table  9  lists  the  coefficients  of  the  response 
equation  for  the  various  designs  for  the  five  factor 
problem. 


TABLE  9. 

COEFFICIENTS  OF  THE  RESPONSE  EQUATION 
FOR  FIVE  FACTOR  DESIGNS 

DESIGN  COEFFICIENTS 


bo 

^2 

b4 

bs 

b  5  5 

(1)  Central 
Comp 

1615.7 

.501 

.  630 

.712 

.755 

.951 

.0000261 

(2)  Box- 

Behnken 

844.0 

.603 

.680 

.729 

.769 

1.04 

. 0000295 

(3)  Box- 

Draper  (H) 

536.6 

.604 

.  683 

.739 

.785 

1.09 

. 0000329 

(4)  Box- 

Draper  (F) 

388.0 

.629 

.678 

.756 

.798 

1.09 

. 0000329 
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SIX  FACTOR  RESULTS.  Equation  (17)  is  the  basic 
response  equation  -for  the  six  -factor  problem. 


DE  =  b^  +  b,(Ml)  +  bjCMZ)  +  b3<M3) 


+  b4(Bl)  +  b5<B2)  +  bgCSl) 


-  b5g(Sl*Sl)  -  bgg  (B2*S1) 


The  terms  Mhich  have  been  omitted  from  the  equation 
account  for  less  than  two  percent  of  the  total  sum  of 
squares,  and  therefore  are  dropped  from  the  final 
equation. 

Table  10  lists  the  coefficients  of  the  response 
equation  for  the  various  designs  for  the  six  factor 
problem. 

TABLE  10. 

COEFFICIENTS  OF  THE  RESPONSE  EQUATION 
FOR  SIX  FACTOR  DESIGNS 


DESIGN 


COEFFICIENTS 


(1)  Central  | 

1  620. 1 

Compos i te 

(2)  Box- 

-642. 1 

Behnken 

(3)  Hybrid 

715.9 

(628A) 

(4)  Hybrid 

695.9 

(628B) 

(5)  Box- 

-101.0 

Draper 

°2 

°3 

D4 

°66 

°56 

.562 

.606 

.673 

.685 

1.24 

1.05 

. 0000235 

. 000057 

.536 

.628 

.654 

.692 

1.39 

1.28 

. 0000331 

. 000073 

.563 

.598 

.626 

.684 

1.23 

1.04 

. 0000228 

. 000057 

.554 

.606 

.614 

.  683 

1.23 

1.04 

. 0000232 

. 000057 

.568 

.645 

.678 

.695 

1.22 

1.21 

. 0000328 

. 000054 

INIiRPBilAIiQN  OF  CQEFFICIENIS.  The  parameter  is 
the  y-intercept  of  the  regression  surface.  This  term  does 
not  have  any  particular  meaning  as  a  separate  term  in  any  of 
the  regression  models  because  the  range  of  the  variables  in 
the  model  does  not  include  the  point  where  lil=Q, .  .  .  ySl=(9. 

If  this  point  had  been  in  the  relevant  range  of  the 
variables,  the  b^  term  would  indicate  the  response  level 
when  all  of  the  independent  variables  are  set  to  zero. 

When  the  regression  model  does  not  include  cross 
product  or  squared  terms  (i.e. ,  a  first  order  equation), 
the  parameters  b  (k^0)  indicate  the  change  in  the  mean 
response  per  unit  increase  in  when  all  other  variables 
included  in  the  model  are  held  constant.  That  is,  b^ 
indicates  the  change  in  the  expected  value  of  the 
response  with  a  unit  increase  in  x,  when  all  other 
variables  in  the  model  are  held  constant.  Since  the 
experimental  designs  which  are  being  considered  in  this 
section  are  nearly  orthogonal,  the  independent  variables 
are  almost  uncorrelated  and  the  variables  have  an 
additive  effect. 

When  the  regression  model  includes  cross  product  or 
squared  terms,  such  as  Eq  (14)  of  the  three  factor  model, 
the  coefficients  of  the  main  effect  terms  no  longer  indicate 
the  change  in  the  mean  response  per  unit  increase  in  x^  for 
any  given  level  of  the  remaining  variables,  if  that  variable 
is  also  included  in  the  cross  product  or  squared  terms.  For 
example,  in  Eq(14),  b,  no  longer  indicates  the  change  in  the 


mean  response  for  a  unit  increase  in  B1  because  of  the 
interaction  of  B1  and  SI.  Similarly,  bg  no  longer  indicates 
the  change  in  the  mean  response  for  a  unit  increase  in  SI 
because  of  the  interaction  of  SI  and  B1  and  because  of  the 
interaction  of  SI  with  itself.  Because  of  the  complexity  of 
a  regression  model  which  includes  several  cross  product  and 
squared  terms,  the  most  desirable  regression  equation  which 
explains  the  response  surface  is  one  in  which  the  main 
effects  dominate  and  the  cross  product  and  squared  terms 
terms  take  minor,  secondary  roles.  In  this  manner,  the 
parameters  b  ^  would  indicate  the  change  in  the  mean  response 
per  unit  increase  in  x,^  when  all  other  variables  included  in 
the  model  are  held  constant. 

Differences  between  Coefficients.  Tables  7 
through  10  indicate  that  no  two  equations  have  the  same 
coefficients  to  estimate  the  same  surface.  That  is,  each 
design  produces  a  unique  response  equation  to  estimate  the 
true  surface.  These  results  are  not  unexpected  and  does  not 
disqualify  the  use  of  multiple  regression  analysis 
techniques  as  an  analytical  tool  when  evaluating  response 
surfaces  of  deterministic  models. 

For  example,  small  differences  between  the  coefficients 
are  to  be  expected  because  the  predicted  surface  is  not  a 
perfect  representation  of  the  actual  surface.  The  errors  in 
the  predicted  surface  are  due  to  two  factors:  first,  each 
design  selects  different  quantities  of  sample  points  to 
estimate  the  true  surface  and  each  design  selects  its  design 


points  at  di-fferent  locations  in  the  region  of  interest;  and 
second,  a  small  amount  of  mul ti col 1  inear ity  is  still  present 
between  some  of  the  variables  in  these  equations. 

The  selection  of  the  number  of  design  points  and  the 
placement  of  the  design  points  have  a  tremendous  impact  on 
the  fitted  regression  function  because  the  response  equation 
is  influenced  by  whatever  design  points  are  selected. 

Because  each  of  the  designs  selects  points  from  different 
areas  of  the  true  surface,  the  response  equation  must 
necessarily  be  different  for  each  design.  Any  change  in  the 
selection  of  the  design  points  influences  the  coefficients 
of  the  response  equation  to  accomodate  this  change. 

For  example,  the  response  equation  of  the  hybrid 
design  (311A)  is 

y  =  .5243  (Ml)  +  1.195  (Bl)  +  .961  (SI) 

-  .00000251  (S1*S1)  -  .0000602  (B1*S1)  +  68.4118  (18) 

If  the  design  is  altered  so  that  the  first  and  third 
variables  in  the  design  are  switched,  then  the  new  design 
matrix  is  coded  as  follows;  the  SI  variable  is  coded  as 
variable  1;  the  Bl  variable  is  still  coded  as  variable  2; 
and  the  Ml  variable  is  coded  as  variable  3.  As  a  result, 
different  response  values  are  obtained  from  the  same  design 
matrix  because  different  uncoded  values  have  been  used  in 
the  MPOS  input  files.  The  resulting  response  equation  for 


the  altered  hybrid  design  is: 


y  =  .4678  (Ml)  +  1.11095  (Bl)  +  .98925  (SI) 

-  .0000275  (S1*S1)  -  .0000512  (B1*S1)  +  209.2299  (19) 

The  predictive  power  for  both  the  original  and  the  altered 
designs  are  equivalent.  That  is,  the  degrees  of  accuracy 
for  both  equations  are  the  same.  Howe'  er,  a  totally 
different  set  of  coefficients  has  been  generated  from  the 
altered  design  to  estimate  the  true  surface. 

A  second  reason  to  account  for  the  small  differences 
between  the  coefficients  is  the  small  amount  of 
mul ticol 1 ineari ty  which  exists  between  the  intercept  term 
and  the  squared  terms  and  between  the  squared  terms 
themselves.  Variables  which  possess  even  a  small  amount  of 
multicollinearity ,  produce  regression  coefficients  which 
tend  to  vary  from  one  design  to  the  next. 

For  example,  a  regression  equation  is  developed  from  a 
representative  sample  of  points  (design  1)  from  a  large 
population.  If  multicollinearity  exists  between  some  or  all 
of  the  variables,  a  second  sampling  (design  2)  from  the 
population  will  necessarily  produce  a  different  regression 
equation  because  of  the  multicollinearity  between  the 
variables.  The  degree  of  variability  between  the 
coefficients  of  the  two  equations  may  be  directly  related  to 
the  degree  of  multicollinearity  between  the  variables. 

Multicollinearity  can  be  eliminated  between  the 
variables  of  an  experimental  design,  if  the  design  is  made 


perfectly  orthogonal.  A  full  factorial  design  is  perfectly 
orthogonal,  if  the  required  transformations  are  performed  on 
the  squared  terms  of  the  design.  (The  required 
transformation  is:  Xj  “  Xj  .)  The  full  factorial  design  is 
not  commonly  used  in  RSM  because  it  requires  a  large  number 
of  design  points.  However,  this  design  serves  as  an 
excellent  standard  to  which  other  designs  can  be  compared. 

By  performing  the  necessary  transformations,  the  full 
factorial  three  factor  design  for  the  determinstic  model 
produces  the  following  response  equation: 

y  =  .4979  (Ml)  +  1.216  (Bl)  +  1.1978  (SI) 

-  .0(9003865  (S1«S1)  -  .0000649  (B1«S1) 

-  964.8424  (20) 

The  response  equation  of  the  full  factorial  design  achieves 
the  same  degree  of  accuracy  as  the  other  nearly  orthogonal 
designs  in  estimating  the  true  surface.  Thus,  the  primary 
difference  between  the  full  factorial  design  and  the  other 
nearly  orthogonal  designs  is  that  mul ticol 1 ineari ty  does  not 
exists  between  the  variables  in  the  full  factorial  design  so 
that  all  of  the  variables  are  uncorrelated. 

In  that  case,  the  full  factorial  design  provides  the 
most  representative  estimate  of  the  true  surface.  That  is, 
the  regression  equation  which  is  developed  from  the  full 
factorial  design  is  strictly  determined  by  the  number  and 
placement  of  its  design  points.  Whereas,  the  response 
equations  of  the  other  designs  are  also  determined  by  the 


number  and  placement  of  the  design  points,  but,  the  effects 
of  mul ticol 1 ineari ty  also  influence  the  coefficients  of  some 
of  the  variables  of  the  response  equations. 


Thus,  the  degree  of  variability  which  exists  between  the 
coefficients  of  the  full  factorial  response  equation  and  the 
coefficients  of  the  other  response  equations  may  be  a 
measure  of  the  degree  of  mul ticol 1 ineari ty  which  exists  in 
those  designs.  For  example,  the  coefficients  of  the 
response  equation  for  the  Box— Behnken  design  compare  more 
favorably  to  the  full  factorial  design  than  some  of  the 
other  three  factor  designs.  This  may  indicate  that  the 
influence  of  mul ticol 1 i nearity  for  three  factor  designs,  is 
not  as  significant  on  the  variables  of  the  Box-Behnken 
design  as  on  the  variables  of  the  other  designs.  However, 
this  statement  is  just  a  hypothesis  because  there  is  no 
measure  of  the  degree  of  orthogonality  for  the  various 
designs. 

In  short,  the  effects  of  mul ticol 1 ineari ty  on  response 
surface  equations  is  of  primary  importance  when  evaluating 
deterministic  models  because  the  analyst  is  attempting  to 
draw  direct  inferences  from  the  coefficients  of  the  response 
equation.  If  the  degree  of  mul ti col 1 i nearity  is  too  great 
between  the  variables  in  the  equation,  any  information 
obtained  from  the  coefficients  may  be  misleading. 

As  stated  earlier,  a  small  amount  of  mul ti col 1 ineari ty 
exists  in  all  of  the  equations  (except  the  full  factorial), 
therefore,  no  two  designs  produce  the  same  set  of 

3-21 


coef -f icients  when  estimating  a  surface.  Although  the 


differences  between  the  coefficients  are  small,  their 
effects  may  be  quite  significant. 

To  further  address  the  subject  matter  of  the  effects  of 
the  differences  between  the  coefficients,  the  designs  and 
the  response  equations  of  the  five  factor  model  are 
evaluated.  The  five  factor  model  is  selected  because  of  the 
simplicity  of  the  response  equations.  These  equations 
include  five  main  effect  terms  and  one  squared  term  in  its 
modified  form,  eq<16).  In  this  manner,  the  coefficients  of 
the  Ml,  Ii2,  B1  and  B2,  terms  should  represent  the  marginal 
effects  of  change  of  the  dependent  variable  when  the 
independent  variables  are  altered.  The  coefficient  of  the 
SI  term  is  excluded  from  the  list  because  the  effects  of 
change  caused  by  SI  are  influenced  by  the  Sl-squared  term. 
Therefore,  the  marginal  effects  of  SI  vary  for  different 
levels  of  SI.  Tables  11  through  14  compare  the  predicted 
change  to  the  actual  change  of  the  dependent  variable  when 
one  independent  variable  is  altered  while  all  other 
independent  variables  are  held  constant. 


TABLE  11 


ACTUAL  AND 
WHEN  Ml  IS 

PREDICTED  CHANGES 
i  DECREASED  BY  100 

CHANGE  IN  RESPONSE 

DESIGN 

ACTUAL 

PREDICTED 

PERCENT 

<1) 

Central 

Composite 

-57.0 

-50.  1 

12.  1 

(2> 

Box- 

Behnken 

-57.0 

—60. 3 

5.  S 

<3) 

Box- 

Draper  (H) 

-57.0 

-60.4 

6.0 

(4) 

Box- 

-57.0 

-62.9 

10.4 

Draper  (F) 


TABLE  12. 

ACTUAL  AND  PREDICTED  CHANCES 
WHEN  M2  IS  INCREASED  BY  400 


CHANGE  IN  RESPONSE 


DESIGN 

ACTUAL 

PREDICTED 

PERCENT 

<1) 

Central 
Composi te 

260.0 

251.9 

3.  1 

(2) 

Box- 

Behnken 

260.0 

271.8 

4.5 

(3) 

Box- 

Draper  (H) 

260.0 

273.0 

5.0 

(4) 

Box- 

Draper  (F) 

260.0 

271.0 

4.2 

TABLE  13 


ACTUAL  AND  PREDICTED  CHANGES 
WHEN  B1  IS  INCREASED  BY  200 


CHANGE  IN  RESPONSE 


DESIGN 

ACTUAL 

PREDICTED 

PERCENT 

(1) 

Central 
Composi te 

140.0 

142.5 

1.8 

(2) 

Box- 

Behnken 

140.0 

145.9 

4.2 

(3) 

Box- 

Draper  (H) 

140.0 

147.9 

5.6 

(4) 

Box- 

Draper  <F) 

140.0 

151. 1 

7.9 

TABLE  14. 

ACTUAL  AND  PREDICTED  CHANGES 
WHEN  B2  la  DECREASED  BY  500 

CHANGE  IN  RESPONSE 

DESIGN 

ACTUAL 

PREDICTED 

PERCENT 

(1) 

Central 
Composi te 

-370.0 

-377.25 

2.0 

(2) 

Box- 

Behnken 

-370.0 

-384.5 

3.9 

<3) 

Box- 

Draper  (H) 

-370.0 

-392.3 

6.0 

<4) 

Box- 

-370. 0 

-398.8 

7.8 

Draper  (F) 


ERROR 


ERROR 


There  are  two  primary  reasons  to  explain  the 
discrepancies  between  the  actual  and  the  predicted  changes 
of  the  various  response  surface  equations. 

The  first  reason  for  these  discrepancies  is  due  to 
mul ticol 1 ineari ty.  As  was  stated  earlier,  the  effects  of 
mul ti col 1 inear ity  influence  the  coefficients  of  the 
correlated  variables,  and  thereby  influence  their  ability 
to  accurately  measure  the  marginal  effects  of  change. 

The  second  reason  for  these  discrepancies  is  due  to  the 
fact  that  several  terms  have  been  omitted  from  the  full 
regression  model.  In  this  five  factor  example,  ten 
interaction  terms  and  four  squared  terms  have  been  removed. 
If  these  variables  are  inserted  into  eq(16),  the  differences 
between  the  predicted  change  and  the  actual  change  of  the 
dependent  variable  is  reduced.  The  full  regression  equation 
for  the  five  factor  model  is 

y  =  bg  +  b,<m)  +  b2(M2)  +  b3(Bl>  +  b^  (B2)  +  b5(Sl) 

+  b,2  (M1*M2)  +  b,3(Ml*Bl)  +  b,4  (M1*B2>  +  b,5(Ml*S2> 

+  bj3  <M2*B1)  +  b  24  <M2*B2)  +  bjg  (li2*Sl)  +  bj^  (B1*B2) 

+  bjg  (B1*S1)  +  b^g  (B2*S1)  +  b„  (M1*M1)  +  bj2(M2*M2) 

4-  b33  (B1*B1)  +  b^  <B2*B2>  +  b55(Sl*Sl)  (21) 

Tables  15  through  18  contain  the  same  information  as 
Tables  11  through  14,  except  for  the  fact  that  the  full 
regression  model  (5  main  effects,  10  interaction  terms,  5 
squared  terms,  and  1  intercept  term)  is  used  to  determine 


the  predicted  change  o-f  the  response  values  in  Tables  15 
through  18. 


TABLE  15. 

ACTUAL  AND  PREDICTED  CHANGES  WHEN  Ml 
IS  DECREASED  BY  100  (FULL  MODEL) 


CHANGE  IN  RESPONSE 


DESIGN 

ACTUAL 

PREDICTED 

PERCENT 

(1) 

Central 

Composite 

-57.0 

-57.5 

0.9 

(2) 

Box- 

Behnken 

-57.0 

-59.9 

5.  1 

(3) 

Box- 

Draper  (H) 

-57.0 

-60.2 

5.6 

(4) 

Box- 

Draper  (F) 

-57.0 

• 

S 

1 

5.8 

TABLE  16. 

ACTUAL  AND  PREDICTED  CHANGES  WHEN  M2 
IS  INCREASED  BY  400  (FULL  MODEL) 


CHANGE  IN  RESPONSE 


DESIGN 

ACTUAL 

PREDICTED 

PERCENT 

(1) 

Central 

Composite 

260.  0 

247.6 

4.8 

(2) 

Box- 

Behnken 

260.0 

267.5 

2.9 

(3) 

Box- 

Draper  (H) 

260.0 

272.5 

ffi 

• 

(4) 

Box- 

260.0 

269.8 

3.8 

Draper  (F) 


TABLE  17 


ACTUAL  AND  PREDICTED  CHANGES  WHEN  B1 
IS  INCREASED  BY  200  <FULL  MODEL) 


CHANGE  IN  RESPONSE 


DESIGN 

ACTUAL 

PREDICTED 

PERCENT 

(1) 

Central 

Composite 

140.0 

140.  1 

0.  1 

(2) 

Box- 

Behnken 

140.0 

145.  1 

3.6 

(3) 

Box- 

Draper  (H) 

140.0 

146.9 

4.9 

(4) 

Box- 

Draper  <F) 

140.0 

148.2 

5.8 

TABLE  18. 

ACTUAL  AND  PREDICTED  CHANGES  WHEN  B2 
IS  DECREASED  BY  500  (FULL  MODEL) 


CHANGE  IN  RESPONSE 


DESIGN 

ACTUAL 

PREDICTED 

PERCENT 

(1) 

Central 

Composite 

-370.0 

-364. 1 

1.6 

(2) 

Box- 

Behnken 

-370.0 

-377.9 

2.  1 

(3) 

Box- 

Draper  <H) 

-370.0 

-387. 1 

4.6 

(4) 

Box- 

Draper  (F) 

-370.0 

-389.5 

5.3 

From  the  -Full  regression  models,  the  analyst  is  able  to 
obtain  very  e-f-fective  measures  of  the  marginal  effects  of 


the  independent  variables.  The  procedure  to  obtain  this 
measure  consists  of  calculating  the  response  value  from  the 
response  equation  twice.  During  the  first  calculation,  all 
variables  are  set  to  a  base  level;  during  the  second 
calculation,  all  variables  remain  fixed,  except  for  the 
variable  whose  marginal  effect  is  to  be  measured. 
Subsequently,  the  difference  in  the  response  values  of  the 
two  calculations  represents  the  marginal  effect  of  that 
variable. 

If  the  full  regression  model  is  used  to  measure  the 
marginal  effects  of  an  independent  variable,  then  the 
requirement  for  low  multicollinearity  between  the  variables 
is  eliminated.  The  analyst  does  not  obtain  the  marginal 
effects  information  from  the  coefficients,  but  from  the 
difference  between  the  estimated  response  values.  Since 
multicollinearity  does  not  affect  a  design's  ability  to 
produce  a  response  equation  that  predicts  effectively,  the 
response  equation  of  a  non— orthogonal  design  is  able  to 
measure  the  marginal  effects  of  an  independent  variable  as 
well  as  the  equation  of  an  orthogonal  design. 

The  drawback  to  using  the  full  regression  model  is  the 
complexity  of  eq<21).  The  ability  to  measure  the  marginal 
effects  is  no  longer  a  simple  proposal  and  the  marginal 
effects  of  an  independent  variable  is  dependent  on  the  level 
at  which  it  is  evaluated.  Essentially,  the  analyst  performs 
a  one  variable  at  a  time  procedure  to  obtain  his 


information 


Therefore,  the  method  involving  the  modified 


regression  model  may  be  a  more  attractive  alternative 
because  o-F  its  simplicity,  provided  that  the  decision  maker 
is  able  to  accept  some  limitations  in  the  magnitudes  o-F  the 
estimates  of  the  predicted  marginal  effects. 

Anal.yti.cal.  Val.ue  of  the  t3Qdi.fi.ed  Response 
Equations.  The  question  of  whether  or  not  any  error  level 
is  too  great  is  dependent  on  the  degree  of  accuracy  which 
the  decision  maker  requires.  An  important  point  which  the 
decision  maker  must  consider,  is  that  the  information 
pertaining  to  the  marginal  effects  of  the  independent 
variables  is  readily  available  from  the  simplified  response 
equation,  but  not  easily  obtainable  from  other  sources. 

In  examiming  the  predicted  marginal  effects  of  the 
modified  equations.  Tables  11  through  14  indicate  that  there 
are  discrepancies  between  the  predicted  and  the  actual 
changes  in  the  response  values.  The  range  of  the  error 
between  the  predicted  and  actual  change  is  l.BV.  to  12.1%. 

In  general,  an  error  of  less  than  27.  is  likely  to  be 
accepted  in  just  about  any  circumstance.  However,  an  error 
of  12%  may  not  be  acceptable  in  certain  situations. 

A  simple  test  to  check  the  magnitude  of  the  error  of 
the  predicted  change  for  the  coefficients  of  the  modified 
response  equation,  is  to  compare  the  predicted  change  of  the 
modified  response  equation  to  the  predicted  change  of  the 
full  regression  equation.  If  the  level  of  error  is 
acceptable  to  the  decision  maker,  the  modified  response 
equation  is  used  and  is  a  very  useful  analytical  tool. 


However,  if  the  level  of  error  is  not  acceptable,  the 
analyst  may  choose  to  use  the  full  regression  equation  to 
measure  the  marginal  effects  of  change  for  the  independent 
variables.  But  as  was  noted  earlier,  this  procedure  may  be 
rather  lengthy,  in  that  the  marginal  effects  of  change  of 
any  variable  is  not  constant  throughout  the  entire  surface. 
Therefore,  the  variables  must  be  evaluated  at  several 
different  levels  in  order  to  determine  the  marginal  effect. 

An  alternative  to  using  the  full  regression  model  when 
the  level  of  error  of  the  modified  equation  is  unacceptable, 
is  to  evaluate  the  ratios  of  the  coefficients  of  the 
modified  equations.  In  regression  analysis,  the  ratio  of 
the  coefficients  indicates  the  contribution  of  one  factor  as 
compared  to  another  factor.  Therefore,  although  the  exact 
magnitude  of  their  impact  on  the  response  value  can  not  be 
determined,  the  relative  importance  of  the  various 
independent  variables  can  be  assessed  by  rank  ordering  the 
ratios  of  the  coefficients. 

For  example,  the  modified  response  equation  for  the 
five  factor  Box-Behnken  design  is 

y  =  .6029  (Ml)  +  .6795  (M2)  +  .7294  (Bl) 

+  .7689  (B2)  +  1.0363  (SI) 

-  .0000295  (S1*S1)  +  843.9903  (22) 

The  mean  contribution  of  B2  to  the  response  value,  y,  is 
estimated  to  be  1.27  times  as  great  as  the  mean  contribution 


o-f  Ml  to  the  response  value.  This  amount  is  obtained  by 
calculating  the  ratio  of  the  coefficients  for  the  two 
independent  variables  (. 7689/ . 6029) . 

In  a  similar  manner,  eq<23)  is  the  response  equation 
for  the  five  factor  Central  Composite  Rotatable  design. 
According  to  this  response  equation,  the  mean  contribution 
of  B2  is  1.5  times  as  great  as  the  mean  contribution  of  Ml 
to  the  response  value  (. 7545/. 501 1 ) . 

y'  =  .5011  <M1)  +  .6289  (M2)  +  .7123  (Bl) 

+  .7545  (B2)  +  .9511  (SI) 

-  .0000261  (S1*S1)  +  1615.6741  (23) 

Since  the  various  designs  which  are  used  in  this  study 
produce  different  coefficients  for  each  design,  the 
ratios  of  coefficients  are  different.  Hence,  the  exact 
magnitude  of  the  contribution  of  each  factor  as  compared  to 
other  factors  can  not  be  determined  with  certainty. 
Nevertheless,  the  relative  contribution  of  the  factors  can 
still  be  obtained  and  these  contributions  can  be  rank 
ordered. 

Therefore,  an  accurate,  but  less  definitive  statement 
about  the  relative  contribution  of  the  two  independent 
variables,  is  that  the  mean  contribution  of  B2  is  greater 
than  the  mean  contribution  of  Ml  to  the  response  value. 

This  conclusion  is  permissible  because  the  relative 
magnitudes  of  the  coefficients  are  maintained  for  the 
various  designs  even  though  the  values  of  the  coefficients 


themselves  change  from  one  design  to  another.  That  is,  no 
matter  what  orthogonal  design  is  Selected,  the  regression 
coefficient  for  B2  always  exceeds  the  coefficient  for  Ml. 

In  this  manner,  the  following  conclusion  is  drawn  for  the 
five  factor  problem:  the  independent  variable  which  causes 
the  greatest  change  is  SI,  provided  the  values  for  Si  remain 
within  the  relevant  range.  This  variable  is  followed  in 
order  of  importance  by:  B2,  B1 ,  M2,  and  Ml.  Although  this 
conclusion  does  not  provide  an  exact  measure  of  the 
difference  between  the  factors,  the  relative  ranking  of  the 
factors  is  still  an  important  piece  of  information  which  the 
decision  maker  needs  in  order  to  evaluate  his  alternatives. 


IV.  CONCLUSION 


overview 

This  research  identifies  several  experimental  designs 
which  are  effective  and  efficient  in  estimating  response 
surface  equations  for  deterministic  models.  The  specific 
deterministic  model  which  is  used  in  this  study  is  a  nuclear 
exchange  linear  programming  model,  but  the  results  which  are 
obtained  in  this  study  apply  to  any  deterministic  model. 

The  purpose  of  Response  Surface  Methodlogy  in  its 
application  to  any  type  of  model  is  to  obtain  a  response 
surface  equation  which  closely  approximates  the  true 
surface.  A  second  order  polynomial  has  been  shown  to  be  an 
excellent  choice  for  estimating  response  surfaces  of 
determininstic ,  nuclear  exchange  models.  Ideally,  the 
experimental  design  should  use  as  few  points  as  possible  to 
obtain  an  accurate  fit  and  the  response  equation  should 
contain  as  few  terms  as  possible.  That  is,  the  equation 
should  contain  as  few  cross  product  and  squared  terms  as 
possible  so  that  direct  inferences  can  be  drawn  from  the 
coefficients  of  the  main  effects. 

iPECIFIC  RESULIi 

In  keeping  with  the  idea  of  simplifying  the  response 
equation,  the  response  surface  equations  which  are  used  to 
evaluate  error  in  fit  of  the  design  are  modified  versions  of 
the  full  regression  equations.  That  is,  selected  terms  from 


the  -Full  regression  equations  have  been  omitted  from  the 
modified  regression  equations  because  these  terms  account 
for  less  than  two  percent  of  the  total  sum  of  squares. 

In  evaluating  the  response  equations  of  these  designs 
for  their  effectiveness,  the  results  indicate  that  the 
accuracy  of  a  modified  response  equation  is  dependent  on  the 
number  of  design  points  and  on  the  degree  of  rotatability  of 
the  design.  The  accuracy  of  a  response  surface  equation  is 
highly  correlated  to  the  number  of  design  points,  but  this 
relationship  is  not  a  perfect  correlation,  and  their 
differences  can  be  explained  by  the  degree  of  rotatability 
of  the  designs. 

In  short,  more  design  points  improves  accuracy, 
provided  the  points  are  spread  throughout  the  entire 
surface.  Similarly,  if  the  region  of  interest  is  spherical, 
a  rotatable  design  improves  the  accuracy  of  the  response 
equation. 

In  evaluating  these  designs  for  their  efficiency  as 
well  as  effectiveness,  the  results  indicate  that  a  design 
which  provides  the  most  accurate  fit  with  the  fewest  number 
of  design  points  for  all  factors  does  not  exist.  In  fact, 
minimizing  the  number  of  design  points  and  maximizing 
accuracy  are  dichotomous  goals.  This  trade-off  between 
accuracy  and  required  number  of  design  points  is  dependent 
on  the  problem  at  hand  and  needs  to  be  defined  seperately 
for  each  situation.  Therefore,  it  remains  for  the  decision 
maker  to  evaluate  the  problem  to  determine  which  is  more 
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important:  a  more  accurate  fit  or  fewer  design  points. 

In  evaluating  the  explanatory  power  of  the  response 
equations,  the  coefficients  of  the  modified  response 
equation  represent  marginal  effects  of  change  of  the 
dependent  variable  when  the  independent  variables  are 
altered  provided  specific  conditions  are  met.  These 
conditions  are:  a  high  degree  of  multi  col linearity  does  not 
exist  between  the  independent  variables;  the  main  effect 
terms  of  the  response  equation  dominate  the  squared  and  the 
interaction  terms;  and  the  modified  response  surface  is  able 
to  accurately  represent  the  true  surface. 

The  results  of  this  project  indicate  that  there  are 
limitations  in  the  ability  of  the  coefficients  of  the 
modified  response  equation  to  predict  the  marginal  effects 
of  change.  The  accuracy  of  these  models  to  predict  the 
marginal  change  can  be  evaluated  by  comparing  the  predicted 
changes  of  the  coefficients  of  the  modified  response 
equations  to  the  coefficients  of  the  full  regression 
equations.  If  this  difference  is  unacceptable,  then  the 
full  regression  equation  can  be  used  to  measure  the  marginal 
effects  or  the  ratios  of  the  coefficients  of  the  modified 
equation  can  be  used  to  determine  the  relative  importance  of 
the  independent  variables. 

RECOMMENDATIONS 

This  study  originally  intended  to  examine  the 
differences  between  bias  minimizing  and  variance 


minimirzing  designs.  As  this  study  evolved,  it  became 
clear  that  there  were  no  significant  differences  between 
the  two  types  of  experimental  designs,  provided  that  the 
designs  maintained  certain  characteristics,  such  as, 
rotatability  and  orthogonality. 

However,  very  few  designs  possess  pure  rotatability 
or  orthogonality.  Most  of  the  designs  are  classified  in 
a  gray  area  called,  nearly  rotatable  and  nearly 
orthogonal.  As  such,  further  research  is  required  in 
specifying  nearly  rotatable  and  nearly  orthogonal.  This 
research  may  find  what  levels  of  rotatability  are 
necessary  to  achieve  a  desired  level  of  predictive  power. 
Also,  this  research  may  find  what  levels  of  orthogonality 
are  necessary  to  achieve  a  desired  level  of  explanatory 
power . 


Ai  MPOS  iDEBt  EllSB 


This  appendix  describes  the  speci-Fic  problems  which  were 
analyzed  for  the  three,  four,  five,  and  six  factor  cases.  In 
addition,  this  appendix  lists  a  sample  of  the  Multi-Purpose 
Optimization  System  <MPOS)  input  files  which  were  used  to 
obtain  the  response  value. 

three  faciqr  problem 

There  are  three  weapon  types  and  five  target  classes. 

The  three  weapon  types  are: 


WEAPON  TYPE 

LOW  VALUE 

high  VALUg 

Ml 

1650 

3000 

B1 

1200 

3000 

SI 

6000 

12000 

Assumptions  Each  system's  reliability,  probability  of  launch 
survivability,  probability  of  arrival,  and  warhead 
reliability  are  all  1.0.  The  single  shot  probability  of  kill 
(which  is  damage  expectancy  on  a  target  because  of  the 
previous  assumptions  of  reliability,  etc.)  for  each  warhead 
type  is  provided  in  the  following  tables 


TABLE  A. 1 


SINGLE  SHOT  PROBABILITY  OF  KILL 
FOR  THREE  FACTOR  PROBLEM 


TARGET  TYPE 

WEAPON 

lYPE 

Ml 

B1 

SI 

TGT 

CLASS 

1 

.84 

.96 

.73 

TGT 

CLASS 

2 

.75 

.88 

.64 

TGT 

CLASS 

3 

.55 

.70 

.32 

TGT 

CLASS 

4 

.88 

.98 

.78 

TGT 

CLASS 

5 

.25 

.48 

.  15 

Target  class  4  represents  time  urgent  targets  and  only 
ICBMs  (Ml)  and  SLBMs  (Si)  are  allowed  to  hit  them.  In 
addition,  the  maximum  allowable  damage  on  target  class  2 
targets  is  .95  of  the  total  target  value  and  at  least  .60  of 
the  total  target  value  of  target  class  5  targets  must  be 
destroyed.  Furthermore,  a  maximum  of  two  warheads  can 
targeted  at  any  individual  target. 

The  number  of  targets  in  each  class  is  6000,  4000, 
2000,  1000,  and  200,  respectively.  The  measure  of 
effectiveness  for  this  problem  is  the  total  damage 
expectancy,  i.e. ,  the  sum  of  the  product  of  the  damage 
expectancy  on  a  target  times  the  number  of  targets. 


FOUR  FACTOR  PROBLEM 


There  are 

four  weapon  types  and 

five  target  classes 

The  four  weapon 

types  are: 

WEAPON  TYPE 

LQW  VALUE 

high  value 

Ml 

1650 

3000 

M2 

500 

2000 

B1 

1200 

3000 

SI 

6000 

12000 

Assumption:  Each  system's  reliability,  probability  of  launch 
survivability,  probability  of  arrival,  and  warhead 
reliability  are  all  1.0.  The  single  shot  probability  of  kill 
(which  is  damage  expectancy  on  a  target  because  of  the 
previous  assumptions  of  reliability,  etc.)  for  each  warhead 
type  is  provided  in  the  following  table: 


TABLE  A. 2 


« 
* 


SINGLE  SHOT  FHOBABILITY  OF  KILL 


FOR  FOUR 

FACTOR  PROBLEM 

TARGET 

TYPE 

WEAPON  TYPE 

Ml 

M2 

B1 

SI 

TGT  CLASS 

1 

.84 

.90 

.96 

.73 

TGT  CLASS 

2 

.75 

.83 

.88 

.64 

TGT  CLASS 

3 

.55 

.62 

.70 

.32 

TGT  CLASS 

4 

.88 

.94 

.98 

.78 

TGT  CLASS 

5 

.25 

.36 

.48 

.  15 

Target  class  4  represents  time  urgent  targets  and  only 
ICBMs  (Ml,  M2)  and  SLBMs  (SI)  are  allowed  to  hit  them.  In 
addition,  the  maximum  allowable  damage  on  target  class  2 
targets  is  .95  o-f  the  total  target  value  and  at  least  .60  of 
the  total  target  value  of  target  class  5  targets  must  be 
destroyed.  Furthermore,  a  maximum  of  two  warheads  can 
targeted  at  any  individual  target. 

The  number  of  targets  in  each  class  is  6000,  4000, 
2000,  1000,  and  200,  respectively.  The  measure  of 
effectiveness  for  this  problem  is  the  total  damage 
expectancy,  i.e. ,  the  sum  of  the  product  of  the  damage 
expectancy  on  a  target  times  the  number  of  targets. 
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FIVE  FACTOR  PROBLEM 


There  are  five  weapon  types  and  five  target  classes. 
The  five  weapon  types  are: 


WEAPON  TYPE 

LOW  VALUE 

HISH  VALUE 

Ml 

1650 

3000 

M2 

500 

2000 

B1 

1200 

3000 

B2 

2000 

4000 

SI 

6000 

12000 

Assumption:  Each  system's  reliability,  probability  of  launch 
survivability,  probability  of  arrival,  and  warhead 
reliability  are  all  1.0.  The  single  shot  probability  of  kill 
(which  is  damage  expectancy  on  a  target  because  of  the 
previous  assumptions  of  reliability,  etc.)  for  each  warhead 
type  is  provided  in  the  following  table: 
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TABLE  A. 3 


SINGLE  SHOT  PROBABILITY  OF  KILL 
FOR  FIVE  FACTOR  PROBLEM 

IARGEI  iype  weapon  iype 


Ml 

M2 

B1 

B2 

SI 

TGT 

CLASS 

1 

.84 

.90 

.96 

.98 

.73 

TGT 

CLASS 

2 

.75 

.83 

.88 

.92 

.64 

TGT 

CLASS 

3 

.55 

.62 

.70 

.74 

.32 

TGT 

CLASS 

4 

.88 

.94 

.98 

.99 

.78 

TGT 

CLASS 

5 

.25 

.36 

.48 

.55 

.  15 

Target  class  4  represents  time  urgent  targets  and  only 
ICBMs  (Ml,  M2)  and  SLBMs  (SI)  are  allowed  to  hit  them.  In 
addition,  the  maximum  allowable  damage  on  target  class  2 
targets  is  .95  o-f  the  total  target  value  and  at  least  .60  of 
the  total  target  value  of  target  class  5  targets  must  be 
destroyed.  Furthermore,  a  maximum  of  two  warheads  can 
targeted  at  any  individual  target. 

The  number  of  targets  in  each  class  is  9000,  6000,  3000, 
1500,  and  300,  respectively.  The  measure  of  effectiveness 
for  this  problem  is  the  total  damage  expectancy,  i.e. ,  the 
sum  of  the  product  of  the  damage  expectancy  on  a  target 
times  the  number  of  targets. 
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SIX  FACTOR  PROBLEM 


There  are  six  weapon 

types  and 

five  target  classes 

six  weapon  types  are: 

WEAPON  TYPE 

LOW  VALUE 

HIGH  VALUE 

Ml 

1650 

3000 

M2 

500 

2000 

M3 

400 

1000 

B1 

1200 

3000 

B2 

2000 

4000 

SI 

6000 

12000 

Assumption:  Each  system's  reliability,  probability  o-f  launch 
survivability,  probability  of  arrival,  and  warhead 
reliability  are  all  1.0.  the  single  shot  probability  of  kill 
(which  is  damage  expectancy  on  a  target  because  of  the 
previous  assumptions  of  reliability,  etc.)  for  each  warhead 


type  is 


TABLE  A. 4 


»2»‘l 

€ 
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A 

15 


SINGLE  SHOT  PROBABILITY  OF  KILL 
FOR  SIX  FACTOR  PROBLEM 


TARGET  TYPE 


WEAPON  TYPE 


Ml 

M2 

M3 

B1 

B2 

SI 

TGT 

CLASS 

1 

.84 

.90 

.92 

.96 

.98 

.73 

TGT 

CLASS 

2 

.75 

.83 

.87 

.88 

.92 

.64 

TGT 

CLASS 

3 

.55 

.62 

.65 

.70 

.74 

.32 

TGT 

CLASS 

4 

.88 

.94 

.95 

.98 

.99 

.78 

TGT 

CLASS 

5 

.25 

• 

.45 

.48 

.55 

.  15 

Target  class  4  represents  time  urgent  targets  and  only 
ICBMs  (Ml,  M2,  M3)  and  SLBMs  (SI)  are  allowed  to  hit  them. 

In  addition,  the  maximum  allowable  damage  on  target  class  2 
targets  is  ,95  o-f  the  total  target  value  and  at  least  .60  of 
the  total  target  value  of  target  class  5  targets  must  be 
destroyed.  Furthermore,  a  maximum  of  two  warheads  can 
targeted  at  any  individual  target. 

The  number  of  targets  in  each  class  is  9000,  6000,  3000, 
1500,  and  300,  respectively.  The  measure  of  effectiveness 
for  this  problem  is  the  total  damage  expectancy,  i.e.  ,  the 
sum  of  the  product  of  the  damage  expectancy  on  a  target 
times  the  number  of  targets. 
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m 

m 

******************************************** 
************  MPOS  INPUT  FILE  ************** 
************  3  FACTORS  ************** 

******************************************* 


REGULAR 

TITLE 

RESPONSE  SURFACE  FOR  ARSENAL  EXCHANGE  MODEL 


VARIABLES 

** 

WEAPON 

** 

** 

SYSTEM 

** 

Xll 

TO 

X15 

** 

Ml 

** 

X21 

TO 

X25 

** 

B1 

** 

X31 

TO 

X35 

** 

SI 

** 

X41 

TO 

X45 

** 

Ml 

+ 

Ml 

** 

X51 

TO 

X55 

** 

B1 

+ 

B1 

** 

X61 

TO 

X65 

** 

SI 

+ 

SI 

** 

X71 

TO 

X75 

** 

Ml 

+ 

B1 

** 

X81 

TO 

X85 

** 

Ml 

+ 

SI 

** 

X91 

TO 

X95 

** 

B1 

+ 

SI 

** 

MAXIMIZE 

** 

DE 

TABLE  ** 

.8400 

Xll 

+ 

.7500 

X12 

.5500 

X13 

+ 

.8800 

X14 

+ 

.2500 

X15 

+ 

.9600 

X21 

+ 

.8800 

X22 

*  . 7000 

X23 

+ 

.  0000 

X24 

+ 

.4800 

X25 

+ 

.7300 

X31 

+ 

.6400 

X32 

+  .3200 

X33 

+ 

.7800 

X34 

+ 

.  1500 

X35 

.9744 

X41 

+ 

.9375 

X42 

+  .7975 

X43 

+ 

.9856 

X44 

+ 

.4375 

X45 

+ 

.9984 

X51 

+ 

.9856 

X52 

+  .9100 

X53 

+ 

.0000 

X54 

+ 

.7296 

X55 

+ 

.9271 

X61 

+ 

.8704 

X62 

+  .5376 

X63 

+ 

.9516 

X64 

+ 

.2775 

X65 

+ 

.9936 

X71 

+ 

.9700 

X72 

+  .8650 

X73 

+ 

.0000 

X74 

+ 

.6100 

X75 

+ 

.9568 

X81 

+ 

.9100 

X82 

+  .6940 

X83 

+ 

.9736 

X84 

+ 

.3625 

X85 

+ 

.9892 

X91 

.9568 

X92 

+  .7960 

X93 

+ 

.  0000 

X94 

+ 

.5580 

X95 

CONSTRAINTS 

**  <  1) 

95  '/ 

L  DE  CRITERIA 

** 

.7500 

X12 

+ 

.8800 

X22 

+  .6400 

X32 

+ 

.9375 

X42 

+ 

.9856 

X52 

.8704 

X62 

+ 

.9700 

X72 

+  .9100 

X82 

+ 

.9568 

X92 

.LE.  3800 

**  (2) 

60  V 

:  DE  CRITERIA 

** 

.2500 

X15 

+ 

.4800 

X25 

+  . 1500 

X35 

+ 

.4375 

X45 

+ 

.7296 

X55 

.2775 

X6S 

+ 

.6100 

X75 

+  .3625 

X85 

+ 

.5580 

X95 

.LE.  120 

**  (3)  TARGET  1  CONSTRAINT  ** 

Xll  +  X21  +  X31  +  X41  +  X51  +  X61  +  X71  + 

X81  +  X91  .LE.  6000 

**  (4)  TARGET  2  CONSTRAINT  ** 

X12  +  X22  +  X32  +  X42  +  X52  +  X62  +  X72  + 

X82  +  X92  .LE.  4000 


**  (5)  TARGET  3  CONSTRAINT  ** 

X13  +  X23  +  X33  +  X43  +  X53  +  X63  +  X73  + 

X83  +  X93  .LE.  2000 


**  (6)  TARGET  4  CONSTRAINT  ** 
X14  +  X24  +  X34  +  X44  +  X54 

X84  +  X94  .LE.  1000 


+  X64  +  X74  + 


**  (7)  TARGET  5  CONSTRAINTS  ** 

X15  +  X25  +  X35  +  X45  +  X55  +  X65  +  ) 

X85  +  X95  .LE.  200 

**  <8>  WEAPON  SYSTEM  1  <M1)  CONSTRAINT  ** 
Xll  +  X12  +  X13  +  X14  +  X15  + 

2  X41  +  2  X42  +  2  X43  +  2  X44  +  2  X45  + 

X71  +  X72  +  X73  +  X74  +  X75  + 

X81  +  X82  +  X83  +  X84  +  X85  .LE.  3000 

**  (9)  WEAPON  SYSTEM  2  <B1)  CONSTRAINT 

X21  +  X22  +  X23  +  X24  +  X25  + 

2  X51  +  2  X52  +  2  X53  +  2  X54  +  2  X55  + 

X71  +  X72  +  X73  +  X74  +  X75  + 

X91  +  X92  +  X93  +  X94  +  X95  .LE.  3000 

**  (10)  WEAPON  SYSTEM  3  (SI)  CONSTRAINT  ** 
X31  +  X32  +  X33  +  X34  +  X35  + 

2  X61  +  2  X62  +  2  X63  +  2  X64  +  2  X65  + 

X81  +  X82  +  X83  +  X84  +•  X85  + 

X91  +  X92  +  X93  +  X94  +  X95  .LE.  9000 

PRINT 

OPTIMIZE 


A  -  10 


«■»*•»*•»«*«**«'  MPQS  INPUT  FILE  ************** 
************  4  FACTOR  ************** 

******************************************* 

REGULAR 

TITLE 

RESPONSE  SURFACE  FOR  ARSENAL  EXCHANGE  MODEL 


VARIABLES 

** 

WEAPON 

** 

** 

SYSTEM 

** 

Xll 

TO 

X15 

** 

Ml 

** 

X21 

TO 

X25 

** 

B1 

** 

X31 

TO 

X35 

** 

SI 

** 

X41 

TO 

X45 

** 

Ml 

+ 

Ml 

** 

X51 

TO 

X55 

** 

B1 

B1 

** 

X61 

TO 

X65 

** 

SI 

+ 

SI 

** 

X71 

TO 

X75 

** 

Ml 

+ 

B1 

** 

X81 

TO 

X85 

** 

Ml 

+ 

SI 

** 

X91 

TO 

X95 

** 

B1 

+ 

SI 

** 

X101 

TO 

X185 

** 

M2 

** 

X1J.1 

TO 

X115 

** 

M2 

+ 

M2 

** 

X12J 

TO 

X125 

** 

Ml 

+ 

M2 

** 

X131 

TO 

X135 

** 

M2 

+ 

B1 

** 

X141 

TO 

X145 

** 

M2 

+ 

SI 

** 

MAXIMIZE 

** 

DE 

TABLE  ** 

.8400 

Xll 

+ 

.7500 

X12 

.  5500 

X13 

+ 

.8800 

X14 

+ 

.2500 

X15  + 

.9600 

X21 

+ 

.8800 

X22 

+  .7000 

X23 

+ 

.0000 

X24 

+ 

.4800 

X25  + 

.7300 

X31 

+ 

.6400 

X32 

+  .3200 

X33 

+ 

.7800 

X34 

+ 

.  1500 

X35  + 

.9744 

X41 

+ 

.9375 

X42 

+  .7975 

X43 

+ 

.9856 

X44 

+ 

.4375 

X45  + 

.9984 

X51 

+ 

.9856 

X52 

+  .9100 

X53 

+ 

.0000 

X54 

.7296 

X55  + 

.9271 

X61 

+ 

.8704 

X62 

+  .5376 

X63 

+ 

.9516 

X64 

+ 

.2775 

X65  + 

.9936 

X71 

.9700 

X72 

+  .8650 

X73 

+ 

.  0000 

X74 

+ 

.6100 

X75  + 

.9568 

X81 

+ 

.9100 

X82 

+  .6940 

X83 

+ 

.9736 

X84 

+ 

.3625 

X85  + 

.9892 

X91 

.9568 

X92 

+  .7960 

X93 

+ 

.  0000 

X94 

+ 

.5580 

X95  + 

.9000 

X101 

+ 

.  8300 

X102 

.6200 

X103 

+ 

.9400 

X104 

+ 

.  3600 

X105  + 

.9900 

Xlll 

+ 

.9711 

X112 

+  .8556 

X113 

+ 

.9964 

X114 

+ 

.5904 

X115  + 

.9840 

X121 

+ 

.9575 

X122 

+  .8290 

X123 

+ 

.9928 

X124 

+ 

.5200 

X125  + 

.9960 

X131 

+ 

.9796 

X132 

+  .8860 

X133 

+ 

.  0000 

X134 

+ 

.6672 

X135  + 

.9730 

X141 

+ 

.9388 

X142 

+  .7416 

X143 

+ 

.9868 

X144 

.4560 

X145 

CONSTRAINTS 

**  (1) 

95  : 

i  DE  CRITERIA 

** 

.7500 

X12 

.8800 

X22 

+  . 6400 

X32 

+ 

.9375 

X42 

+ 

.9856 

X52  + 

.8704 

X62 

+ 

.9700 

X72 

+  .9100 

X82 

+ 

.9568 

X92 

+ 

.8300 

X102  + 

.9711 

X112 

+ 

.9575 

X122 

+  .9796 

X132 

+ 

.9388 

X142 

.LE.  3800 

**  (2) 

60  */.  DE 

:  CRITERIA  ** 

.2500  X15  +  . 

4800  X25  +  .1500  X35  +  .4375 

X45 

+  . 

7296  X55  + 

.2775  X65  +  . 

6100  X75  +  .3625  X85  +  .5580 

X95 

+  . 

3600  XI 05  + 

.5904  XI 15  +  . 

5200  X125  +  .6672  X135  +  .4560 

X145 

.GE 

.  120 

♦*  (3) 

TARGET 

1  CONSTRAINT  ** 

Xll  + 

X21  + 

X31  +  X41  +  X51 

+  X61  +  X71 

+ 

X81 

+ 

X91  + 

X101  + 

Xlll  +  X121  +  X131 

+  X141  .LE. 

6000 

**  (4) 

TARGET 

2  CONSTRAINT  *♦ 

X12  + 

X22  + 

X32  +  X42  +  X52 

+  X62  +  X72 

+ 

X82 

+ 

X92  + 

X102  + 

X112  +  X122  +  X132 

+  X142  .LE. 

4000 

**  (5) 

TARGET 

3  CONSTRAINT  ** 

X13  + 

X23  + 

X33  +  X43  +  X53 

+  X63  +  X73 

+ 

X83 

+ 

X93  + 

X103  + 

X113  +  X123  +  X133 

+  X143  .LE. 

2000 

*•»  <6) 

TARGET 

4  CONSTRAINT  ** 

X14  + 

X24  + 

X34  +  X44  +  X54 

+  X64  +  X74 

+ 

X84 

+ 

X94  + 

X104  + 

X114  +  X124  +  X134 

+  X144  .LE. 

1000 

*♦  (7) 

TARGET 

5  CONSTRAINTS  *♦ 

X15  + 

X25  + 

X35  +  X45  X55 

+  X65  +  X75 

1  + 

X85 

+ 

X95  + 

X105  + 

X115  +  X125  +  X135 

+  X145  .LE. 

200 

##  (8) 

WEAPON 

SYSTEM  1  (Ml)  CONSTRAINT  ** 

Xll  + 

X12  +  X13  +  X14  +  X15 

+ 

2  X41 

+  2  X42  +  2  X43  +  2  X44  + 

2  X45  + 

X71  + 

X72  +  X73  +  X74  +  X75  + 

X81  + 

X82  +  XB3  +  X84  +  X85  + 

X121  + 

X122  + 

X123  +  X124  +  X125 

.LE.  3000 

**  (9) 

WEAPON 

SYSTEM  2  (M2)  CONSTRAINT  ** 

X101  -i- 

X102  + 

X103  +  X104  +  X105 

+ 

2  Xlll 

+  2  X112  +  2  X113  +  2  X114  +  2  X115  + 

X121  + 

X122  + 

X123  +  X124  +  X125 

+ 

X131  + 

X132  + 

X133  +  X134  +  X135 

+ 

X141  + 

X142  + 

X143  +  X144  +  X145 

.LE.  2000 

(10) 

'  WEAPON  SYSTEM  3  (Bl)  CONSTRAINT 

X21  + 

X22  + 

X23  +  X24  +  X25 

+ 

2  X51  +  2  X52 

+  2  X53  +  2  X54  +  I 

2  X55  + 

X71  + 

X72  + 

X73  +  X74  +  X75 

+ 

X91  + 

X92  + 

X93  +  X94  +  X95 

+ 

X131  + 

X132  + 

X133  +  X134  +  X135 

.LE.  3000 

A  -  12 

**  (11)  WEAPON  SYSTEM  4  <S1)  CONSTRAINT  ** 
X31  +  X32  +  X33  +  X34  +  X35  + 

2  X61  +  2  X62  +  2  X63  +  2  X64  +  2  X65  + 

X81  +  X82  +  X83  +  X84  +  X85  + 

X91  +  X92  +  X93  +  X94  +  X95  + 

X141  +  X142  +  X143  +  X144  +  X145  .LE.  9000 

PRINT 

OPTIMIZE 


MPQS  INPUT  PILE  ************** 
************  5  FACTORS  ************** 
******************************************* 


REGULAR 

TITLE 

RESPONSE  SURFACE  FOR  ARSENAL  EXCHANGE  MODEL 


VARIABLES 

** 

WEAPON 

*« 

** 

SYSTEM 

*« 

Xll 

TO 

X15 

** 

Ml 

** 

X21 

TO 

X25 

** 

B1 

*« 

X31 

TO 

X35 

** 

SI 

** 

X41 

TO 

X45 

** 

Ml 

+ 

Ml 

*« 

X51 

TO 

X55 

** 

B1 

B1 

** 

X61 

TO 

X6S 

** 

SI 

+ 

SI 

** 

X71 

TO 

X75 

** 

Ml 

+ 

B1 

*« 

X81 

TO 

X85 

** 

Ml 

+ 

SI 

** 

X91 

TO 

X95 

** 

B1 

+ 

SI 

** 

XlOl 

TO 

X105 

** 

M2 

** 

Xlll 

TO 

X115 

** 

M2 

+ 

M2 

** 

X121 

TO 

X125 

** 

Ml 

+ 

M2 

** 

X131 

TO 

X135 

** 

M2 

+ 

B1 

** 

X141 

TO 

X145 

** 

M2 

SI 

** 

XlSl 

TO 

X155 

** 

B2 

** 

X161 

TO 

X165 

** 

B2 

+ 

B2 

** 

X171 

TO 

X175 

** 

Ml 

B2 

** 

X181 

TO 

X185 

** 

M2 

+ 

B2 

** 

X191 

TO 

X195 

** 

B1 

+ 

B2 

** 

X201 

TO 

X2t95 

** 

B2 

+ 

SI 

** 

MAXIMIZE 

** 

DE 

TABLE  ** 

.8480 

Xll 

+ 

.7500 

X12 

+ 

.5500 

X13 

4- 

.8800 

X14 

4- 

.2500 

X15  4- 

.9600 

X21 

+ 

.8800 

X22 

+ 

.7000 

X23 

4- 

.  0000 

X24 

4- 

.4800 

X25  + 

.7300 

X31 

+ 

.6400 

X32 

+ 

.3200 

X33 

4- 

.7800 

X34 

4- 

.  1500 

X35  4- 

.9744 

X41 

+ 

.9375 

X42 

+ 

.7975 

X43 

4- 

.9856 

X44 

4- 

.4375 

X45  4- 

.9984 

X51 

+ 

.9856 

X52 

+ 

.9100 

X53 

4- 

.  0000 

X54 

4- 

.7296 

X55  4- 

.9271 

X61 

+ 

.8704 

X62 

+ 

.5376 

X63 

4- 

.9516 

X64 

4- 

.2775 

X65  4- 

.9936 

X71 

+ 

.9700 

X72 

+ 

.8650 

X73 

4- 

.  0000 

X74 

4- 

.6100 

X75  + 

.9568 

X81 

+ 

.9100 

X82 

+ 

.6940 

X83 

4- 

.9736 

X84 

4- 

.3625 

X85  4- 

.9892 

X91 

4- 

.9568 

X92 

+ 

.7960 

X93 

+ 

.  0000 

X94 

4- 

.5580 

X95  4- 

.9000 

X101 

+ 

.  8300 

X102 

+ 

.6200 

X103 

4- 

.9400 

X104 

4- 

.3600 

X105  + 

.9900 

Xlll 

+ 

.9711 

X112 

+ 

.8556 

X113 

4- 

.9964 

X114 

4- 

.5904 

X115  4- 

.9840 

X121 

+ 

.9575 

X122 

+ 

.8290 

X123 

4- 

.9928 

X124 

4- 

.5200 

X125  4- 

.9960 

X131 

+ 

.9796 

X132 

+ 

.8860 

X133 

4- 

.  0000 

X134 

4- 

.6672 

X135  + 

.9730 

X141 

+ 

.9388 

X142 

+ 

.7416 

X143 

4- 

.9868 

X144 

4- 

.4560 

X145  + 

.9800 

XlSl 

+ 

.9200 

X152 

+ 

.7400 

X153 

4- 

.  0000 

X154 

4- 

.5500 

X155  + 

.9996 

X161 

+ 

.9936 

X162 

+ 

.9324 

X163 

4- 

.0000 

X164 

+ 

.7975 

X165  4- 

.9968 

X171 

+ 

.9936 

X172 

4- 

.8830 

X173 

4- 

.  0000 

X174 

4- 

.6625 

X175  + 

.9980 

X181 

+ 

.9864 

X182 

4- 

.9012 

X183 

4- 

.0000 

X184 

4- 

.7120 

X185  + 

.9992 

X191 

+ 

.9904 

X192 

4- 

.9220 

X193 

4- 

.  0000 

X194 

4- 

.7660 

X195  4- 

.9946 

X201 

+ 

.9712 

X202 

4- 

.8232 

X203 

4- 

.  0000 

X204 

4- 

.6175 

X205 

CONSTRAINTS 

**  (1)  95  7.  DE  CRITERIA  ** 


.7500 

X12 

+  .8800 

X22 

+  .6400 

X32 

+  .9375 

X42 

+ 

.9856 

X52 

+  .8704 

X62 

+  .9700 

X72 

+  .9100 

X82 

+ 

.9568 

X92 

.  8300 

X102 

+  .9711 

X112 

+  .9575 

X122 

+ 

.9796 

X132 

+  .9388 

X142 

+  .9200 

X152 

+  .9936 

X162 

+ 

.9936 

X172 

+  .9864 

X182 

+  .9904 

X192 

+  .9712 

X202 

• 

r 

m 

• 

5700 

**  (2) 
.2500 

60  7 
X15 

;  OE  CRITERIA 
■«-  .4800  X25 

** 

-t-  .1500 

X35 

+  .4375 

X45 

+ 

.7296 

X55 

+  .2775 

X65 

.6100 

X75 

.3625 

X85 

+ 

.5580 

X95 

+  .3600 

X105 

+  .5904 

X115 

.5200 

X125 

+ 

.6672 

X135 

+  .4560 

X145 

+  .5500 

X155 

+  .7975 

X165 

+ 

.6625 

X175 

+  .7120 

X185 

+  .7660 

X195 

+  .6175 

X205 

iij 

0 

■ 

180 

**  (3) 
Xll  + 

TARGET  1  CONSTRAINT  ** 

X21  +  X31  +  X41  +  X51 

+  X61  +  X71 

+  ; 

X81 

+  X91 

XlOl 

Xlll 

+  X121 

+  X131  +  X141 

+  X151  +  X161  +  : 

X171 

+  X181 

X191  +  X201  .LE.  9000 

**  (4)  TARGET  2  CONSTRAINT  ** 

X12  +  X22  +  X32  +  X42  +  X52  +  X62  +  X72  +  X82  +  X92  + 

X102  +  X112  +  X122  +  X132  +  X142  +  X152  +  X162  +  X172  +  X182  + 

X192  +  X202  .LE.  6000 

**  (5)  TARGET  3  CONSTRAINT  ** 

X13  +  X23  +  X33  +  X43  +  X53  +  X63  +  X73  +  X83  +  X93  + 

X103  +  X113  +  X123  +  X133  +  X143  +  X153  +  X163  +  X173  X183  + 

X193  +  X203  .LE.  3000 

**  (6)  TARGET  4  CONSTRAINT  ** 

X14  +  X24  +  X34  +  X44  +  X54  +  X64  +  X74  +  X84  +  X94  + 

X104  +  X114  +  X124  +  X134  +  X144  +  X154  +  X164  +  X174  +  X184  + 

X194  +  X204  .LE.  1500 

**  (7)  TARGET  5  CONSTRAINTS  ** 

X15  +  X25  +  X35  +  X45  +  X55  +  X65  +  X75  +  X85  +  X95  + 

X105  +  X115  +  X125  +  X135  +  X145  +  X155  +  X165  +  X175  +  X185  + 

X195  +  X205  .LE.  300 

**  (8)  WEAPON  SYSTEM  1  (Ml)  CONSTRAINT  ♦♦ 

Xll  +  X12  +  X13  +  X14  +  X15  + 

2  X41  +  2  X42  +  2  X43  +  2  X44  +  2  X45  + 

X71  +  X72  +  X73  +  X74  +  X75  + 

X81  +  X82  +  X83  +  X84  +  X85  + 

X121  +  X122  +  X123  +  X124  +  X125  + 

X171  +  X172  +  X173  +  X174  +  X175  .LE.  2325 


**  (9)  WEAPON  SYSTEM  2  (M2)  CONSTRAINT  ** 
X101  +  X102  +  X103  +  X104  +  X105  + 

2  Xlll  +  2  X112  +  2  X113  +  2  X114  +  2  XH5  + 
X121  +  X122  +  X123  +  X124  +  X125  + 

X131  +  X132  +  X133  +  X134  +  X135  + 

X141  +  X142  +  X143  +  X144  +  X145  + 

X181  +  X182  +  X183  +  X184  +  X185  .  LE.  1250 

**  (10)  WEAPON  SYSTEM  3  (Bl)  CONSTRAINT 
X21  +  X22  +  X23  +  X24  +  X25  + 

2  X51  +  2  X52  +  2  X53  +  2  X54  +  2  X55  + 

X71  +  X72  +  X73  +  X74  +  X75  + 

X91  +  X92  +  X93  +  X94  +  X95  + 

X131  +  X132  +  X133  +  X134  +  X135  + 

X191  +  X192  +  X193  +  X194  +  X195  .LE.  2100 

**  (11)  WEAPON  SYSTEM  4  (B2)  CONSTRAINT  ** 
X151  +  X152  +  X153  +  X154  +  X155  + 

2  X161  +  2  X162  +  2  X163  +  2  X164  +  2  X165  + 
X171  +  X172  +  X173  +  X174  +  X175  + 

X181  +  X182  +  X183  +  X184  +  X185  + 

X191  +  X192  +  X193  +  X194  +  X195  + 

X201  +  X202  X203  X204  +  X205  .LE.  3000 

**  (12)  WEAPON  SYSTEM  5  (SI)  CONSTRAINT  *♦ 
X31  +  X32  +  X33  +  X34  +  X35  + 

2  X61  +  2  X62  +  2  X63  +  2  X64  +  2  X65  + 
xai  +  xa2  +  xa3  +  X84  +  xas  + 

X91  +  X92  +  X93  +  X94  +  X95  + 

X141  +  X142  +  X143  +  X144  +  X145  + 

X201  +  X202  +  X203  +  X204  +  X205  .LE.  9000 

PRINT 

OPTIMIZE 


************  MPOS  INPUT  FILE  ************** 
************  6  FACTORS  ************** 

******************************************* 

REGULAR 

TITLE 

RESPONSE  SURFACE  FOR  ARSENAL  EXCHANGE  MODEL 


VARIABLES  **  WEAPON  ** 

**  SYSTEM  ** 


Xll 

TO 

X15 

** 

Ml 

*« 

X2l 

TO 

X25 

** 

B1 

** 

X31 

TO 

X35 

** 

SI 

** 

X41 

TO 

X45 

** 

Ml 

+ 

Ml 

** 

X51 

TO 

X55 

** 

B1 

+ 

B1 

** 

X61 

TO 

X65 

** 

SI 

+ 

SI 

** 

X71 

TO 

X75 

** 

Ml 

+ 

B1 

** 

X81 

TO 

X8S 

** 

Ml 

+ 

SI 

** 

X91 

TO 

X95 

** 

B1 

+ 

SI 

** 

XlOl 

TO 

X10S 

** 

M2 

** 

Xlll 

TO 

X115 

** 

M2 

+ 

M2 

** 

X121 

TO 

X125 

** 

Ml 

+ 

M2 

** 

X131 

TO 

X135 

** 

M2 

+ 

B1 

*« 

X141 

TO 

X145 

** 

M2 

+ 

SI 

*« 

X151 

TO 

X155 

** 

B2 

** 

X161 

TO 

X165 

** 

B2 

+ 

B2 

** 

X171 

TO 

X175 

** 

Ml 

+ 

B2 

** 

X181 

TO 

X185 

** 

M2 

+ 

B2 

*♦ 

X191 

TO 

X195 

** 

B1 

+ 

B2 

** 

X201 

TO 

X205 

** 

B2 

+ 

SI 

** 

X211 

TO 

X215 

** 

M3 

** 

X221 

TO 

X225 

** 

M3 

+ 

M3 

** 

X231 

TO 

X235 

** 

Ml 

+ 

M3 

** 

X241 

TO 

X245 

** 

M2 

+ 

M3 

** 

X251 

TO 

X255 

** 

B1 

+ 

M3 

** 

X261 

TO 

X265 

** 

B2 

+ 

M3 

** 

X271 

TO 

X275 

** 

M3 

+ 

SI 

** 

A  -  17 


MAXIMIZE 


**  DE  TABLE  ** 


8400 

Xll 

+ 

.7500 

X12 

+ 

.5500 

X13 

+ 

.8800 

X14 

+ 

.2500 

X15  + 

9600 

X21 

+ 

.8800 

X22 

+ 

.7000 

X23 

+ 

.  0000 

X24 

+ 

.4800 

X25  + 

7300 

X31 

+ 

.6400 

X32 

+ 

.3200 

X33 

+ 

.7800 

X34 

+ 

.  1500 

X35  + 

9744 

X41 

+ 

.9375 

X42 

+ 

.7975 

X43 

+ 

.9856 

X44 

+ 

.4375 

X45  + 

9984 

X51 

+ 

.9856 

X52 

+ 

.9100 

X53 

+ 

.0000 

X54 

+ 

.7296 

X55  + 

9271 

X61 

+ 

.8704 

X62 

+ 

.5376 

X63 

+ 

.9516 

X64 

+ 

.2775 

X65  + 

9936 

X71 

+ 

.9700 

X72 

+ 

.8650 

X73 

+ 

.  0000 

X74 

+ 

.6100 

X75  + 

9568 

X81 

+ 

.9100 

X82 

+ 

.6940 

X83 

+ 

.9736 

X84 

+ 

.3625 

X85  + 

9892 

X91 

+ 

.9568 

X92 

+ 

.7960 

X93 

+ 

.  0000 

X94 

+ 

.  5580 

X95  + 

9000 

X101 

+ 

.8300 

X102 

+ 

.6200 

X103 

+ 

.9400 

X104 

+ 

.3600 

X105  + 

9900 

Xlll 

+ 

.9711 

X112 

+ 

.8556 

X113 

+ 

.9964 

X114 

+ 

.5904 

X115  + 

9840 

X121 

+ 

.9575 

X122 

-*• 

.8290 

X123 

+ 

.9928 

X124 

+ 

.5200 

X125  + 

9960 

X131 

+ 

.9796 

X132 

+ 

.8860 

X133 

+ 

.0000 

X134 

.6672 

X135  + 

9730 

X141 

+ 

.9388 

X142 

+ 

.7416 

X143 

.9868 

X144 

+ 

.4560 

X145  + 

9800 

X151 

.9200 

X152 

+ 

.7400 

X153 

+ 

.  0000 

X154 

+ 

.5500 

X155  + 

9996 

X161 

+ 

.9936 

X162 

+ 

.9324 

X163 

+ 

.0000 

X164 

+ 

.7975 

X165  + 

9968 

X171 

+ 

.9936 

X172 

+ 

.  8830 

X173 

+ 

.  0000 

X174 

+ 

.6625 

X175  + 

9980 

X181 

+ 

.9864 

X182 

+ 

.9012 

X183 

+ 

.  0000 

X184 

+ 

.7120 

X185  + 

9992 

X191 

+ 

.9904 

X192 

.9220 

X193 

+ 

.  0000 

X194 

+ 

.7660 

X195  + 

9946 

X201 

+ 

.9712 

X202 

+ 

.8232 

X203 

+ 

.0000 

X204 

+ 

.6175 

X205  + 

9200 

X211 

+ 

.8700 

X212 

+ 

.6500 

X213 

+ 

.9500 

X214 

+ 

.4500 

X215  + 

9936 

X221 

+ 

.9831 

X222 

+ 

.8775 

X223 

.9975 

X224 

+ 

.6975 

X225 

9872 

X231 

+ 

.9675 

X232 

+ 

.8425 

X233 

.9940 

X234 

+ 

.5875 

X235 

9920 

X241 

+ 

.9779 

X242 

+ 

.8670 

X243 

+ 

.9970 

X244 

+• 

.6480 

X245  + 

9968 

X251 

+ 

.9844 

X252 

+ 

.8950 

X253 

+ 

.0000 

X254 

+ 

.7140 

X255  + 

9984 

X261 

+ 

.9896 

X262 

+ 

.9090 

X263 

+ 

.  0000 

X264 

.7525 

X265  + 

9784 

X271 

+ 

.9532 

X272 

+ 

.7620 

X273 

+ 

.9890 

X274 

+ 

.5325 

X275 

CONSTRAINTS 


**  (1) 

95  *> 

i  DE  CRITERIA 

** 

.7500 

X12 

+ 

.8800 

X22 

-1-  .6400 

X32 

+ 

.9375 

X42 

+ 

.9856 

X52 

+ 

.8704 

X62 

+ 

.9700 

X72 

+  .9100 

X82 

+ 

.9568 

X92 

+ 

.  8300 

X102 

+ 

.9711 

X112 

+ 

.9575 

X122 

+  .9796 

X132 

+ 

.9388 

X142 

+ 

.9200 

X152 

+ 

.9936 

X162 

+ 

.9936 

X172 

+  .9864 

X182 

+ 

.9904 

X192 

+ 

.9712 

X202 

+ 

.8700 

X212 

+ 

.9831 

X222 

+  .9675 

X232 

+ 

.9779 

X242 

+ 

.9844 

X252 

+ 

.9896 

X262 

+ 

.9532 

X272 

.LE.  5700 

♦♦  (2) 

60 

i  DE  CRITERIA 

«* 

.2500 

X15 

+ 

.4800 

X25 

+  . 1500 

X35 

+ 

.4375 

X45 

+ 

.7296 

X55 

+ 

.2775 

X65 

.6100 

X75 

-t-  .  3625 

X85 

+ 

.5580 

X95 

+ 

.3600 

X105 

+ 

.5904 

X115 

+ 

.5200 

X125 

+  .6672 

X135 

+ 

.4560 

X145 

+ 

.5500 

X155 

+ 

.7975 

X165 

+ 

.  6625 

X175 

+  .7120 

X185 

+ 

.7660 

X195 

+ 

.6175 

X205 

+ 

.  4500 

X215 

+ 

.6975 

X225 

+  . 5875 

X235 

+ 

.6480 

X245 

+ 

.7140 

X255 

+ 

.7525 

X265 

+ 

.5325 

X275 

.GE.  180 

A  -  18 


**  (3)  TARGET  1  CONSTRAINT  ** 

Xll  +  X21  +  X31  +  X41  +  X51  +  X61  +  X71  +  X81  + 

X91  +  X101  +  Xlll  +  X121  +  X131  +  X141  +  X151  +  X161  + 

X171  +  X181  +  X191  +  X201  +  X211  +  X221  +  X231  +  X241  + 

X251  +  X261  +  X271  .LE.  9000 

**  (4)  TARGET  2  CONSTRAINT  ** 

X12  +  X22  +  X32  +  X42  +  X52  +  X62  +  X72  +  X82  + 

X92  +  X102  +  X112  +  X122  +  X132  +  X142  +  X152  +  X162  + 

X172  +  X182  +  X192  +  X202  +  X212  +  X222  +  X232  +  X242  + 

X252  +  X262  +  X272  .LE.  6000 

**  (5)  TARGET  3  CONSTRAINT  ** 

X13  +  X23  +  X33  +  X43  +  X53  +  X63  +  X73  +  X83  + 

X93  +  X103  +  X113  +  X123  +  X133  +  X143  +  X153  +  X163  + 

X173  +  X183  +  X193  +  X203  +  X213  +  X223  +  X233  +  X243  + 

X253  +  X263  +  X273  .LE.  3000 

**  (6)  TARGET  4  CONSTRAINT  ** 

X14  +  X24  +  X34  +  X44  +  X54  +  X64  +  X74  +  X84  + 

X94  +  X104  +  X114  +  X124  +  X134  +  X144  +  X154  +  X164  + 

X174  +  X184  +  X194  +  X204  +  X214  +  X224  +  X234  +  X244  + 

X254  +  X264  +  X274  .LE.  1500 

**  (7)  TARGET  5  CONSTRAINTS  ** 

X15  +  X25  +  X35  +  X45  +  X55  +  X65  +  X75  +  X85  + 

X95  +  X105  +  X115  +  X125  +  X135  +  X145  +  X155  +  X165  + 

X175  +  X185  +  X195  +  X205  +  X215  +  X225  +  X235  +  X245  + 

X255  +  X265  +  X275  .LE.  300 

**  (8)  WEAPON  SYSTEM  1  (Ml)  CONSTRAINT  ** 

Xll  +  X12  +  X13  +  X14  +  X15  + 

2  X41  +  2  X42  +  2  X43  +  2  X44  +  2  X45  + 

X71  +  X72  +  X73  +  X74  +  X75  + 

X81  +  X82  +  X83  +  X84  +  X85  + 

X121  +  X122  +  X123  +  X124  +  X125  + 

X171  +  X172  +  X173  +  X174  +  X175  + 

X231  +  X232  +  X233  +  X234  +  X235  .LE.  2325 

**  (9)  WEAPON  SYSTEM  2  <M2)  CONSTRAINT  ** 

X101  +  X102  +  X103  +  X104  +  X105  + 

2  Xlll  +  2  X112  +  2  X113  +  2  X114  +  2  X115  + 

X121  +  X122  +  X123  +  X124  +  X125  + 

X131  +  X132  +  X133  +  X134  +  X135  + 

X141  +  X142  +  X143  +  X144  +  X145  + 

X181  +  X182  +  X183  +  X184  +  X185  + 

X241  +  X242  +  X243  +  X244  +  X245  .LE.  1250 


**  (10)  WEAPON  SYSTEM  3  (M3)  CONSTRAINT  ** 
X211  +  X212  +  X213  +  X214  +  X215  + 

2  X221  +  2  X222  +  2  X223  +  2  X224  +  2  X225  + 
X231  +  X232  +  X233  +  X234  +  X235  + 

X241  +  X242  +  X243  +  X244  +  X245  + 

X251  +  X252  +  X253  +  X254  +  X255  + 

X261  +  X262  +  X263  +  X264  +  X265  + 

X271  +  X272  +  X273  +  X274  +  X275  .LE.  700 

**  (11)  WEAPON  SYSTEM  4  (Bl)  CONSTRAINT 
X21  +  X22  +  X23  +  X24  +  X25  + 

2  X51  +  2  X52  +  2  X53  +  2  X54  +  2  X55  + 

X71  +  X72  +  X73  +  X74  +  X75  + 

X91  +  X92  +  X93  +  X94  +  X95  + 

X131  +  X132  +  X133  +  X134  +  X135  + 

X191  +  X192  +  X193  +  X194  +  X195  + 

X251  +  X252  +  X253  +  X254  +  X255  .LE.  2100 

*■*  (12)  WEAPON  SYSTEM  5  (B2)  CONSTRAINT  ** 
X151  +  X152  +  X153  +  X154  +  X155  + 

2  X161  +  2  X162  +  2  X163  +  2  X164  +  2  X165  + 
X171  +  X172  +  X173  +  X174  +  X175  + 

xiai  +  xia2  +  xia3  +  xia4  +  xias  + 

X191  +  X192  +  X193  +  X194  +  X195  + 

X201  +  X202  +  X203  +  X204  +  X205  + 

X261  +  X262  +  X263  +  X264  X265  .LE.  3000 

**  (13)  WEAPON  SYSTEM  6  (SI)  CONSTRAINT  ** 
X31  +  X32  +  X33  +  X34  +  X35  + 

2  X61  +  2  X62  +  2  X63  +  2  X64  +  2  X65  + 
xai  +  xa2  +  X83  +  xa4  +  xas  + 

X91  +  X92  +  X93  +  X94  +  X95  + 

X141  +  X142  +  X143  +  X144  +  X145  + 

X201  +  X202  +  X203  +  X204  +  X205  + 

X271  +  X272  +  X273  +  X274  +  X275  .LE.  9000 

PRINT 

OPTIMIZE 


1l  Experimental  Desi.ans 


TABLE  B. 1 

BOX-BEHNKEN  EXPERIMENTAL  DESIGN  FOR  THREE  FACTORS 
NON-CODED  CODED 


Ml 

B1 

SI 

Ml 

B1 

SI 

RESPONSE 

3000 

3000 

9000 

1 

1 

0 

10131.4892 

3000 

1200 

9000 

1 

-1 

0 

9242.6923 

1650 

3000 

9000 

-1 

1 

0 

9671.9492 

1650 

1200 

9000 

-1 

-1 

0 

8493. 1923 

3000 

2100 

12000 

1 

0 

1 

10395.4062 

3000 

2100 

6000 

1 

0 

-1 

8477.3846 

1650 

2100 

12000 

-1 

0 

1 

9935.8662 

1650 

2100 

6000 

-1 

0 

-1 

7464.8846 

2325 

3000 

12000 

0 

1 

1 

10584.9202 

2325 

3000 

6000 

0 

1 

-1 

8757.9423 

2325 

1200 

12000 

0 

-1 

1 

9724.2762 

2325 

1200 

6000 

0 

-1 

-1 

7179. 1346 

2325 

2100 

9000 

0 

0 

0 

9474.4362 

TABLE  B.2 


BOX-DRAPER  EXPERIMENTAL  DESISN  FOR  THREE  FACTORS 
NON-CQDED  CODED 


RUN* 

Ml 

B1 

SI 

Ml 

B1 

SI 

RESPONSE 

1 

2914 

2886 

11619 

0.873 

0.873 

0.873 

10642.0176 

2 

2914 

2886 

6381 

0.873 

0.873 

-0.873 

9189. 1323 

3 

2914 

1314 

11619 

0.873 

-0.873 

0.873 

9892.8949 

4 

2914 

1314 

6381 

0.873 

-0.873 

-0.873 

7965.0446 

5 

1736 

2886 

11619 

-0. 873 

0.873 

0.873 

10251.0156 

b 

1736 

2886 

6381 

-0.873 

0.873 

-0.873 

8464.9046 

7 

1736 

1314 

11619 

-0. 873 

-0.873 

0.873 

9491.9037 

8 

1736 

1314 

6381 

-0. 873 

-0.873 

-0.873 

7081.5446 

9 

1334 

2100 

9000 

-1.468 

0 

0 

8943.0723 

10 

3316 

2100 

9000 

1.468 

0 

0 

9811.7725 

11 

2325 

779 

9000 

0 

-1.468 

0 

8576.7423 

12 

2325 

3421 

9000 

0 

1.468 

0 

10099.7576 

13 

2325 

2100 

4596 

0 

0 

-1.468 

7071.7846 

14 

2325 

2100 

13404 

0 

0 

1.468 

10471.3739 

15 

2325 

2100 

9000 

0 

0 

0 

9474.4362 

B  -  2 
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TABLE  B.3 


MINIMUM  POINT  EXPERIMENTAL  DESIGN  FOR  THREE  FACTORS 
NON-CODED  CODED 


RUN# 

Ml 

B1 

SI 

Ml 

B1 

SI 

RESPONSE 

1 

1650 

1200 

6000 

-1 

-1 

-1 

6672.8846 

2 

3000 

1200 

6000 

1 

-1 

-1 

7685.3846 

3 

1650 

3000 

6000 

-1 

1 

-1 

8256.8846 

4 

1650 

1200 

12000 

-1 

-1 

1 

9494.5061 

5 

2457 

2280 

6000 

0. 1925 

0. 1925 

-1 

8228.5346 

h 

2457 

1200 

9600 

0. 1925 

-1 

0. 1925 

9208.0423 

7 

1650 

2280 

9600 

-1 

0. 1925 

0. 1925 

9471.1782 

8 

2127 

3000 

12000 

-0.2912 

1 

1 

10520. 1731 

9 

3000 

1836 

12000 

1 

-0.2912 

1 

10265.9406 

10 

3000 

3000 

8100 

1 

1 

-0.2912 

9924. 1292 

B  - 
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NATiONftl  BUREAU  OF  STANDARDS* 


TABLE  B.4 


KOSHAL  EXPERIMENTAL  DESIGN  FOR  THREE  FACTORS 
NQNzQQDED  coded 


RUN# 

Ml 

B1 

SI 

Ml 

B1 

SI 

RESPONSE 

1 

1988 

1650 

7500 

~0>  5 

-0.5 

-0.5 

8282.3846 

2 

2662 

1650 

7500 

0.5 

-0.5 

—0.  5 

8695.0323 

3 

1988 

2550 

7500 

-0.5 

0.5 

-0.5 

8940. 8523 

4 

1988 

1650 

10500 

-0.5 

-0.5 

0.5 

9484.6413 

S 

3337 

1650 

7500 

1.5 

-0.5 

-0.5 

9079.7823 

6 

1988 

3450 

7500 

—0.  5 

1.5 

-0.5 

9570.8523 

7 

1988 

1650 

13500 

-0.5 

-0.5 

in 

• 

10175.8413 

8 

2662 

2550 

7500 

0.5 

0.5 

-0.5 

9325.0323 

9 

2662 

1650 

10500 

0.5 

-0.5 

in 

• 

8 

9714.0709 

10 

1988 

2550 

10500 

-0.5 

0.5 

0.5 

9920.9244 

TABLE  B.5 


CENTRAL  COMPOSITE  EXPERIMENTAL  DESIGN  FOR  THREE  FACTORS 


NON-CODED 

CODED 

RUN* 

Ml 

B1 

SI 

Ml 

B1 

SI 

RESPONSE 

1 

3000 

3000 

12000 

1 

1 

1 

10792.2127 

2 

3000 

3000 

6000 

1 

1 

-1 

9142.6923 

3 

3000 

1200 

12000 

1 

-1 

1 

9954.0461 

4 

3000 

1200 

6000 

1 

-1 

-1 

7685.3846 

5 

1650 

3000 

12000 

-1 

1 

1 

10362.5723 

6 

1650 

3000 

6000 

-1 

1 

-1 

8256.8846 

7 

1650 

1200 

12000 

-1 

-1 

1 

9494.5061 

8 

1650 

1200 

6000 

-1 

-1 

-1 

6672.8846 

9 

3460 

2100 

9000 

1.682 

0 

0 

9860.7901 

10 

1190 

2100 

9000 

-1.682 

0 

0 

8860.9923 

11 

2325 

3614 

9000 

0 

1.682 

.  0 

9978.8648 

12 

2325 

586 

9000 

0 

-1.682 

0 

8441.6423 

13 

2325 

2100 

14046 

0 

0 

1.682 

10597.9122 

14 

2325 

2100 

3954 

0 

0 

-1.682 

6654.4846 

15 

2325 

2100 

9000 

0 

0 

0 

9474.4362 

TABLE  B.6 


HYBRID  (31 lA)  EXPERIMENTAL  DESIGN  FOR  THREE  FACTORS 
NON-CODED  CODED 


RUN# 

Ml 

B1 

SI 

Ml 

B1 

SI 

RESPONSE 

1 

2325 

2100 

15000 

0 

0 

2 

10785.9456 

2 

2325 

2100 

3000 

0 

0 

-2 

6034.3846 

3 

1370 

827 

12000 

-1.414 

-1.414 

1 

9204. 3765 

4 

3280 

827 

12000 

1.414 

-1.414 

1 

9866.4389 

5 

1370 

3373 

12000 

-1.414 

1.414 

1 

10438.0595 

6 

3280 

3373 

12000 

1.414 

1.414 

1 

11037.5089 

7 

3675 

2100 

6000 

2 

0 

-1 

8897.4423 

e 

975 

2100 

6000 

-2 

0 

-1 

6958.3846 

9 

2325 

3900 

6000 

0 

2 

-1 

9387.9423 

10 

2325 

300 

6000 

0 

-2 

-1 

6387. 1346 

11 

2325 

2100 

9000 

0 

0 

0 

9474.4362 

TABLE  B.7 


HYBRID  (310+CP)  EXPERIMENTAL  DESIGN  FOR  THREE  FACTORS 
NON-CODED  CODED 


RUN# 

Ml 

B1 

SI 

Ml 

B1 

SI 

RESPONSE 

1 

2325 

2100 

12872 

0 

0 

1 . 2906 

10366.5168 

2 

2325 

2100 

8592 

0 

0 

-0. 136 

9320.2623 

3 

1650 

1200 

10916 

-1 

-1 

0.6386 

9244.7525 

4 

3000 

1200 

10916 

1 

-1 

0. 6386 

9704.2925 

5 

1650 

3000 

10916 

-1 

1 

0.6386 

10113.3956 

6 

3000 

3000 

10916 

1 

1 

0. 6386 

10569.6987 

7 

3117 

2100 

6218 

1.1736 

0 

-0.9273 

8679.6623 

8 

1533 

2100 

6218 

-1.1736 

0 

-0.9273 

7516.6546 

9 

2325 

3156 

6218 

0 

1.1736 

-0.9273 

8967.4223 

10 

2325 

1044 

6218 

0 

-1. 1736 

-0.9273 

7181.3746 

11 

2325 

2100 

9000 

0 

0 

0 

9474.4361 

TABLE  B.8 


FULL  FACTORIAL  EXPERIMENTAL  DESIGN  FOR  THREE  FACTORS 
NON-CODED  CODED 


RUN* 

Ml 

B1 

SI 

Ml 

B1 

SI 

RESPONSE 

1 

3000 

3000 

12000 

1 

1 

1 

10792.2129 

2 

3000 

3000 

9000 

1 

1 

0 

10131.4893 

3 

3000 

3000 

6000 

1 

1 

-1 

9142.6924 

4 

3000 

2100 

12000 

1 

0 

1 

10395.4063 

5 

3000 

2100 

9000 

1 

0 

0 

9704.2061 

6 

3000 

2100 

6000 

1 

0 

-1 

8477.3848 

7 

3000 

1200 

12000 

1 

-1 

1 

9954.0459 

8 

3000 

1200 

9000 

1 

-1 

0 

9242.6924 

9 

3000 

1200 

6000 

1 

-1 

-1 

7685.3848 

10 

2325 

3000 

12000 

0 

1 

1 

10534.9199 

11 

2325 

3000 

9000 

0 

1 

0 

9901.7188 

12 

2325 

3000 

6000 

0 

1 

-1 

8757.9424 

13 

2325 

2100 

12000 

0 

0 

1 

10165.6357 

14 

2325 

2100 

9000 

0 

0 

0 

9474.4365 

15 

2325 

2100 

6000 

0 

0 

-1 

7971.1348 

16 

2325 

1200 

12000 

0 

-1 

1 

9724.2764 

17 

2325 

1200 

9000 

0 

-1 

0 

8871.4424 

18 

2325 

1200 

6000 

0 

-1 

-1 

7179. 1348 

B  -  8 


TABLE  B.8  (continued) 


FULL  FACTORIAL  EXPERIMENTAL 
NQN-CODED 


RUN# 

Ml 

B1 

SI 

Ml 

19 

1650 

3000 

12000 

-1 

20 

1650 

3000 

9000 

-1 

21 

1650 

3000 

6000 

-1 

22 

1650 

2100 

12000 

-1 

23 

1650 

2100 

9000 

-1 

24 

1650 

2100 

6000 

-1 

25 

1650 

1200 

12000 

-1 

26 

1650 

1200 

9000 

-1 

27 

1650 

1200 

6000 

-1 

DESIGN  FOR 

THREE 

FACTORS 

CODED 

B1 

SI 

RESPONSE 

1 

1 

10362.5723 

1 

0 

9671.9492 

1 

-1 

8256.8848 

0 

1 

9935.8662 

0 

0 

9123. 1924 

0 

-1 

7464.8848 

-1 

1 

9494.5059 

-1 

0 

8493. 1924 

-1 

-1 

6672.8848 

B  -  9 


TABLE 

B.9 

BOX-BEHNKEN  EXPERIMENTAL 

DESIGN 

FOR  FOUR  FACTORS 

NON-CODED 

CODED 

RUN#  Ml 

M2 

B1 

SI 

Ml 

M2 

B1 

SI 

RESPONSE 

1 

3000 

2000 

2100 

9000 

1 

1 

0 

0 

10530.8908 

2 

3000 

500 

2100 

9000 

1 

-1 

0 

0 

9914.4062 

3 

1650 

2000 

2100 

9000 

-1 

1 

0 

0 

10085.4662 

4 

1650 

500 

2100 

9000 

-1 

-1 

0 

0 

9445. 1923 

5 

2325 

1250 

3000 

12000 

0 

0 

1 

1 

11043.7952 

6 

2325 

1250 

3000 

6000 

0 

0 

1 

-1 

9570.4423 

7 

2325 

1250 

1200 

12000 

0 

0 

-1 

1 

10242.4877 

8 

2325 

1250 

1200 

6000 

0 

0 

-1 

-1 

8216.6346 

9 

3000 

1250 

2100 

12000 

1 

0 

0 

1 

10862.7746 

0 

3000 

1250 

2100 

6000 

1 

0 

0 

-1 

9325. 1923 

1 

1650 

1250 

2100 

12000 

-1 

0 

0 

1 

10444.2508 

2 

1650 

1250 

2100 

6000 

-1 

0 

0 

-1 

8502.3846 

3 

2325 

2000 

3000 

9000 

0 

1 

1 

0 

10714.9038 

4 

2325 

2000 

1200 

9000 

0 

1 

-1 

0 

9873.8762 

5 

2325 

500 

3000 

9000 

0 

-1 

1 

0 

10111.9192 

6 

2325 

500 

1200 

9000 

0 

-1 

-1 

0 

9186.4423 

7 

3000 

1250 

3000 

9000 

1 

0 

1 

0 

10645.8738 

8 

3000 

1250 

1200 

9000 

1 

0 

-1 

0 

9788.3462 

9 

1650 

1250 

3000 

9000 

-1 

0 

1 

0 

10197.4492 

0 

1650 

1250 

1200 

9000 

-1 

0 

-1 

0 

9287.6923 

1 

2325 

2000 

2100 

12000 

0 

1 

0 

1 

10930.8071 

2 

2325 

2000 

2100 

6000 

0 

1 

0 

-1 

9427.9423 

3 

2325 

500 

2100 

12000 

0 

-1 

0 

1 

10374.5477 

4 

2325 

500 

2100 

6000 

0 

-1 

0 

-1 

8386. 1346 

5 

2325 

1250 

2100 

9000 

0 

0 

0 

0 

9999.9361 

B 

-  10 

j*  ^  y*  ^ 

TABLE  B. 10 


I 

I 


> 


V 


BOX-DRAPER  EXPERIMENTAL  DESIGN  FOR  FOUR  FACTORS 


NON-CODED  CODED 


RUN# 

Ml 

M2 

B1 

SI 

Ml 

M2 

B1 

SI 

RESPONSE 

1 

2781 

1756 

2707 

11025 

.675 

.675 

.675 

.675 

11052.2726 

2 

2781 

1756 

2707 

6975 

.675 

.675 

.675  - 

.675 

10190.7768 

3 

2781 

1756 

1493 

11025 

.675 

.675 

-.675 

.675 

10514.5027 

4 

2781 

1756 

1493 

6975 

.675 

.675 

-.675  - 

.675 

9542.6223 

5 

2781 

744 

2707 

11025 

.675  - 

.675 

.675 

.675 

10680.2466 

6 

2781 

744 

2707 

6975 

.675  - 

.675 

.675  - 

.675 

9744.8623 

7 

2781 

744 

1493 

11025 

.675  - 

.675 

-.675 

.675 

10111.3233 

8 

2781 

744 

1493 

6975 

.675  - 

.675 

-.675  - 

.675 

8895.0623 

9 

1869 

1756 

2707 

11025 

-.675 

.675 

.675 

.675 

10772. 1974 

10 

1869 

1756 

2707 

6975 

-.675 

.675 

.675  - 

.675 

9880.3320 

11 

1869 

1756 

1493 

11025 

-.675 

.675 

-.675 

.675 

10220.2949 

12 

1869 

1756 

1493 

6975 

-.675 

.675 

-.675  - 

.675 

9033.0223 

13 

1869 

744 

2707 

11025 

-.675  - 

.675 

.675 

.675 

10385.9703 

14 

1869 

744 

2707 

6975 

-.675  - 

.675 

.675  - 

.675 

9225.0223 

15 

1869 

744 

1493 

11025 

-.675  - 

.675 

-.675 

.675 

9800. 8786 

16 

1869 

744 

1493 

6975 

-.675  - 

.675 

-.675  - 

.675 

8336.4946 

17 

1414 

1250 

2100 

9000 

-1.349 

0 

0 

0 

9689.8318 

18 

3236 

1250 

2100 

9000 

1.349 

0 

0 

0 

10310.0405 

19 

2325 

238 

2100 

9000 

0 

-1.349  0 

0 

9574.4913 

20 

2325 

2262 

2100 

9000 

0 

1 . 349  0 

0 

10414.2856 

21 

2325 

1250 

886 

9000 

0 

0 

-1.349 

0 

9404.5905 

22 

2325 

1250 

3314 

9000 

0 

0 

1.349 

0 

10566.3179 

23 

2325 

1250 

2100 

4953 

0 

0 

0-1 

.349 

8338.5546 

24 

2325 

1250 

2100 

13047 

0 

0 

0  1 

.349 

10861.8458 

25 

2325 

1250 

2100 

9000 

0 

0 

0 

0 

9999.9361 

TABLE  B. 11 


MINIMUM  POINT  EXPERIMENTAL  DESIGN  FOR  FOUR  FACTORS 


NDN-CODED 

CODED 

UN« 

Ml 

M2 

B1 

SI 

Ml 

M2 

B1 

SI 

RESPONSE 

1 

1650 

500 

1200 

6000 

-1 

-1 

-1 

-1 

7087.8846 

2 

3000 

500 

1200 

6000 

1 

-1 

-1 

-1 

8100.3846 

3 

1650 

2000 

1200 

6000 

-1 

1 

-1 

-1 

8332.8846 

4 

1650 

500 

3000 

6000 

-1 

-1 

1 

-1 

8671.8846 

5 

1650 

500 

1200 

12000 

-1 

-1 

-1 

1 

9704.7061 

6 

2604 

1560 

1200 

6000 

.4114  . 

4114 

-1 

-1 

8690.9723 

7 

2604 

500 

2472 

6000 

.4114 

-1 

.4114 

-1 

8872.3723 

8 

2604 

500 

1200 

10500 

.4114 

-1 

-1 

.4114 

9683.8477 

9 

1650 

1560 

2472 

6000 

-1  . 

4114 

.4114 

-1 

9017.5923 

10 

1650 

1560 

1200 

10500 

-1  . 

4114 

-1 

.4114 

9804.7301 

11 

1650 

500 

2472 

10500 

-1 

-1 

.4114 

.4114 

9979.3780 

12 

1886 

2000 

3000 

12000 

-.6502 

1 

1 

1 

11181.9322 

13 

3000 

760 

3000 

12000 

1  -. 

6502 

1 

1 

11071.2087 

14 

3000 

2000 

1512 

12000 

1 

1  - 

.6502 

1 

10875.2048 

15 

3000 

2000 

3000 

6900 

1 

1 

1  -. 

6502 

10485.3723 

TABLE  B. 12 


KOSHAL  EXPERIMENTAL  DESIGN  FOR  FOUR  FACTORS 


NON-CODED 

CODED 

Ml 

M2 

B1 

SI 

Ml 

M2 

B1 

SI 

RESPONSE 

2055 

950 

1740 

7800 

-.4 

-.4 

-.4 

-.4 

9167.5423 

2730 

950 

1740 

7800 

.6 

-.4 

-.4 

-.4 

9542.6923 

2055 

1700 

1740 

7800 

-.4 

.6 

-.4 

-.4 

9643.9423 

2055 

950 

2640 

7800 

-.4 

-.4 

.6 

-.4 

9763.3672 

2055 

950 

1740 

10800 

-.4 

-.4 

-.4 

.6 

10020.0841 

3405 

950 

1740 

7800 

1.6 

-.4 

-.4 

-.4 

9788.4241 

2055 

2450 

1740 

7800 

-.4 

1.6 

-.4 

-.4 

9959.4841 

2055 

950 

3540 

7800 

-.4 

-.4 

1.6 

-.4 

10186.7272 

2055 

950 

1740 

13800 

-.4 

-.4 

-.4 

1.6 

10656.2591 

2730 

1700 

1740 

7800 

•  6 

.6 

-.4 

-.4 

9873.9541 

2730 

950 

2640 

7800 

•  6 

-.4 

•  6 

-.4 

9993. 1372 

2730 

950 

1740 

10800 

.6 

-.4 

-.4 

.6 

10249.8541 

2055 

1700 

2640 

7800 

-.4 

.6 

.6 

-.4 

10078.6672 

2055 

1700 

1740 

10800 

-.4 

.6 

-.4 

.6 

10326.9688 

2055 

950 

2640 

10800 

-.4 

-.4 

.6 

.6 

10449. 1518 

TABLE  B. 13 


CENTRAL  COMPOSITE  EXPERIMENTAL  DESIGN  FOR  FOUR  FACTORS 


NON-CODED 

CODED 

RUN#  Ml 

M2 

B1 

SI 

Ml 

M2 

B1 

SI 

RESPONSE 

1 

3000 

2000 

3000 

12000 

1 

1 

1 

1 

11453.3481 

2 

3000 

2000 

3000 

6000 

1 

1 

1 

-1 

10281.0892 

3 

3000 

2000 

1200 

12000 

1 

1 

-1 

1 

10735.7096 

4 

3000 

2000 

1200 

6000 

1 

1 

-1 

-1 

9182.6923 

5 

3000 

500 

3000 

12000 

1 

-1 

1 

1 

10975.7627 

6 

3000 

500 

3000 

6000 

1 

-1 

1 

-1 

9467.6923 

7 

3000 

500 

1200 

12000 

1 

-1 

-1 

1 

10164.2461 

8 

3000 

500 

1200 

6000 

1 

-1 

-1 

-1 

8100.3846 

9 

1650 

2000 

3000 

12000 

-1 

1 

1 

1 

11111. 8277 

10 

1650 

2000 

3000 

6000 

-1 

1 

1 

-1 

9673. 1923 

11 

1650 

2000 

1200 

12000 

-1 

1 

-1 

1 

10315. 1908 

12 

1650 

2000 

1200 

6000 

-1 

1 

-1 

-1 

8332. 8846 

13 

1650 

500 

3000 

12000 

-1 

-1 

1 

1 

10559.2338 

14 

1650 

500 

3000 

6000 

-1 

-1 

1 

-1 

8671.8846 

15 

1650 

500 

1200 

12000 

-1 

-1 

-1 

1 

9704.7061 

16 

1650 

500 

1200 

6000 

-1 

-1 

-1 

-1 

7087.8846 

17 

975 

1250 

2100 

9000 

-2 

0 

0 

0 

9540.3961 

18 

3675 

1250 

2100 

9000 

2 

0 

0 

0 

10459.4762 

19 

2325 

0 

2100 

9000 

0 

-2 

0 

0 

9474.4362 

20 

2325 

2750 

2100 

9000 

0 

2 

0 

0 

10599.9208 

21 

2325 

1250 

300 

9000 

0 

0 

-2 

0 

9028.9423 

22 

2325 

1250 

3900 

9000 

0 

0 

2 

0 

10825.9638 

23 

2325 

1250 

2100 

3000 

0 

0 

0 

-2 

7071.8846 

24 

2325 

1250 

2100 

15000 

0 

0 

0 

2 

11239.2052 

25 

2325 

1250 

2100 

9000 

0 

0 

0 

0 

9999.9361 

B  -  14 


TABLE  B. 14 


HYBRID  (416A+CP)  EXPERIMENTAL  DESIGN  FOR  FOUR  FACTORS 
NON-CODED  CODED 


RUN# 

Ml 

M2 

B1 

SI 

Ml 

M2 

B1 

SI 

RESPONSE 

1 

2325 

1250 

2100 

14353 

0 

0 

0 

1 . 7844 

11118. 1515 

2 

2325 

1250 

2100 

4517 

0 

0 

0  - 

1.4945 

8057.9946 

3 

1650 

500 

1200 

10932 

-1 

-1 

-1 

.644 

9458.6389 

4 

3000 

500 

1200 

10932 

1 

-1 

-1 

.644 

9918. 1789 

5 

1650 

2000 

1200 

10932 

-1 

1 

-1 

.644 

10084.3805 

6 

3000 

2000 

1200 

10932 

1 

1 

-1 

.644 

10523.0236 

7 

1650 

500 

3000 

10932 

-1 

-1 

1 

.644 

10327.2820 

8 

3000 

500 

3000 

10932 

1 

-1 

1 

.644 

10765.2599 

9 

1650 

8 

8 

8 

CM 

3000 

10932 

-1 

1 

.644 

10901.3249 

10 

3000 

2000 

3000 

10932 

1 

1 

1 

.644 

11290.6917 

11 

3463 

1250 

2100 

6278 

1 . 6853 

0 

0 

-.9075 

9716.9823 

12 

1187 

1250 

2100 

6278  - 

1.6853 

0 

0 

-.9075 

8333. 0546 

13 

2325 

2514 

2100 

6278 

0  1 

.6853 

0 

-.9075 

9887.5823 

14 

2325 

0 

2100 

6278 

0  -1 

.  6853 

0 

-.9075 

8149.0546 

15 

2325 

1250 

3617 

6278 

0 

0  1. 

6853 

-.9075 

10090.3072 

16 

2325 

1250 

583 

6278 

0 

0-1. 

6853 

-.9075 

7851.5946 

17 

2325 

1250 

2100 

9000 

0 

0 

0 

0 

9999.9361 

B  -  15 


TABLE  B.15 


HYBRID  (4160  EXPERIMENTAL  DESIGN  FOR  FOUR  FACTORS 
NON-CODED  CODED 


RUN# 

Ml 

M2 

B1 

SI 

Ml  M2 

B1 

SI 

RESPONSE 

1 

2325 

1250 

2100 

14296 

0  0 

0 

1 . 7654 

11107.4868 

2 

2325 

1250 

2100 

9000 

0  0 

0 

0 

9999.9361 

3 

1650 

500 

1200 

10702 

1 

1 

-1 

.5675 

9405.6469 

4 

3000 

500 

1200 

10702 

1  -1 

-1 

.5675 

9865. 1869 

5 

1650 

2000 

1200 

10702 

-1  1 

-1 

.5675 

1 0033 . 6885 

6 

3000 

2000 

1200 

10702 

1  1 

-1 

.5675 

10474.6316 

7 

1650 

500 

3000 

10702 

-1  -1 

1 

.5675 

10274.2900 

8 

3000 

500 

3000 

10702 

1  -1 

1 

.5675 

10718.6746 

9 

1650 

to 

8 

8 

8 

3000 

10702 

-1  1 

1 

.5675 

10855.9919 

10 

3000 

2000 

3000 

10702 

1  1 

1 

.5675 

11255.6627 

11 

3317 

1250 

2100 

5847 

1 . 4697  0 

0 

-1.0509 

9435.5023 

12 

1383 

1250 

2100 

5847 

-1.4697  0 

0 

-1.0509 

8166.7146 

13 

2325 

2352 

2100 

5847 

0  1 . 4697 

0 

-1.0509 

9586.3623 

14 

2325 

148 

2100 

5847 

0  -1.4697 

0 

-1.0509 

7996.0546 

15 

2325 

1250 

3423 

5847 

0  0  1. 

4697 

-1.0509 

9796. 1623 

16  2325  1250  111  5847 


0 


0  -1.4697  -1.0509  7746.4746 


TABLE  B.16 


BOI-BEHNKEN  EIPERIHENTAL  DESIBN  FOR  FIVE  FACTORS 
NON-CODED  CODED 


RUN  i 

Hi 

n2 

B1 

B2 

SI 

HI 

H2 

B1 

B2 

SI 

RESPONSE 

1 

3111 

2111 

21BB 

3BBB 

9BN 

1 

1 

B 

B 

B 

14743.1514 

2 

3111 

511 

2in 

3BBB 

9BBB 

1 

-1 

B 

B 

B 

13768.1514 

3 

1651 

2111 

21BB 

3BN 

9BBB 

-1 

1 

B 

B 

B 

13973.6514 

4 

1651 

SH 

21BB 

3BN 

9BN 

-1 

-1 

B 

B 

B 

1288B.7422 

5 

2325 

1251 

3BN 

4BBB 

9BH 

B 

B 

1 

1 

B 

1524B.9B14 

& 

2325 

1251 

3BBB 

2BN 

9SBB 

B 

B 

1 

-1 

B 

1376B.9B14 

7 

2325 

1251 

12BB 

4BBB 

9BBB 

B 

B 

-1 

1 

B 

1398B.9B14 

8 

2325 

1251 

12BB 

2BBB 

9BBB 

B 

B 

-1 

-1 

B 

12385.4922 

9 

2325 

2111 

21BB 

3BBI 

12BN 

B 

1 

B 

B 

1 

15379.2168 

IB 

2325 

2111 

21BB 

3BBB 

6BBB 

B 

1 

B 

B 

-1 

12697.4922 

11 

2325 

SBI 

21BB 

3BBB 

12BBB 

B 

-1 

B 

B 

1 

14748.6172 

12 

2325 

SIB 

21BB 

3BBB 

6BBB 

B 

‘1 

B 

B 

-1 

11452.4922 

13 

3BII 

125B 

3BBB 

3BN 

9BBB 

1 

B 

1 

B 

B 

14885.6514 

14 

3BIB 

12SB 

12BB 

3BBB 

9BBB 

1 

B 

-1 

B 

B 

13625.6514 

IS 

165B 

12SB 

3BBB 

3BN 

9BBB 

-1 

B 

1 

B 

B 

14116.1514 

16 

1651 

125B 

12BB 

3BBB 

9BBB 

-1 

B 

-I 

B 

B 

12719.2422 

17 

2325 

12SB 

21BB 

4BN 

12BBB 

B 

B 

B 

1 

1 

15574.3164 

IB 

2325 

12SB 

21BB 

4BBB 

6NB 

B 

B 

B 

1 

-1 

12994.9922 

8  -  17 


TAn.E  B.16  (continued) 


NON-CDDEO  CODED 


ni 

N2 

B1 

B2 

SI 

HI 

H2 

Bl 

B2 

SI 

RESPONSE 

2325 

1251 

21BB 

2BN 

12BBB 

8 

8 

8 

-1 

1 

14518.9814 

2325 

1251 

21M 

2111 

6111 

8 

8 

8 

-1 

-1 

11154.9922 

2325 

2IN 

3NI 

3IN 

9111 

8 

1 

1 

8 

8 

14988.4814 

2325 

2aN 

1211 

3111 

9BN 

8 

1 

-1 

8 

8 

13728.4814 

2325 

5N 

3HI 

3BII 

9Bfl 

8 

-1 

1 

8 

8 

14813.4814 

2325 

5N 

12N 

3IN 

9111 

8 

-1 

-1 

8 

8 

12682.9922 

3ni 

1251 

2111 

4SII 

9881 

1 

8 

8 

1 

8 

14995.6514 

3111 

1251 

2111 

2Bn 

9888 

1 

8 

8 

-1 

8 

13515.6514 

1651 

1251 

2111 

4IN 

9818 

-t 

8 

8 

1 

8 

14226.1514 

1651 

1251 

21N 

2111 

9BN 

-1 

8 

8 

-1 

8 

12591.2422 

2325 

1251 

3111 

3IM 

12888 

8 

8 

1 

8 

1 

15487.2764 

2325 

1251 

3IM 

3811 

6H8 

8 

8 

1 

8 

-1 

12867.7422 

2325 

1251 

1211 

3BII 

12888 

8 

8 

-1 

8 

1 

14628.9114 

2325 

1251 

i2n 

3111 

6888 

8 

8 

-1 

8 

-1 

11273.9922 

3NI 

1251 

2111 

3111 

12888 

1 

8 

8 

8 

1 

15293.6865 

3BII 

1251 

2111 

3111 

68N 

1 

8 

8 

8 

-1 

12588.7422 

1651 

1251 

2111 

3111 

12888 

-1 

8 

8 

8 

1 

14B34.1465 

1651 

1251 

2111 

3IN 

68H 

-1 

8 

8 

8 

-1 

11555.2422 

2325 

2IN 

2111 

4NI 

98H 

8 

1 

8 

1 

8 

15898.4814 

2325 

2IN 

21BI 

2IN 

9888 

8 

1 

8 

-1 

8 

13618.4814 

2325 

5N 

21BI 

4IN 

9888 

8 

-1 

8 

1 

8 

14123.4814 

2325 

5H 

2111 

2111 

9N8 

8 

-1 

8 

-1 

8 

12474.9922 

2325 

1251 

21N 

3IN 

98N 

8 

8 

8 

8 

8 

13878.9814 

TABLE  B.17 

BOX-DRAPER  EXPOtlHENTAL  DESI6N  FOR  FIVE  FACTORS 


NON-CODED 


HI 

82 

Bi 

B2 

SI 

81 

82 

Bl 

B2 

SI 

RESPOUK 

2610 

1564 

2477 

3419 

18257 

.419 

.419 

.419 

.419 

.419 

15281.8667 

2618 

1564 

2477 

3419 

7743 

.419 

.419 

.419 

.419 

-.419 

14238.5615 

2618 

1564 

2477 

2581 

18257 

.419 

.419 

.419 

-.419 

.419 

14768.3715 

2688 

1564 

2477 

2581 

7743 

.419 

.419 

.419 

-.419 

-.419 

13618.4415 

2688 

1564 

1723 

3419 

18257 

.419 

.419 

-.419 

.419 

.419 

14868.6915 

2688 

1564 

1723 

3419 

7743 

.419 

.419 

-.419 

.419 

-.419 

13702.7615 

2688 

1564 

1723 

2581 

18257 

.419 

.419 

-.419 

-.419 

.419 

14248.5715 

2688 

1564 

1723 

2581 

7743 

.419 

.419 

-.419 

-.419 

-.419 

12967.1524 

2688 

936 

2477 

3419 

18257 

.419 

-.419 

.419 

.419 

.419 

14988.2915 

2688 

936 

2477 

3419 

7743 

.419 

-.419 

.419 

.419 

-.419 

13822.3615 

2688 

936 

2477 

2581 

18257 

.419 

-.419 

.419 

-.419 

.419 

14360.1715 

2688 

936 

2477 

2581 

7743 

.419 

-.419 

.419 

-.419 

-.419 

13189.4324 

2688 

936 

1723 

3419 

18257 

.419 

-.419 

-.419 

.419 

.419 

14452.4915 

2688 

936 

1723 

3419 

7743 

.419 

-.419 

-.419 

.419 

-.419 

13216.8724 

2618 

936 

1723 

2581 

18257 

.419 

-.419 

-.419 

-.419 

.419 

13832.3715 

2618 

936 

1723 

2581 

7743 

.419 

-.419 

-.419 

-.419 

-.419 

12445.9124 

2842 

1564 

2477 

3419 

18257 

-.419 

.419 

.419 

.419 

.419 

15865.8715 

2842 

1564 

2477 

3419 

7743 

-.419 

.419 

.419 

.419 

-.419 

13902.2815 

2842 

1564 

2477 

2581 

11257 

-.419 

.419 

.419 

-.419 

.419 

14445.7515 

2842 

1564 

2477 

2581 

7743 

-.419 

.419 

.419 

-.419 

-.419 

13288.5124 
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TABLE  B.t8 


BOX-DRAPER  (HALF-REPLICATE)  EXPERIHENTAL  DESIGN  FOR  FIVE  FACTORS 


NOH-CODED 


HI 

n2 

Dl 

B2 

SI 

HI 

H2 

B1 

B2 

SI 

RESPONSE 

2159 

955 

1745 

2616 

7B1B 

-.394 

-.394 

-.394 

-.394 

-.394 

12135.5524 

2591 

1545 

1745 

2616 

7B1B 

.394 

.394 

-.394 

-.394 

-.394 

13129.5724 

2591 

955 

2455 

2616 

7B1B 

.394 

-.394 

.394 

-.394 

-.394 

13164.6724 

2159 

1545 

2455 

2616 

7B1B 

-.394 

.394 

.394 

-.394 

-.394 

13251.1524 

2591 

955 

1745 

3394 

7B1B 

.394 

-.394 

-.394 

.394 

-.394 

13264.8324 

2159 

1545 

1745 

3394 

7B1B 

-.394 

.394 

-.394 

.394 

-.394 

13350.2124 

2159 

955 

2455 

3394 

7B1B 

-.394 

-.394 

.394 

.394 

-.394 

13485.3124 

2591 

1545 

2455 

3394 

7B1B 

.394 

.394 

.394 

.394 

-.394 

14209.7015 

2591 

955 

1745 

2616 

1I1B2 

.394 

-.394 

-.394 

-.394 

.394 

13B34.4315 

2159 

1545 

1745 

2616 

101B2 

-.394 

.394 

-.394 

-.394 

.394 

13914.6915 

2B59 

955 

2455 

2616 

1I1B2 

-.394 

-.394 

.394 

-.394 

.394 

14028.1915 

2591 

1545 

2455 

2616 

1I1B2 

.394 

.394 

.394 

-.394 

.394 

14714.9315 

2B59 

955 

1745 

3394 

1I1B2 

-.394 

-.394 

-.394 

.394 

.394 

14114.3115 

2591 

1545 

1745 

3394 

1I1B2 

.394 

.394 

-.394 

.394 

.394 

14801.0515 

2591 

955 

2455 

3394 

1I1B2 

.394 

-.394 

.394 

.394 

.394 

14914. 5515 

2059 

1545 

2455 

3394 

1I1B2 

-.394 

.394 

.394 

.394 

.394 

14994.8115 

•,  4“  in 


TABLE  B.19 


HINIHUIt  POINT  EXPERINENTAL  DESIGN  FOR  FIVE  FACTORS 

NON-COOED  CODED 


RUN  i 

HI 

N2 

81 

82 

SI 

HI 

H2 

81 

82 

SI 

RESPONSE 

1 

1651 

511 

12N 

28N 

68N 

-1 

-1 

-1 

-1 

-1 

92N.2424 

2 

3111 

SN 

1211 

28N 

6N8 

1 

-1 

-1 

-1 

-1 

18251.2424 

3 

1651 

2111 

1211 

28N 

6N8 

-1 

1 

-1 

-1 

-1 

18456.7424 

4 

1651 

SH 

3NI 

28N 

6N8 

-1 

-1 

1 

-1 

-1 

1N12.2424 

5 

1651 

sn 

1211 

48N 

6888 

-1 

-1 

-1 

1 

-1 

11848.2424 

& 

1651 

SN 

1211 

2888 

12N8 

-1 

-1 

-1 

-1 

1 

13888.6515 

7 

1778 

642 

12N 

28N 

6N8 

-.8188 

-.81N 

-1 

-1 

-1 

9419.3624 

8 

1778 

SN 

1371 

28N 

68N 

-.aiN 

-1 

-.8188 

-1 

9453.8824 

9 

1778 

5N 

12N 

2189 

68N 

-.81N 

-1 

-1 

-.BIN 

-1 

9473.9624 

11 

1778 

Sll 

1281 

28N 

6568 

-.8188 

-1 

-1 

-1 

-.BIN 

9688.6424 

11 

1651 

642 

1378 

2888 

6888 

-.8188 

-.8188 

-1 

-1 

9472.5224 

12 

1651 

642 

1288 

2189 

68H 

-.81N 

-1 

-.8188 

9493.4824 

13 

1651 

642 

1288 

28N 

6568 

-1 

-.8188 

-1 

-1 

-.8188 

9788.  N24 

14 

1651 

SN 

1378 

2189 

6888 

-1 

-1 

-.8188 

-.8188 

-1 

9527.1224 

IS 

1651 

SN 

1378 

28N 

6568 

-1 

-1 

-.8188 

-1 

-.8188 

9733.8824 

16 

1651 

5N 

1288 

2189 

6568 

-1 

-1 

-1 

-.BIN 

-.8188 

9754.6824 

17 

2686 

2IN 

3888 

4888 

12N8 

.5355 

1 

1 

1 

1 

16396.5653 

18 

3181 

1652 

38N 

48N 

128H 

1 

.5355 

1 

1 

1 

16364.7917 

19 

3111 

2IN 

2582 

48N 

12N8 

1 

1 

.5355 

1 

1 

16313.8837 

21 

3811 

2IN 

3888 

3535 

12N8 

1 

1 

1 

.5355 

1 

16276.9228 

21 

3111 

2IN 

38N 

48N 

18686 

1 

1 

1 

1 

.5355 

16197.7212 

TABLE  B.21 


I 


CENTRAL  COHPOSITE  ROTATABLE  (HALF-REPLICATE)  EXPERIMENTAL  DESIGN  FOR  FIVE  FACTORS 

NON-CODED  CODED 


N1 

H2 

B1 

B2 

SI 

HI 

H2 

Bi 

B2 

SI 

RESPONSE 

1651 

511 

12N 

28N 

68N 

-I 

-1 

-1 

-1 

-1 

92N.2424 

3BII 

2IN 

12N 

28N 

68N 

1 

1 

-1 

-1 

-1 

11496.2424 

3NI 

SN 

3111 

28N 

6BN 

1 

-1 

1 

-1 

-1 

11838.2424 

1651 

2111 

3IN 

2N8 

68N 

-1 

1 

1 

-1 

-1 

12857.2424 

3BII 

5IB 

1218 

4888 

6N8 

1 

-1 

-1 

1 

-1 

12891.2424 

1651 

2111 

1288 

48N 

6888 

-1 

1 

-1 

1 

-1 

12296.7424 

1651 

SN 

3N8 

48N 

68N 

-1 

1 

1 

-1 

12652.2424 

3NI 

2IH 

3NB 

4SN 

68N 

1 

1 

1 

1 

-1 

1471B.1S1S 

3IU 

SN 

12N 

28N 

128N 

I 

1 

13778.1515 

1651 

2IN 

1288 

2888 

128N 

1 

-1 

-1 

1 

13983.6515 

1651 

581 

3888 

2888 

128N 

-1 

-1 

1 

-1 

1 

1426B.6515 

3BN 

2IN 

3888 

2888 

128N 

1 

1 

1 

-1 

1 

15521.9467 

1651 

SN 

1288 

4888 

12N8 

-1 

-1 

-1 

1 

1 

1448B.6S1S 

3BII 

2BN 

1288 

4888 

128N 

1 

1 

-1 

1 

1 

15694. 2388 

3BH 

5N 

3888 

4N8 

128N 

1 

-1 

1 

1 

1 

15984.3588 

1651 

2IN 

3N8 

4888 

12N8 

.1 

1 

1 

1 

1 

16858.9188 
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TABLE  B.22 


BOI-BEHNKEN  EXPERIHENTAL  DESI6N  FOR  SIX  FACTORS 


NON-CODED  CODED 


RIM  t 

HI 

N2 

N3 

B1 

B2 

SI 

Ml 

112 

N3 

B1 

B2 

SI 

RESPONSE 

1 

2325 

2IN 

711 

2111 

4IN 

12NI 

1 

1 

1 

1 

1 

1 

16271.1111 

2 

2325 

2111 

711 

2111 

4NI 

6NI 

1 

1 

1 

1 

1 

-1 

14213. S6B2 

3 

2325 

2IH 

7N 

21N 

2IN 

12NI 

1 

1 

1 

1 

-1 

1 

15246. SSII 

4 

2323 

2NI 

711 

21N 

2111 

6NI 

1 

1 

1 

1 

-1 

-1 

12417.1591 

S 

2325 

Sll 

711 

2111 

4IN 

12IH 

1 

-1 

1 

1 

1 

1 

15644. 75N 

b 

2325 

SIB 

7N 

2111 

4IN 

6111 

1 

-1 

1 

1 

1 

-1 

11172.1591 

7 

2325 

SBI 

711 

2111 

2IN 

12IN 

1 

-1 

1 

1 

-1 

1 

14531.1682 

B 

2325 

Sll 

7N 

2111 

2HI 

6NI 

1 

-1 

1 

1 

-1 

-1 

13112.1591 

9 

3111 

1251 

nil 

2111 

3111 

12NI 

1 

1 

1 

1 

1 

1 

15817.5411 

11 

3aH 

1251 

IIH 

2111 

3111 

6111 

1 

1 

1 

1 

I 

-1 

13489.4191 

11 

3111 

1251 

4H 

2111 

3111 

12III 

1 

1 

-1 

1 

1 

1 

15541. 3BBI 

12 

3aN 

1251 

411 

2111 

3IH 

6111 

1 

1 

-1 

1 

1 

12967.4191 

13 

1651 

1251 

IIN 

2111 

3UI 

12NI 

-1 

1 

1 

1 

1 

1 

15358.  NN 

14 

1651 

1251 

IIN 

21N 

3IN 

6IH 

-1 

1 

1 

1 

1 

-1 

12455.9191 

IS 

1651 

1251 

411 

21N 

3IN 

12NI 

-1 

1 

-1 

1 

1 

1 

11933.9191 

16 

1651 

1251 

411 

21N 

3MI 

6NI 

-1 

1 

-1 

1 

1 

1SIS8.NN 

17 

3BN 

2111 

4N 

3NI 

3IN 

9NI 

1 

1 

1 

1 

1 

1 

15715.  INI 

18 

3111 

2111 

411 

1211 

3IN 

9111 

1 

1 

I 

1 

1 

14621.8182 

19 

3IBI 

Sll 

411 

3111 

3IN 

9111 

1 

-1 

1 

1 

1 

1 

14915.8182 

21 

3IN 

Sll 

411 

1211 

3111 

9111 

1 

1 

-1 

I 

I 

13645.8182 

21 

1651 

2111 

411 

3111 

3IN 

9111 

-1 

1 

1 

1 

1 

1 

15111.3182 

22 

1651 

2111 

411 

1211 

3111 

9NI 

-1 

1 

1 

-1 

1 

1 

13851.3182 

23 

1651 

Sll 

4N 

3NI 

3111 

9NI 

-1 

1 

1 

I 

1 

14136.3182 

24 

1651 

Sll 

411 

1211 

3111 

9NI 

-1 

1 

-1 

1 

1 

12736.4191 

23 

2325 

2111 

IIN 

2111 

4111 

9111 

1 

1 

1 

1 

1 

1 

15699.SIN 

26 

2325 

2111 

IIN 

2111 

2111 

9NI 

1 

1 

1 

1 

-1 

1 

14333.0682 

27 

2325 

2HI 

411 

2111 

4IN 

9IN 

1 

1 

-I 

1 

1 

1 

15399.1682 

2B 

2323 

2111 

411 

2111 

2111 

9111 

1 

1 

-I 

1 

-1 

1 

13919.1682 

29 

2323 

Sll 

nil 

2111 

4HI 

9NI 

1 

-1 

1 

1 

1 

1 

14838.1682 

31 

2325 

Sll 

nil 

2111 

2111 

9111 

1 

-1 

1 

1 

-1 

1 

13358.1682 

31 

2323 

511 

411 

2111 

4111 

9IN 

1 

-1 

-1 

1 

1 

1 

14424.1682 

32 

2325 

Sll 

411 

2111 

2111 

9111 

1 

-1 

1 

1 

12853.6591 

TABLE  B.22  (continued) 


NON-CODED  CODED 


RUN  t 

N1 

n2 

N3 

B1 

B2 

SI 

HI 

N2 

H3 

B1 

B2 

SI 

RESPONSE 

33 

2325 

1258 

1888 

3888 

3888 

12888 

8 

8 

1 

8 

1 

16811.1388 

34 

2325 

1258 

1888 

3888 

3888 

6888 

8 

8 

1 

1 

8 

-1 

13768.4891 

35 

2325 

1258 

1888 

1288 

3888 

12888 

8 

8 

1 

-1 

8 

1 

15154.8888 

36 

2325 

1258 

1881 

1288 

3888 

6n8 

8 

8 

1 

-1 

8 

-1 

12174.6591 

37 

2325 

1258 

488 

3888 

3888 

12888 

8 

8 

-1 

1 

8 

1 

15734.B9N 

38 

2325 

1258 

488 

3888 

3888 

68N 

8 

8 

-1 

1 

8 

-1 

13246.4891 

39 

2325 

1258 

488 

1288 

38N 

128H 

8 

8 

-1 

-1 

8 

1 

14855. 58N 

41 

2325 

1258 

488 

1288 

3888 

68N 

8 

8 

-1 

-1 

8 

-1 

11652.6591 

41 

3111 

1258 

788 

3888 

48H 

9888 

1 

8 

8 

1 

1 

8 

15928.8888 

42 

3881 

1258 

788 

3888 

2888 

9888 

1 

8 

8 

1 

-1 

8 

14653.3182 

43 

3888 

1258 

788 

1288 

4888 

9888 

1 

8 

8 

-1 

1 

8 

14873.3182 

44 

3888 

1258 

788 

1288 

2888 

9888 

1 

8 

8 

-1 

-1 

8 

13393.3182 

45 

1658 

1258 

788 

3888 

4888 

9888 

-1 

8 

8 

1 

1 

8 

15363.8182 

46 

1658 

1258 

788 

3888 

2888 

9888 

-1 

8 

8 

1 

-1 

8 

13883.8182 

47 

1658 

1258 

788 

1288 

4888 

9888 

-1 

8 

8 

-1 

1 

8 

14183.8182 

48 

1658 

1258 

7N 

1288 

28H 

9N8 

-1 

8 

8 

-1 

-1 

8 

12438.9891 

TABLE  B.23 


BOX-DRAPER  (HALF-REPLICATE)  EXPERIHENTAL  DESI6N  FOR  SIX  FACTORS 


NON-CDDED 


ni 

N2 

N3 

B1 

B2 

SI 

HI 

N2 

N3 

B1 

B2 

SI 

RESPONSE 

2142 

936 

574 

1723 

2581 

18257 

-.419 

-.419 

-.419 

-.419 

-.419 

.419 

13938.5553 

2142 

936 

574 

1723 

3419 

7743 

-.419 

-.419 

-.419 

-.419 

.419 

-.419 

13316.8614 

2B42 

936 

574 

2477 

2581 

7743 

-.419 

-.419 

-.419 

.419 

-.419 

-.419 

13289.5974 

2149 

936 

574 

2477 

3419 

18257 

-.419 

-.419 

-.419 

.419 

.419 

.419 

15878.6153 

2142 

936 

826 

1723 

2581 

7743 

-.419 

-.419 

.419 

-.419 

-.419 

-.419 

12764.6194 

2142 

936 

826 

1723 

3419 

18257 

-.419 

-.419 

.419 

-.419 

.419 

.419 

14724.1413 

2842 

936 

826 

2477 

2581 

18257 

-.419 

-.419 

.419 

.419 

-.419 

.419 

14631.9613 

2842 

936 

826 

2477 

3419 

7743 

-.419 

-.419 

.419 

.419 

.419 

-.419 

14888.4931 

2842 

1564 

574 

1723 

2581 

7743 

-.419 

.419 

-.419 

-.419 

-.419 

-.419 

13867.5564 

2842 

1564 

574 

1723 

3419 

18257 

-.419 

.419 

-.419 

-.419 

.419 

.419 

14959.2883 

2842 

1564 

574 

2477 

2581 

18257 

-.419 

.419 

-.419 

.419 

-.419 

.419 

14867.8283 

2842 

1564 

574 

2477 

3419 

7743 

-.419 

.419 

-.419 

.419 

.419 

-.419 

14323.5521 

2842 

1564 

826 

1723 

2581 

18257 

-.419 

.419 

.419 

-.419 

-.419 

.419 

14512.5463 

2842 

1564 

826 

1723 

3419 

7743 

-.419 

.419 

.419 

-.419 

.419 

-.419 

13969.8781 

2842 

1564 

826 

2477 

2581 

7743 

-.419 

.419 

.419 

.419 

-.419 

-.419 

13876.8981 

2842 

1564 

826 

2477 

3419 

18257 

-.419 

.419 

.419 

.419 

.419 

.419 

15516.5944 

2A88 

936 

574 

1723 

2581 

7743 

.419 

-.419 

-.419 

-.419 

-.419 

-.419 

12975.7954 

2688 

936 

574 

1723 

3419 

18257 

.419 

-.419 

-.419 

-.419 

.419 

.419 

14873.8958 

2688 

936 

574 

2477 

2581 

18257 

.419 

-.419 

-.419 

.419 

-.419 

.419 

14788.9158 

2688 

936 

574 

2477 

3419 

7743 

.419 

-.419 

-.419 

.419 

.419 

-.419 

14243.1841 

2688 

936 

826 

1723 

2581 

18257 

.419 

-.419 

.419 

-.419 

-.419 

.419 

14426.4418 

2688 

936 

826 

1723 

3419 

7743 

.419 

-.419 

.419 

-.419 

.419 

-.419 

13888.6381 

2688 

936 

826 

2477 

2581 

7743 

.419 

-.419 

.419 

.419 

-.419 

-.419 

13796.4581 

2688 

936 

826 

2477 

3419 

18257 

.419 

-.419 

,419 

.419 

.419 

.419 

15442.4314 

2686 

1564 

574 

1723 

2581 

18257 

.419 

.419 

-.419 

-.419 

-.419 

.419 

14661.5888 

2688 

1564 

574 

1723 

3419 

7743 

.419 

.419 

-.419 

-.419 

.419 

-.419 

14123.6891 

2688 

1564 

574 

2477 

2581 

7743 

.419 

.419 

-.419 

.419 

-.419 

-.419 

14831.5891 

2688 

1564 

574 

2477 

3419 

18257 

.419 

.419 

-.419 

.419 

.419 

.419 

15685.8414 

2688 

1564 

826 

1723 

2581 

7743 

.419 

.419 

.419 

-.419 

-.419 

-.419 

13677.8351 

2688 

1564 

826 

1723 

3419 

18257 

.419 

.419 

.419 

-.419 

.419 

.419 

15346.8994 

2688 

1564 

826 

2477 

2581 

18257 

.419 

.419 

.419 

.419 

-.419 

.419 

15278.6414 

2688 

1564 

826 

2477 

3419 

7743 

.419 

.419 

.419 

.419 

.419 

-.419 

14825.8951 

TABLE  B.23  (continued) 


HON-CODEB 


CODED 


RUN  t 

HI 

N2 

(13 

B1 

B2 

SI 

HI 

n2 

n3 

B1 

B2 

SI 

RESPONSE 

33 

2998 

1251 

788 

2188 

3888 

9888 

.997 

8 

8 

8 

8 

8 

14762.1668 

34 

1652 

1231 

718 

2188 

3888 

9N8 

-.997 

8 

8 

8 

8 

8 

13994.9753 

33 

2323 

1998 

718 

2188 

38N 

9888 

8 

.997 

8 

8 

8 

8 

14864.6857 

36 

2323 

582 

7H 

2188 

38N 

9888 

8 

-.997 

8 

8 

8 

8 

13892.5387 

37 

2325 

1258 

999 

2188 

3888 

9N8 

8 

8 

.997 

8 

8 

8 

14584.9472 

3B 

2323 

1251 

481 

2188 

3888 

9888 

8 

8 

-.997 

8 

8 

8 

14172.1892 

39 

2325 

1258 

788 

2997 

38n 

9888 

8 

8 

8 

.997 

8 

8 

15886.6782 

41 

2325 

1258 

788 

1283 

3888 

9888 

8 

8 

8 

-.997 

8 

8 

13758.4582 

41 

2325 

1251 

788 

2188 

3997 

9888 

8 

8 

8 

8 

.997 

8 

15116.3482 

42 

2325 

1258 

788 

2188 

2883 

9n8 

8 

8 

8 

8 

-.997 

8 

13648.7882 

43 

2325 

1251 

788 

2188 

3888 

11991 

8 

8 

8 

8 

8 

.997 

15447.5764 

44 

2325 

1258 

788 

2188 

3888 

6889 

8 

8 

8 

8 

8 

-.997 

12728.5991 

45 

2325 

1258 

788 

2188 

3888 

9888 

8 

8 

8 

8 

8 

8 

14378.5682 

TABLE  B.24 

HYBRID  (&28A)  EXPERIHENTAL  DESIGN  FDR  SIX  FACTORS 


NON-CODED 


N2 

N3 

B1 

B2 

SI 

12SI 

7BI 

2iBB 

3BN 

1&1I3 

12SI 

7BB 

21BB 

3BBB 

35A7 

SN 

4BI 

12BB 

2BBB 

lBB2f 

2111 

4N 

12BB 

2BBB 

1B829 

SBI 

IBBB 

12BB 

2BN 

tBS29 

2NI 

IBBB 

12BB 

2BBB 

1B829 

Sll 

4Ba 

3BH 

2BN 

18829 

2ni 

4BI 

SBBB 

2IBB 

1B829 

12SB 

IBBB 

3BIB 

2flBfl 

18829 

2BBB 

IBBB 

3BBB 

2BBB 

1BB29 

SBB 

4BB 

12BB 

4BBB 

1B829 

2BBB 

4BB 

12BB 

4BN 

1B829 

SBB 

IBBB 

12BB 

4BBB 

1B829 

2BBB 

IBBB 

12BB 

4BBB 

1BB29 

SBB 

4BB 

3BBB 

4BBB 

1B829 

2BBB 

4BB 

3BBB 

4BBB 

18829 

SBB 

IBBB 

3BBB 

4BBB 

18829 

2BBB 

IBBB 

3BBB 

4BBB 

1B829 

B  -  31 


CODED 
N3  B1 

•  I 

I  I 


B2 

SI 

RESPONSE 

8 

2.3677 

16388.4422 

B 

-1.8118 

11B96.2991 

-1 

.6896 

1272S.B411 

-1 

.6896 

14515. B662 

-1 

.6896 

13954. B662 

-1 

.6896 

14159.5662 

-1 

.6B96 

14881. B662 

-1 

.6B96 

14995. 776B 

-1 

.6B96 

14444.5662 

-1 

.6B96 

15775.9555 

1 

.6B96 

15BBB.7766 

1 

.6896 

15177. 776B 

1 

.6B96 

14664.5662 

1 

.6B96 

1595B.B355 

1 

.6896 

15485.7768 

1 

.6896 

16589.2275 

1 

.6896 

16166.1555 

1 

.6896 

16327.4275 

TABLE  B.24  (continued) 


NON-CODED  CODED 


n2 

H3 

B1 

B2 

SI 

HI 

N2 

N3 

B1 

B2 

SI 

RESPONSE 

1251 

711 

2181 

3888 

5987 

2.1749 

8 

8 

8 

8 

-1.831 

13778.6847 

125B 

7BI 

2111 

3888 

5987 

-2.1749 

8 

8 

8 

8 

-1.831 

11519.8122 

2BB1 

7N 

2181 

3888 

5987 

8 

2.1749 

8 

8 

8 

-1.831 

13991.6952 

1 

711 

2188 

3888 

5987 

8 

■2.1749 

8 

8 

8 

-1.831 

11614.8991 

125B 

1352 

2188 

3n8 

5987 

8 

8 

2.1749 

8 

8 

-1.831 

13228.9288 

12SB 

48 

2188 

3888 

5987 

8 

8  • 

•2.1749 

8 

8 

-1.831 

12885.6382 

12SB 

7BI 

4857 

3888 

5987 

8 

8 

8 

2.1749 

8 

-1.831 

14383.2952 

12SB 

7BI 

143 

3888 

5987 

8 

8 

8  ■ 

■2.1749 

8 

-1.831 

18982.4381 

12SB 

711 

2188 

5175 

5987 

8 

8 

8 

8 

2.1749 

-1.831 

14548.8542 

1251 

7BI 

2188 

825 

5987 

8 

8 

8 

8  - 

■2.1749 

-1.831 

18652.3711 

TABLE  B.2S 


HYBRID  (628B)  EXPERIHENTAL  DESI6N  FOR  SIX  FACTORS 


NON-CODED  CODED 


RUN  t 

Hi 

N2 

H3 

B1 

B2 

SI 

HI  H2 

N3 

B1 

B2 

SI 

RESPONSE 

1 

232S 

1258 

7N 

2188 

38N 

15928 

8  8 

8 

8 

8 

2.3894 

16343.6433 

2 

2325 

1258 

788 

2188 

3888 

9881 

8  8 

8 

8 

8 

8 

14378.5682 

3 

1651 

SM 

488 

1288 

2N8 

18732 

-1  -1 

-1 

-1 

.5773 

12663.8251 

4 

3881 

2888 

488 

1288 

2888 

18732 

1  1 

-1 

-1 

-1 

.5773 

14478.4922 

5 

3888 

588 

1888 

1288 

2Bn 

18732 

1  -1 

1 

-1 

-1 

.5773 

13989.4922 

6 

1658 

2888 

1881 

1288 

28n 

18732 

-1  1 

1 

-1 

-1 

.5773 

14114.9922 

7 

3818 

588 

4N 

3888 

2888 

18732 

1  -1 

-1 

1 

-1 

.5773 

14755.4922 

8 

1658 

2n8 

481 

38n 

2888 

18732 

1 

-1 

.5773 

14968.9922 

9 

1658 

1258 

1888 

3888 

2888 

18732 

1 

1 

-1 

.5773 

14399.9922 

18 

3888 

2888 

1888 

3888 

2888 

18732 

1  1 

1 

1 

-1 

.5773 

15753.6298 

11 

3888 

588 

481 

1288 

4888 

18732 

“1 

1 

.5773 

14974.6138 

12 

1658 

2888 

488 

1288 

4888 

18732 

-1  1 

-1 

-1 

1 

.5773 

15151.6138 

13 

1658 

588 

1888 

1288 

4888 

18732 

-1  -1 

1 

-1 

1 

.5773 

14619.9922 

14 

3888 

2888 

1888 

1288 

4888 

18732 

1  1 

1 

-1 

1 

.5773 

15927.7898 

IS 

1658 

588 

488 

3888 

4888 

18732 

-1  -1 

-1 

1 

1 

.5773 

15379.6138 

16 

3888 

2888 

488 

3888 

4888 

18732 

1  1 

-1 

1 

1 

.5773 

16487.8788 

17 

3888 

588 

1888 

3888 

4888 

18732 

1  -1 

1 

1 

1 

.5773 

16143.8298 

18 

1658 

2888 

1888 

3888 

4888 

18732 

~1  1 

1 

1 

1 

.5773 

16386.8788 

TABLE  B.26 


17 


niNinUH  POINT  EXPERIHENTAL  DESI6N  FOR  SIX  FACTORS 


NON-CODED 


RtW  I  N1 


N1 

N2 

N3 

81 

B2 

SI 

HI 

H2 

ns  Bi 

B2 

SI 

RESPONSE 

2742 

1714 

4N 

12N 

2888 

6888 

.6183 

.6183 

-1  -1 

-1 

-1 

11438.9145 

2742 

511 

8B5 

1288 

2888 

6881 

.6183 

-1 

.6183  -1 

-1 

-1 

18853.3184 

2742 

511 

411 

2656 

28N 

68N 

.6183 

-1 

-1  .6183 

-1 

-1 

11718.7887 

2742 

511 

411 

12M 

3618 

68N 

.6183 

-1 

-1  -1 

.6183 

-1 

11919.7781 

2742 

5fll 

4N 

1288 

28N 

18855 

.6183 

-1 

-1  -1 

-1 

.6183 

13485.2117 

i6sa 

1714 

885 

1288 

2888 

6888 

-1 

.6183 

.6183  -1 

-1 

-1 

11828.1813 

16SI 

1714 

481 

2656 

2888 

6888 

-1 

.6183 

-1  .6183 

-1 

-1 

11894.8633 

16SI 

1714 

488 

12N 

3618 

6881 

-1 

.6183 

-1  -1 

.6183 

-1 

12886.6418 

1651 

1714 

488 

1288 

28M 

18855 

.6183 

-1  -1 

-1 

.6183 

13571.4967 

1651 

511 

885 

2656 

28H 

6M8 

-1 

.6183  .6183 

-1 

-1 

11389.8437 

1651 

511 

885 

1288 

3618 

68N 

-1 

.6183  -1 

.6183 

-1 

11498.1214 

1651 

511 

885 

1288 

2888 

18^5 

-1 

-1 

.6183  -1 

-1 

.6183  ■ 

13117.5683 

1651 

511 

4N 

2656 

3618 

6N8 

-1 

-1 

-1  .6183 

.6183 

-1 

12375.5834 

1651 

511 

488 

2656 

2H8 

18885 

-1 

-1 

-1  .6183 

-1 

.6183 

13882.1812 

1651 

511 

488 

12N 

3618 

18885 

-1 

-1 

-1  -1 

.6183 

.6183 

13988.1142 

1727 

2NI 

18N 

3888 

4888 

12N8 

-.8854 

1 

1  1 

1 

1 

16611.1197 

3111 

SB6 

18H 

3888 

4888 

12888 

1 

-.8854 

1  1 

1 

1 

16465.4144 

3BN 

2BII 

434 

3888 

4888 

128n 

1 

1 

-.8854  1 

1 

1 

16782.5889 

3IN 

2NI 

1888 

1383 

4888 

12N8 

1 

1 

1  -.8854 

1 

1 

16267.3343 

3IN 

2IM 

1888 

3N8 

2115 

12881 

1 

1 

1  1 

-.8854 

1 

16184.2918 

TABLE  B.2&  (continued) 


NON-CODED  CODED 


RUN  t 

N1 

R2 

N3 

B1 

B2 

Si 

N1 

N2 

N3 

B1 

B2 

SI 

RESPONSE 

2! 

3NI 

2111 

IIM 

3IN 

4111 

6344 

1 

1 

1 

1 

1 

-.BB54 

15594.4142 

22 

3HI 

511 

4N 

1211 

2IN 

6NI 

1 

-1 

-1 

-1 

-1 

-1 

11629.9191 

23 

1651 

2ni 

411 

1211 

2IN 

6BII 

-1 

1 

-1 

-1 

-1 

-1 

1IB35.4I91 

24 

1651 

5N 

INI 

12M 

2IN 

6NI 

-1 

-1 

1 

-1 

-1 

-1 

1B1BI.9I91 

2S 

1651 

511 

4n 

3111 

2IN 

6NI 

-1 

-1 

-1 

1 

-1 

-1 

lil9B.9l91 

26 

1651 

sn 

4N 

1211 

4IN 

6MI 

-1 

-1 

-1 

-1 

1 

-1 

li41B.9B91 

27 

1651 

511 

4M 

12N 

2IN 

12IM 

-1 

-1 

-1 

-1 

-1 

1 

13319. 31B2 

2B 

1651 

511 

4M 

1211 

2IN 

6NI 

-I 

-1 

-1 

957B.9I91 

B  -  36 


Urn 


r  -  -  'j 


‘"..A, 


•r»  ^ ,  • 


TABLE  B.Z7 


KOSHAL  EXPERlflENTAL  DESI6N  FDR  SH  FACTORS 


NON-CODEO 


HI 

HZ 

N3 

B1 

B2 

SI 

HI 

2887 

1784 

414 

1843 

2714 

8143 

.7143 

2887 

1834 

914 

1843 

2714 

8143 

.7143 

2887 

1834 

414 

2743 

2714 

8143 

.7143 

2887 

1834 

414 

1843 

3714 

8143 

.7143 

2887 

1834 

414 

1843 

2714 

11143 

.7143 

2132 

1784 

914 

1843 

2714 

8143 

-.2857 

2132 

1784 

414 

2743 

2714 

8143 

-.2857 

2132 

1784 

414 

1843 

3714 

8143 

-.2857 

2132 

1784 

414 

1843 

2714 

11143 

-.2857 

2132 

1836 

914 

2743 

2714 

8143 

-.2857 

2132 

1836 

914 

1843 

3714 

8143 

-.2857 

2132 

1836 

914 

1843 

2714 

11143 

-.2857 

2132 

1834 

414 

2743 

3714 

8143 

-.2857 

2132 

1834 

414 

2743 

2714 

11143 

-.2857 

2132 

1834 

414 

2743 

3714 

11143 

-.2857 

3482 

1836 

414 

1843 

2714 

8143 

1.7143 

2132 

2534 

414 

1843 

2714 

8143 

-.2857 

2132 

1834 

1214 

1843 

2714 

8143 

-.2857 

2132 

1834 

414 

3443 

2714 

8143 

-.2857 

2132 

1834 

414 

1843 

4714 

8143 

-.2857 

CODED 

N3  B1  B2  SI  RESPONSE 
-.2857  -.2857  -.2857  -.2857  1415&.8B33 
.7143  -.2857  -.2857  -.2857  13874.3133 
-.2857  .7143  -.2857  -.2857  14299.3133 
-.2857  -.2857  .7143  -.2857  14419.3833 
-.2857  -.2857  -.2857  .7143  14984.6895 
.7143  -.2857  -.2857  -.2857  13974.8838 
-.2857  .7143  -.2857  -.2857  14399.8838 
-.2857  -.2857  .7143  -.2857  14519.8138 
-.2857  -.2857  -.2857  .7143  15175.1195 
.7143  .7143  -.2857  -.2857  14119.3838 
.7143  -.2857  .7143  -.2857  14229.3138 
.7143  -.2857  -.2857  .7143  14871.5533 
-.2857  .7143  .7143  -.2857  14452.3838 
-.2857  .7143  -.2857  .7143  15189.1895 
-.2857  -.2857  .7143  .7143  15288.1895 
-.2857  -.2857  -.2857  -.2857  14854.8533 
-.2857  -.2857  -.2857  -.2857  14257.3838 
1.7143  -.2857  -.2857  -.2857  13496.3838 
-.2857  1.7143  -.2857  -.2857  14542.3838 
-.2857  -.2857  1.7143  -.2857  14742.3838 


B  -  37 


TABLE  B.28 


CENTRAL  CDNPOSITE  ROTATABLE  (HALF-REPLICATE)  EXPERINENTAL  DESIGN  FOR  SIX  FACTORS 


NON-COOED 


CODED 


N1 

H2 

N3 

B1 

B2 

SI 

N1  N2 

N3 

B1 

B2 

SI 

RESPONSE 

16SI 

511 

481 

1288 

2888 

128N 

-1 

-1 

•1 

-1 

1 

13389.3182 

1651 

511 

4N 

1288 

4888 

6N8 

-1  -1 

-1 

-1 

1 

-1 

1141B.9891 

1651 

511 

488 

38N 

2888 

68N 

-1  -1 

-1 

1 

-1 

-1 

11198.9891 

1651 

511 

488 

3888 

4888 

12888 

-1  -1 

-1 

1 

1 

1 

15788.2288 

1651 

511 

1888 

1288 

28N 

6888 

-1  -1 

1 

-1 

-1 

-1 

18188.9891 

1651 

511 

1888 

1288 

4888 

128n 

-1  -1 

1 

-1 

1 

1 

15184.  NN 

1651 

511 

18N 

38n 

2888 

12888 

-1  -1 

1 

1 

-1 

1 

14922.8888 

1651 

511 

1888 

3888 

4888 

6888 

-1  -1 

1 

1 

1 

-1 

13552.9891 

1651 

2in 

488 

1288 

2888 

6888 

-1  1 

-1 

-1 

-1 

-1 

18835.4891 

1651 

2NI 

4N 

1288 

4888 

12888 

-1  1 

-1 

-1 

1 

1 

15484. 18n 

1651 

2111 

4H 

3888 

2BN 

12888 

-1  1 

-1 

1 

-1 

1 

15318.8288 

1651 

2111 

488 

3N8 

4888 

6N8 

-1  1 

-1 

1 

1 

-1 

14235.8182 

1651 

2in 

1888 

1288 

2888 

128N 

-1  1 

1 

-1 

-1 

1 

14691.8888 

1651 

2111 

1888 

1288 

4888 

6888 

-1  1 

1 

-1 

1 

-1 

13197.4891 

1651 

2in 

1888 

3888 

2888 

6888 

1 

1 

-1 

12957.9891 

1651 

2fN 

1888 

3818 

4888 

12888 

1 

1 

1 

1 

16585.5688 

3111 

511 

4N 

12N 

28N 

68H 

1  -1 

-1 

18629.9891 

3111 

SIB 

4N 

1288 

4888 

12N8 

-1 

-1 

1 

1 

15313.  B4N 

3111 

Sll 

4N 

3888 

2888 

12888 

-1 

1 

-I 

1 

15135.8888 

3011 

511 

4n 

38H 

48N 

6H8 

1  -1 

-1 

1 

1 

-1 

14843.8182 

3IN 

511 

1888 

12H 

28N 

12888 

1  -1 

1 

-1 

-1 

1 

14492.8182 

3IN 

511 

1888 

1288 

4888 

6888 

1 

-1 

1 

-1 

12991.9891 

3111 

511 

1888 

3888 

2888 

6888 

1 

1 

-1 

-1 

12738.9891 

3111 

581 

1888 

3888 

4888 

12888 

1  -1 

1 

1 

1 

1 

16431.8888 

3IN 

2811 

488 

1288 

28H 

12888 

1  1 

-1 

-1 

-1 

1 

14889.8888 

3IH 

2888 

4N 

1288 

4888 

6888 

1  1 

-1 

-1 

1 

-1 

13714.9891 

3IN 

2888 

488 

3N8 

2888 

6888 

1  1 

-1 

1 

-1 

-1 

13461.9891 

3111 

2888 

4N 

3888 

4888 

12888 

1  1 

-1 

1 

1 

1 

16767.3688 

3IN 

2888 

18H 

1288 

2888 

6888 

1  1 

1 

-1 

-1 

12396.9891 

3BN 

2888 

1888 

1288 

4888 

128n 

1  1 

1 

-1 

1 

1 

16219.8888 

3IN 

2888 

1888 

3888 

2888 

12888 

1  1 

1 

1 

-1 

1 

16845.8888 

3BN 

2888 

1888 

3888 

4888 

6888 

1  1 

1 

1 

1 

-1 

15432.8182 

TABLE  B.28  (continued) 


NON-CODED  CODED 


RUN  « 

HI 

n2 

n3 

B1 

B2 

SI 

ttl  N2 

R3 

B1 

B2 

SI 

RESPONSE 

33 

3931 

1251 

7N 

21N 

3BN 

9BII  2. 

378  B 

B 

B 

B 

B 

15264. 4S7B 

34 

721 

1251 

711 

21N 

3Bfll 

9BBB  -2. 

378  B 

B 

B 

B 

B 

13445.4951 

33 

2325 

3134 

7N 

2181 

3BII 

9fllB 

B  2.378 

B 

B 

B 

B 

1S474.91BB 

36 

2323 

1 

711 

21BB 

3BBS 

9BH 

B  -2.378 

B 

B 

B 

B 

13566.8682 

37 

2325 

1251 

1413 

21BB 

3BN 

9NB 

B  B 

2.378 

B 

B 

B 

14878.8142 

38 

2323 

1251 

1 

21BB 

3BB8 

9SflB 

B  B 

-2.378 

B 

B 

B 

13895.5682 

39 

2325 

1251 

711 

424B 

3BBB 

9BaB 

8  B 

B 

2.378 

B 

B 

15746.  BB2B 

41 

2325 

1251 

711 

B 

3BBB 

9BSB 

B  B 

B 

-2.378 

B 

B 

128B9.1591 

41 

2325 

1251 

7BI 

21BB 

5378 

9HB 

B  B 

B 

a 

2.378 

B 

15962. 4BBB 

42 

2325 

1231 

7BI 

21BB 

622 

9Bn 

B  B 

B 

B 

-2.378 

B 

12469.3991 

43 

2325 

125B 

711 

21BB 

3nB 

16134 

B  B 

B 

B 

B 

2.378 

16386.9436 

44 

2325 

12SB 

7BI 

21BB 

3BBB 

1866 

B  B 

B 

8 

B 

-2.378 

9953.7664 

45 

2325 

1258 

711 

21flB 

3BBB 

9BBB 

B  B 

B 

a 

B 

8 

14378.5682 

ABBeQdi.x  Qs.  Linear  RBaceisi.on  ResuLts 


This  section  contains  the  results  of  the  stepwise  linear 
regression  run  of  a  sample  input  for  the  three  factor 
problem.  The  experimental  design  which  was  used  to  complete 
this  sample  run  is  the  Box -Draper  bias  minimizing  design  for 
three  factors. 
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PROGRAM  CONTROL  INFORMATION 

/problei  title  is  'saeple  for  three  factor  problee' 

/input  variables  are  4. 

foreat  is  free, 
file  is  bd3a. 

/variable  naees  are  M1,B1,S1,0E, 

Mlsq(l,Blsq(l,Slsqd, 

MIDI, NISI, BlSl. 

add  -  6. 

/transfora  Nlsqd  =  Ml  *  Ml. 

Blsqd  =  Bl-  f  Bl. 

Slsqd  -  SI  t  SI. 

MlBl  =  Ml  i  Bl. 

MlSl  =  Ml  *  SI. 

BlSl  =  Bl  *  SI. 

/regress  depend  is  DE. 

independent  are  Mi,Bl,St, 

Mlsqd, Blsqd, Slsqd, 

MlBl, MlSl, BlSl. 

tol  =  .111. 

/print 

covariance. 

correlation. 


PROBLEH  TITLE  IS 

saiple  input  Tor  3  factor  problea 


NUnSER  OF  VARIABLES  TO  READ  IN .  4 

NUMBER  OF  VARIABLES  ADDED  BY  TRANSFORMATIONS.  .  6 

TOTAL  NWBER  OF  VARIABLES .  IB 

NUMBER  OF  CASES  TO  READ  IN .  TO  END 

CASE  LABELIN6  VARIABLES  . 

MISSIN6  VALUES  CHECKED  BEFORE  OR  AFTER  TRANS.  .  NEITHER 

BLANKS  ARE . MISSINS 

INPUT  FILE.  .......  UNIT  7 . bd3a 

REMIND  INPUT  UNIT  PRIM  TO  READINS.  .  DATA.  .  .  YES 
NUMBER  OF  NORDS  OF  DYNAMIC  STORAGE .  25598 

VARIABLES  TO  BE  USED 

1  Ml  2  B1  3  SI  4  DE  5  Mlsqd 

6  Blsqd  7  Slsqd  8  MlBl  9  MlSl  IB  BlSl 


INPUT  FORMAT  IS 
FREE 

MAXIMUM  LENGTH  DATA  RECORD  IS  SB  CHARACTERS. 


REGRESSION  INTERCEPT . .  .  .NON  ZERO 

HEIGHT  VARIABLE  .  .  . 

PRINT  COVARIANCE  MATRIX  .  YES 

PRINT  CORRELATION  MATRIX .  YES 

PRINT  ANOVA  AT  EACH  STEP .  YES 

PRINT  STEP  OUTPUT  .  YES 

PRINT  REGRESSION  COEFFICIENT  SUMMARY  TABLE.  .  .  YES 

PRINT  PARTIAL  CORRELATION  SUMMARY  TABLE  ....  NO 

PRINT  F-RATIO  SUMMARY  TABLE  .  NO 

PRINT  SUMMARY  TABLE  .  YES 

PRINT  RESIDUALS  AND  DATA .  NO 

PRINT  CORRELATION  OF  REGRESSION  COEFFICIENTS.  .  YES 

PRINT  NORMAL  PROBABILITY  PLOT  .  NO 

PRINT  DETRENDED  NORMAL  PROBABILITY  PLOT  ...  .  NO 

PRINT  DIAGNOSTIC  PLOTS .  NO 

PRINT  PLOTS  FOR  XVAR  AND  YVAR .  NO 

PRINT  PLOTS  AT  EACH  STEP .  NO 

PRINT  PLOTS  AND  DATA .  NO 

PRINT  PLOTS  NITH  STATISTICS  .  MO 

NUMBER  OF  CASES  READ .  15 


C  -  3 


COVARIANCE  HATRII 


HI 

Bl 

SI  Ntsqd  Blsqd 

Slsqd  HlBl 

HlSl 

1 

2 

3  5  6 

7  8 

( 

HI 

1 

.7434 

BI 

2 

.4257e-B8 

.7434 

SI 

3 

.NNem 

.88880+88 

.7434 

nisqd 

S 

-.3692e-a8 

.88880+88 

.K1S0-B8  .  4796 

Blsqd 

6 

-.33620-18 

.88880+88 

.88880+88  -.1838  .4796 

Slsqd 

7 

-.33620-18 

-.42570-88 

.88810+88  -.1838  -.1838 

.4796 

NlBl 

8 

.85150-18 

.25540-87 

.88880+N  .29890-89  .29890-89 

.29890-89  .3319 

NISI 

9 

.42570-88 

.8N80+88 

.12770-87  .  88880+88  .88N0+88 

.88880+88  -.12770-87 

.3319 

BlSl 

11 

.42570-88 

.85150-88 

.12770-87  .88880+88  .88880+88 

.88880+88  .12770-87 

.8n80+8B 

CORRELATION  HATRIX 


HI 

Bl 

SI 

Nisqd 

Blsqd 

Slsqd 

NlBl 

NISI 

BlSl 

1 

2 

3 

5 

6 

7 

8 

9 

18 

HI 

1 

1.8888 

Bl 

2 

.8888 

1.B8N 

SI 

3 

.8888 

.8888 

1.8888 

Nisqd 

5 

.8888 

.8888 

.8888 

i.8n8 

Blsqd 

6 

.8888 

.8888 

.8888 

-.3833 

l.BBH 

Slsqd 

7 

.8888 

.8888 

.8888 

-.3833 

-.3833 

1.8888 

HlBl 

8 

.8888 

.8888 

.8888 

.8888 

.8888 

.8888 

1.8888 

HlSl 

9 

.8888 

.8888 

.8888 

.8888 

.8888 

.8888 

.8888 

1.8888 

BlSl 

18 

.8888 

.8H8 

.8888 

.8888 

.8888 

.8888 

.8888 

.8n8 

1.88N 

STEPPIN6  ALGORITHM . 

NAXinUN  NUMBER  OF  STEPS  . 

•  •  • 

■  •  •  • 

•  •  ■  •  • 

21 

.F 

DEPENDENT  VARIABLE . 

«  •  •  • 

4  DE 

HINIHUH  ACCEPTABLE  F  TO  ENTER  .  .  . 

■  •  ■  ■ 

.  .  4.IIB, 

4.  HI 

NAXINUN  ACCEPTABLE  F  TO  RENOVE.  .  . 

•  •  •  • 

.  .  3.911, 

3.911 

NININUN  ACCEPTABLE  TOLERANCE.  .  .  . 

■  •  «  • 

.  .  .11111 

SUBSCRIPTS  OF  THE  INDEPENDENT  VARIABLES  .  . 

.  .  12 

3 

3  6 

7  8 

9 

11 

STEP  NO.  1 

STD.  ERROR  OF  EST.  1138.3123 

ANALYSIS  OF  VARIANCE 

SUN  OF  SQUARES 

DF 

NEAN  SQUARE 

RESIDUAL  18141236. 

14 

1293733. 

VARIABLES  IN  EQUATION  FOR  DE 

STB. 

ERROR 

STD  RE6 

F 

VARIAn.E  COEFFICIENT  OF  COEFF 

COEFF  T1B.ERANCE  TO  RENOVE  LEVEL 

(Y-INTERCEPT  9161.82713  ) 

.  VARIABLES  NOT  IN  EQUATION 

PARTIAL 

F 

VARIABLE  CORR.  TOLERANCE 

TO  ENTER  LEVEL 

.  Ni  1  .24328  i.nni 

.83 

1 

.  B1  2  .42422  l.BIIII 

2.83 

1 

.  SI  3  .84463  l.HIN 

32.37 

1 

.  Nlsqd  5  .16236  l.lllll 

.13 

1 

.  Blsqd  6  .14787  l.NHI 

.13 

1 

.  Slsqd  7  -.17424  l.BIIII 

.41 

1 

.  NlBl  8  -.11413  l.NHI 

.11 

1 

.  NISI  9  -.16771  I.IINI 

.16 

1 

.  BlSl  11  -.19123  l.NHI 

.11 

1 

VARIABLE  ENTERED  3  SI 


HULTIPLE  R  .8447 
INJLTIPLE  R-SBUARE  .7134 
ADJUSTED  R-SOUARE  .6914 
STD.  ERROR  OF  EST.  632.334S 


ANALYSIS  OF  VARIANCE 

SW  OF  SQUARES 
RE6RESSI0N  12941919. 
RESIDUAL  519B338.S 


DF  HEAR  SQUARE 
1  .12941926^ 
13  399872.2 


F  RATIO 
32.37 


VARIABLES  IN  EQUATION  FOR  DE 


VARIABLE 
(Y-INTERCEPT 
SI  3 


COEFFICIENT 

9161.82715 

1115.15415 


STD.  ERROR  STD  RES  F 

OF  COEFF  COEFF  TOLERANCE  TO  REHOVE  LEVa 

196.1181  .845  l.HNI  32.37  1 


VARIABLES  NOT  IN  EQUATION 


PARTIAL  F 

VARIABLE  CORR.  TOLERANCE  TO  ENTER  LEVEL 


• 

.  N1 

1 

.45821 

l.INN 

3.19 

.  B1 

2 

.79247 

i.mN 

21.26 

.  Nlsqd 

5 

.11651 

l.INN 

.17 

.  Blsqd 

6 

.18792 

l.INN 

.19 

.  Slsqd 

7 

-.32548 

l.INN 

1.42 

.  NlBl 

8 

-.12625 

l.INN 

.11 

.  NISI 

9 

-.12649 

l.INN 

.21 

.  BlSl 

11 

-.17145 

l.INN 

.36 

C  -  6 


.*•«  . 


I 


VARIABLE  ENTERED  2  B1 

RULTIPLE  R  .9452 
HULTIPLE  R-SQUARE  .8934 
ADJUSTED  R-S8UARE  .8756 
STD.  ERROR  OF  EST.  411.4315 


ANALYSIS  OF  VARIANCE 

SUN  SOUARES 
RE6RESSI0N  16206489. 
RESIDUAL  1933766.9 


DF  NEAN  SQUARE  F  RATIO 
2  8113245.  51.28 
12  161147.2 


VARIABLES  IN  EQUATION  FOR  DE 


STD.  ERROR  STD  REG  F 


VARIABLE 

(Y-INTERCEPT 

COEFFICIENT 
9161.82715  ) 

OF  COEFF 

COEFF 

TOLERANCE 

TO  REHOVE 

LEVEL 

B1 

2 

560.87813 

124.4363 

.424 

l.BOOH 

20.26 

1 

SI 

3 

1115.15405 

124.4363 

.845 

l.NOOO 

80.31 

1 

VARIABLES  NOT  IN  EQUATION 


PARTIAL 

F 

VARIABLE 

CORR. 

TOLERANCE 

TO  ENTER 

LEVEL 

HI 

1 

.75126 

1.00000 

14.25 

1 

Nlsqd 

5 

.19101 

l.OHOO 

.42 

1 

Blsqd 

6 

.14416 

1.00000 

.23 

1 

Slsqd 

7 

-.53365 

i.ooon 

4.38 

1 

HlBl 

8 

-.04303 

l.OOON 

.02 

1 

NISI 

9 

-.20740 

I.OOON 

.49 

1 

BlSl 

10 

-.27947 

1.00000 

.93 

1 

C  -  7 


s 


// 


VARIAn.E  ENTERED  9  NISI 


IttILTIPLE  R  .9984 

NULTIPLE  R-SQIMRE  .9968 

ADJUSTED  R-SaUARE  .9943 

STD.  ERROR  OF  EST.  84.7437 

ANALYSIS  OF  VARIANCE 

SUN  OF  SQUARES  DF  NEAN  SQUARE  F  RATIO 
RE6RESSI0N  18882814.  6  3813881.  419.66 

RESIDUAL  37431.914  8  7181.489 

VARIABLES  IN  EQUATION  FOR  DE 

STD.  ERROR  STD  RE6 


VARIABLE 

(Y-INTERCEPT 

COEFFICIENT 
9368.52344  ) 

OF  COEFF 

COEFF 

TIH.ERANCE 

TO  REHOVE 

LEVEL 

HI 

1 

323.83828 

26.2698 

.245 

1.B8888 

151.97 

1 

B1 

2 

368.87813 

26.2698 

.424 

1.S88H 

454.58 

1 

SI 

3 

1113.13483 

26.2698 

.843 

1.B8888 

1882.12 

1 

Slsqd 

7 

-286.38565 

32.7848 

-.174 

1.BI8N 

76.68 

1 

NISI 

9 

-133.79181 

39.3128 

-.868 

l.NBH 

11.38 

1 

BlSl 

IB 

-188.28658 

39.3128 

-.891 

1.NB8B 

21.83 

1 

VARIABLES  NOT  IN  EQUATION 


VARIABLE 

PARTIAL 

CORR. 

TOLERANCE 

F 

TO  ENTER 

LEVa 

Nlsqd 

5 

-.B8S88 

.85387 

.85 

1 

Blsqd 

6 

-.37948 

.85387 

1.18 

1 

NlBl 

8 

-.24967 

l.BHN 

.47 

1 

tH«f  F  LEVELS!  4.NB,  3.988)  OR  TOLERANCE  INSUFFICIENT  FOR  FURTHER  STEPPIN6 


Appendix  Dj.  ResBQQse  Eguati_ons  and  Random  Check  Points 


TABLE  D. 1 

RESPONSE  SURFACE  EQUATIONS  FOR  THREE  FACTORS 

DESIGN  EQUATION 


1) 

BOX  -  BEHNKEN 

Y 

= 

0.4965 (Ml)  +  1.2245<B1)  +‘1.1802 (SI) 

-  .0000375 (SI *S1) 

-  . 0000665 (B1*S1) 

-  648.3319 

2) 

BOX  -  DRAPER 

Y 

= 

0.4798 (Ml)  +  1.2232(B1)  +  1.0847 (SI) 

-  .0000318 (S1*S1) 

-  . 0000667 (B1*S1) 

-  285.6179 

3) 

HYBRID 

Y 

s 

0.5243(M1)  +  1.1947(B1)  +  0.9613(81) 

<311A) 

-  .0000251(81*81) 

-  . 0000602 (B1*S1) 

+  68.4118 

4) 

HYBRID 

Y 

s 

0.5002 (Ml)  +  1.327(B1)  +  1.163 (SI) 

(310+CP) 

-  .0000352(81*81) 

-  . 0000774 (B1*S1) 

-  626. 1627 

5) 

MIN  POINT 

Y 

a 

0.4881 (Ml)  +  1.423 

(Bl)  +  0.969(S1) 

-  .000025(31*81)  - 

. 000088 (B1*S1) 

+  210.5637 


TABLE  D. 1  (continued) 


RESPONSE  SURFACE  EQUATIONS  FOR  THREE  FACTORS 


DESIGN 


EQUATION 


6)  CENTRAL  COMP 
ROTATABLE 


=  0.4848 (Ml)  +  1.153(B1)  +  1.004 (SI) 


-  . 0000276 ( SI »S1)  -  . 0000618 (B1*S1) 


+  124.9285 


7)  KOSHAL 


=  0.5453 (Ml)  +  0.6747(B1)  +  0.6047(81) 


-  .0000142(81*81)  +  2396.592 


D  -  2 


TABLE  D.2 


RESPONSE  SURFACE  EQUATIONS  FOR  FOUR  FACTORS 


desibn 

1)  BOX  -  BEHNKEN 


2)  BOX  -  DRAPER 


3)  HYBRID 
<4160 


4)  HYBRID 

(416A+CP> 


5)  MIN  POINT 


iQUAIIQN 

Y  =  0.3833<M1)  +  0.9461(M2>  +  0.9871(B1> 

+  0.8775 (SI)  -  . 0000229 <S1*S1) 

-  . 0000539 (M2«S1)  -  . 0000512 <B1«S1) 

+  1319.9931 

Y  =  0.3950 (Ml)  +  0.9441 (M2)  +  0.9879(B1) 

+  0.8868 (SI)  -  . 000024 (S1*S1) 

-  . 0000527 <M2*S1)  -  . 0000503 (B1*S1) 

-t-  1332.3008 

Y  =  0.434 (Ml)  +  1.1 172 (M2)  +  1.0793(B1) 

+  0.7926(81)  -  .0000149(81*81) 

-  . 0000677 (M2*S1)  -  .0000645(61*31) 

+  1129.0738 

Y  =  0.4394(M1)  +  1.0915(M2)  +  1.1122(61) 

+  0.823 (SI)  -  .0000174(31*31) 

-  . 0000643 (M2*S1)  -  .0000601(61*31) 

*  1060.9086 

Y  =  0.4088 (Ml)  +  1.0122 (M2)  +  1.055(B1) 

+  0.9707(31)  -  .0000265(31*31) 

-  . 0000582 (M2*S1)  -  .000055(61*31) 

+  810.7629 


D  -  3 


•i,  '  .  •  .'•v- 


TABLE  D.2  (continued) 


RESPONSE  SURFACE  EQUATIONS  FOR  FOUR  FACTORS 

design  eduaiion 

6)  CENTRAL  COMP  Y  =  0.416 (Ml)  +  0.9609 (M2)  +  1.0236(B1) 

ROTATABLE 

+  0.907 (SI)  -  . 0000233 (S1*S1) 

-  . 0000536 (M2«S1)  -  . 0000512 (B1*S1) 

+  939.3 

7)  KOSHAL  Y  =  0.4087 (Ml)  +  0.4833 (M2)  +  0.5322(B1) 

<  '  '•  +  0.2406 (SI)  +5142.2911 


D  -  4 


TABLE  D.3 


RESPONSE  SURFACE  EQUATIONS  FOR  FIVE  FACTORS 

design  equaiion 

1)  BOX  -  BEHNKEN  Y  =  0.6029<M1)  +  0.6795<M2)  0.7294(B1) 

+  0.7689<B2)  +  1.0363 (SI) 

-  .0000295 (SI *S1)  +  843.9903 

2)  BOX  -  DRAPER  Y  =  0.6289 (Ml)  +  0.6776 (M2)  +  0.7556(B1) 

(FULL) 

+  0.7976(B2)  +  1.0942 (SI) 

-  .0000329 (SI *S1)  +  387.9629 

3)  BOX  -  DRAPER  Y  =  0.6037 (Ml)  +  0.6825 (M2)  +  0.7393(B1) 

(HALF-REPLICATE) 

+  0.7846(B2)  +  1.0926(81) 

-  .0000329(81*31)  +  536.6394 

4)  MIN  POINT  Y  =  1.2691 (Ml)  +  2.3605(M2)  +  0.9176(B1) 

+  0. 938(B2)  +  0.6281(31) 

-  .001015  (M1*M2)  -  640.3131 

5)  CENTRAL  COMP  Y  =  0.5011 (Ml)  +  0.6289 (M2)  +  0.7123(B1) 

ROTATABLE 

+  0.7545(B2)  +  0.9511 (SI) 

-  .0000261(31*31)  +  1615.6741 

6)  KOSHAL  Y  =  .7010 (Ml)  +  0.7668 (M2)  +  0.8011 (Bl) 

+  0.8332(B2)  +  0.9541(31) 


000024 (S1*S1)  +  465.7382 


TABLE  D.4  <continued) 


RESPONSE  SURFACE  EOUATIONS  FOR  SIX  FACTORS 


design 

6)  CENTRAL  COMP 
ROTATABLE 


7)  KOSHAL 


iQUAIiON 

Y  =  Q.  5622  (Mi)  0.6062  (M2)  +  0.6728  (M3) 

0.6852(B1)  -t-  1.2372(B2)  1.0462 (Si) 

-  .0000235(81*81)  -  . 0000569 (B2*S1) 
620.122 

Y  =  0.5837 (Ml)  +  0.6582 (M2)  +  0.7360 (M3) 
+  0.7033(B1)  +  0.9630(B2)  +  1.0337(81) 

-  .0000262(81*81)  -  . 0000262 (B2*S1) 


+  956.1621 


TABLE  0.5 


*«* 


RANDOMLY 

SELECTED 

CHECK  POINTS 

FOR  THREE 

HI 

ii 

ii 

1 

1691 

1747 

6014 

2 

1660 

1332 

10359 

3 

2031 

2624 

9151 

4 

2527 

1435 

7752 

5 

2888 

2594 

8166 

6 

1831 

1787 

6474 

7 

1985 

2011 

7846 

8 

2955 

1415 

7076 

9 

2411 

1959 

6152 

10 

1757 

2614 

8822 

11 

2724 

2224 

9792 

12 

2326 

2643 

11871 

13 

2980 

1743 

11841 

14 

1750 

2523 

9492 

15 

1729 

2243 

9347 

16 

2768 

2812 

8627 

17 

2666 

1462 

8818 

18 

2157 

2223 

8082 

19 

2292 

2201 

9033 

20 

1771 

2802 

7756 

21 

1650 

1200 

6000 

♦**  22 


3000 


3000 


12000 


TABLE  D.6 


RANDOMLY  SELECTED 

CHECK 

POINTS  FOR 

FOUR  FACTORS 

Ml 

M2 

il 

il 

1 

2257 

643 

1739 

6585 

2 

1665 

1443 

1400 

7407 

3 

2063 

1087 

1242 

8236 

4 

2645 

1954 

2986 

9142 

5 

1881 

540 

2441 

11017 

6 

2521 

1306 

1352 

8755 

7 

2496 

973 

1283 

7598 

8 

2247 

1943 

1613 

7274 

9 

2779 

553 

2957 

11808 

10 

2550 

1120 

1741 

11147 

11 

1739 

1056 

2371 

11175 

12 

2242 

770 

1539 

11044 

13 

2608 

759 

2499 

8595 

14 

2515 

628 

1278 

10883 

15 

2134 

1876 

1921 

8122 

16 

2401 

1174 

1837 

11724 

17 

1800 

1819 

1373 

11811 

18 

2996 

1972 

1949 

6088 

19 

1707 

1972 

1958 

8375 

20 

1712 

1870 

2027 

11200 

21 

1650 

500 

1200 

6000 

22 


30B0 


2000 


3000 


12000 


TABLE  D.7 


RANDOMLY  SELECTED  CHECK  POINTS  FOR  FIVE  FACTORS 


Mi 

M2 

ii 

12 

SI 

1691 

645 

1747 

3911 

6014 

1660 

535 

1332 

2462 

10359 

152 

822 

208 

598 


11543 


♦*  21 
**  22 


TABLE 

D.8 

RANDOMLY 

SELECTED 

CHECK 

POINTS  FOR 

SIX  FACTORS 

Ml 

M2 

!33 

il 

§2 

§1 

2295 

1249 

561 

2371 

2244 

11147 

1782 

1833 

674 

2733 

2698 

11940 

2316 

768 

845 

2553 

2149 

11165 

11751 


11017 


2 


**  21 
**  22 


11808 
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